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NOTATION 


(Engineering Symbols) 



Symbol 

F0RTRAN1 

name 

Description 

A 


A 

cross-sectional area of local section, m^ (ft^) 

Ap 



total cross-sectional flow area at drive 
fan(s) , m^ (ft^) 

Aflow 


AL 

cross-sectional area of local flow, m^ (ft^) 

Ao 


AO 

cross-sectional flow area of test section at 
upstream end, m^ (ft^) 

Al 


Al 

cross-sectional flow area of section at 
upstream end, m^ (ft^) 

3. 

A2 


A2 

cross-sectional flow area of section at 
downstream end, m^ (ft^) 

A* 


ASTAR 

cross-sectional area for sonic flow at 
specified flow conditions, m^ (ft^) 

AR 


AR 

cross-sectional flow area ratio of upstream 
and downstream ends of section 

3t 


AT 

speed of sound in still gas, computed at total 
(stagnation) conditions, m/sec (ft/sec) 



ASO 

speed of sound in moving flow at upstream end 
of test section, m/sec (ft/sec) 

B 



dummy constraint used in defining the friction 
term of turning vane loss function 



CD 

drag coefficient of flow obstructions: 
drag/qS 

Cy 


CH0RD 

chord of turning vanes, m (ft) , 

D 


D 

cross-sectional diameter of circular duct, 
m (ft) 

Dei 



cross-section diameter at the upstream end of 
an equivalent circular duct with equal area, 
m (ft) 

^Note that in this section 
occurring in F0RTRAN names are 
shovm without slashes. 

, as throughout the report, all letter O's 
shown with slashes, as 0; all number zeros are 


V 


Symbol 

F0RTRAN 

name 

Description 



cross-section diameter at the downstream end of 
an equivalent circular duct with equal area, 
m (ft) 

Dh 

DH 

, , , , ^ 4 X (cross-sectional area) 

hydraulic diameter: . , 

perimeter 

m (ft) 

ER 

ER 

energy ratio: ratio of energy of flow at the 

test section to the output energy of the fans 

f(0) 

FKTVl 

FKTV2 

function defining turning vane loss parameter 
Ktv 

K 

EK 

local total pressure loss coefficient of sec- 
tion: 

q 

^CONTRACTION 

EKCNTR 

local total pressure loss coefficient from con- 
tracting portion of thick-airfoil flow 
straighteners 

^DIFFUSION 

EKD 

local total pressure loss coefficient from 

diffusing portion of multi-loss-type sections 

*^EXP 

EKEXP 

net expansion loss coefficient for diffusers 

^EXPAddltional 

EKADD 

additional diffuser expansion loss coefficient 
due only to more diffusion in one plane than 
the other 

*^EXPBaslc 

EKBASE 

basic diffuser expansion loss factor coeffi- 
cient for three-dimensional diffusion 

^EXPcircular 

EKC 

expansion loss coefficient for conical 
diffusers 

^^Rectangular 

EK2DR 

expansion loss coefficient for a two- 
dimensional, rectangular cross-section 
diffuser 

^EXPsquare 

EKS 

expansion loss coefficient for three- 

dimensional expansion in square cross-section 
diffusers 


EK2DCS 

estimated expansion loss value for a two- 
dimensional diffuser (one with expansion in 
only one plane) with cross-section shape of 
some square/circular hybrid 


vl 


Symbol 

■'EXI-2Dciroul« 


■'EXPjD. 

^*^Average 


F0RTRAN 

name Description 

EK2DC estimated expansion loss value for a hypotheti- 
cal two-dimensional diffuser with circular 
sides: 



EKCSAV estimated expansion loss coefficient for three- 
dimensional, combination circular and square 
cross-section diffuser 


^EXPANSION 


*^FRICTION 
^FRICTION (CONICAL) 


diffuser loss coefficient due to expansion: 



turning vane loss due to friction 

diffuser loss due to friction for the equiva- 
lent conical diffuser 


Kmesh 

%ef. 9 


l^RN 

Dotation 

%v 

Ktvqo 

Kv 

*^VANED DIFFUSER 


K 


o 


EKMESH mesh screen-type loss parameter 

diffuser loss factor presented In reference 9: 

^p/q 

[(AR - 1)/AR]^ 

mesh screen Reynolds number sensitivity factor 

turning vane loss coefficient due to rotation 

EKTV turning vane loss coefficient 

EKTV90 turning vane loss parameter for given vanes at 
a 90® turn 

EKV local total pressure loss coefficient for vaned 

diffusers 

local total pressure loss coefficient for vaned 
diffuser, Ky 

EKO section total pressure loss coefficient 

referred to test section conditions: 

Apx 


vll 


F0RTRAN 

Symbol name 

*^°DRAG 


L 

I 

M 

Mo 

N 

^DRAG 

^INPUT 

riNPUToRAO 

^REQUIRED 


EL 

AMACH 

EMO 

N 


PWRIP 


PWR0P 


^ATM 


PT 


sc 



PT 

PATM 


QO 

R 


Description 

flow obstruction (drag item) total pressure 
loss coefficient referred to test section 
conditions 

centerline length of section, m (ft) 

characteristic dimension on which Reynolds 
number is based 

local Mach number 

Mach number at upstream end of test section 

section assigned sequence number for order of 
occurrence in circuit 

power loss due to drag of flow obstruction, 

W (hp) 

tunnel drive power required to be input to flow 
by the fans, W (hp) 

power input required to overcome drag of flow 
obstruction, W (hp) 

total fan motor output power required to drive 
wind tunnel at specified speed, W (hp) 

local static pressure, N/m^ (Ib/ft^) 

tunnel total (stagnation) pressure, ’ N/m^ 
(lb/ft2) 

atmospheric (barometric) pressure, N/m^ 

(lb/ft2) 

total (stagnation) pressure in the circuit 
settling chamber, N/m^ (Ib/ft^) 

local dynamic pressure: , N/m^ (Ib/ft^) 

PoVq^ , o 

test section dynamic pressure: — ^ — , N/m^ 

(Ib/ft^) 

gas constant, m^/sec^ (ft^/sec^ ®R) 

equivalent radius: /A/tt , m (ft) 


viri 


Symbol 

F0RTRAN 

name 

Description 

RN 

RN 

Reynolds number: 

RNref 

RNREF 

reference Reynolds number at which turning vane 
90° loss parameter, determined 

s 


drag area of flow obstruction (i.e., area for 
which Cj) is determined), m^ (ft^) 

s 


distance along diffuser wall, m (ft) 

®2 


length of diffuser, taken along wall, m (ft) 

T 


tunnel temperature in moving flow, °K (°R) 

Tx 

TT 

tunnel total (stagnation) temperature, °K (°R) 

V 

V 

local flow velocity, m/sec (ft/sec) 



flow velocity at the drive fan(s) , m/sec 
(ft/sec) 

'^SYSTEM 


flow velocity in a multiple-duct section, m/sec 
(ft/sec) 

Vo 

VO 

test section upstream-end flow velocity, m/sec 
(ft/sec) 

X 


location of inflection point in contraction • 
wall (distance from upstream end) , m (ft) 

Y 

G 

specific heat ratio of gas 

A 

RUFNES . 

surface roughness in honeycomb cells, m (ft) 

AER 


difference between estimated and true circuit 
energy ratios; i.e., error in energy ratio 
estimate 

apf 


static pressure rise across the fan(s), N/m^ 
(Ib/ft^) 

APt 


total pressure drop through a section, N/m^ 
(lb/ft2) 

^^TduCT 


total pressure drop through a single duct of a 
multiple-duct section, N/m^ (Ib/ft^) 



total pressure rise across the fan(s) , N/m^ 
(ib/ft2) 


ix 


Symbol 


^PT 


^DUCT 


F0RTRAN 

name Description 

average total pressure rise across a single fan, 
N/m2 (Ib/ft^) 


apt 


SYSTEM 


total pressure drop through a multiple-duct 
section, N/m^ (Ib/ft^) 


Ap.p 


TOTAL 


summation of all total pressure drops through 
the wind tunnel circuit, N/m^ (Ib/ft^) 


^^REQUIRED 


AV 


o 


Ae 


6 


s 


G 


0 

X 


y 

^std 


Pt 


DELP local pressure difference across wind tunnel 

wall, N/m^ (Ib/ft^) 

difference between true and estimated required 
drive power levels for given levels of oper- 
ating velocity and fan efficiency; i.e., 
error in required power estimate, W (hp) 

difference between estimated and true test sec- 
tion operating velocity for given power and 
fan efficiency levels; i.e., error in oper- 
ating velocity estimate, m/sec (ft/sec) 

increment of flow-obstruction downstream influ- 
ence factor greater than unity: e - 1 

(greater than or equal to zero) 

SLR diffuser side length ratio: ratio of change in 

height to change in width from upstream to 
downstream end, or its inverse, whichever is 
less than or equal to unity 

EPS flow-obstruction downstream influence coeffi- 

cient (greater than or equal to unity) 

drive motor electrical efficiency, percent 

ETAFAN fan aerodynamic efficiency, percent 

TH diffuser half-angle, rad 

SLAMDA friction coefficient for smooth pipes 

EMU flow viscosity, N sec/m^ (lb sec/ft^) 

EMUSTD standard-day value of viscosity, N sec/m^ 

(lb sec/ft^) 

EMUT reference viscosity at a known temperature, com- 

puted for still gas (stagnation conditions) , 

N sec/m^ (lb sec/ft^) 


X 


Symbol 

F0RTRAN 

name 

Description 

V 

ENU 

kinetic viscosity of gas, m^/sec (ft^/sec) 

P 

RH0S 

local static density, N sec^/m^ (lb sec^/ft^) 

Pf 

RH0SF 

static density at the fan(s), N sec^/m^ 
(lb sec^/ft^) 

Pt 

RH0T 

density computed for total (stagnation) condi- 
tions, N sec^/m^ (lb sec^/ft^) 

Po 

RH0SO 

static density at upstream end of test section, 
N sec^/m^ (lb sec^/ft^) 

N 

i=l ^ 

SUMEKO 

summation of section total pressure losses 
referenced to test section conditions 

N 


SUMEL 

summation of section centerline lengths, m (ft) 

1=1 

<l> 

PHI 

corner flow turning angle, deg 

20 

TH2 

diffuser equivalent cone angle: 


2 tan-> 


xi 
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Ames Directorate, U.S. Army Air Mobility R&D Laboratory 


SUMMARY 


This report brings together and refines the previously scattered and over~ 
simplified techniques for the aerodynamic design and loss prediction of the 
components of subsonic wind tunnels. General guidelines are given for the 
design of diffusers, contractions, corners, and the Inlets and exits of non- 
return tunnels. A system of equations, reflecting the current technology, has 
been compiled and assembled Into a computer program (a user's manual for this 
program Is Included) for determining the total pressure losses. The formula- 
tion presented Is applicable to compressible flow through most closed- or open- 
throat, single-, double-, or non-retum wind tunnels. A comparison of 
estimated performance with that actually achieved by several existing 
facilities produced generally good agreement. 


INTRODUCTION 


In the past) most of the work on the design of ducts and wind tunnels and 
on the determination of their pressure (and power) losses has been either 
highly specialized, considering only one type of component, or over-simplified, 
covering several types of components but giving only a superficial Idea of 
what parameters are Important. However, for the recent NASA studies directed 
toward new and modified wind tunnel facilities. It has been necessary to do a 
careful job of estimating, easily and quickly, the performance of all circuit 
components. This report brings together, revises, and updates the techniques 
for the aerodynamic design and performance prediction of subsonic wind tunnels. 

The basic procedures and guidelines for the aerodynamic design of 
critical wind tunnel components, as presented In references 1 through 3, have 
been revised and updated, as required. The diffuser and contraction design 
curves developed and suggested herein show the relative design points for 
several existing facilities. Also provided are recommendations derived from 
recent NASA studies on end treatments for non-return wind tunnels. 

The method of loss analysis presented is a synthesis of theoretical and 
empirical techniques. Generally, the algorithms used were those substantiated 
by experimental results. The methods of references 4 through 11 for predic- 
ting component losses have been refined and incorporated. The performance 


calculations, based on user-selected flow conditions at the test section, 
assume that the circuit geometry has been predetermined. 

The comparison of the actual and predicted performance for several 
existing wind tunnel facilities shows generally good agreement. 


CAUTIONARY DESIGN GUIDELINES 


This report presents the means for rapidly estimating the performance of 
a wind tunnel circuit after its geometry has been determined. However, an 
improper design of any of its several components (diffusers, contractions, or 
corners, for example) could result in performance penalties caused by inter- 
action with the flow in other components; such penalties cannot be predicted. 
In addition, improper design could cause poor test-section flow quality which 
would not be indicated by the performance analysis. Therefore, the purposes 
of this section are to point out critical areas of concern in wind tunnel 
design and to attempt to establish proper design criteria. 


Diffusers 

Diffusers, especially those just downstream of high-speed sections, are 
very sensitive to design errors which may cause flow separation. The equiva- 
lent cone angle and area ratio must be properly selected to avoid steady-state 
or intermittent separation of the flow from the diffuser walls. (This separa- 
tion can cause vibration, oscillatory fan loading, oscillations in test sec- 
tion velocity, and higher losses in downstream components.) Generally, proper 
diffuser design requires that, for a given area ratio, the equivalent cone 
angle be constrained below a certain value. ("Equivalent** denotes an imagi- 
nary conical section with length and with inlet and exit areas identical to 
the actual section.) This cone angle should probably be held 0.5° to 1° lower 
for diffusers with sharp corners than for those with a rounded cross section. 

Since the portion of the wind tunnel between the test section and the 
fans is usually the higher-loss segment, it is the most critical in affecting 
circuit performance. Therefore, it was used as a basis for establishing 
recommended design limits as a guide to diffuser selection. It was assumed 
that the fans serve to reenergize the boundary layer of downstream sections 
and that the fans and the upstream and downstream components have no inter- 
action that affects their losses; this may or may not be true (see ref. 12). 

The overall area ratio and cone angle between the test section and fan contrac- 
tion were examined for several wind tunnels. This analysis used the centerline 
lengths of all intervening components, including corners. (The actual effect 
of corners is unknown: they may alter the onset of separation somewhat.) 

Figure 1(a) compares curves for the first appreciable stall for flows with 
thin inlet boundary layers, from references 1 and 2, with the design points of 
selected existing wind tunnels. These curves were used to aid in defining the 
separation trend; good correlation with the symbols is not necessarily 
expected. Figure 1(a) shows that most of these wind tunnels were designed 
beyond (above) the two-dimensional stall curve but below the conical stall 
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curve. (Some of these diffusers are far from conical.) The recommended 
design region, shown in figure 1(b), was positioned with the prior knowledge 
that the NASA-Ames 7- by 10-Foot Wind Tunnel has a partially separated dif- 
fuser just downstream of the test section, and that the NASA-Ames 40- by 
80-Foot Wind Tunnel has some local separation in the corners of the primary 
diffuser. The upper portion of the design region is recommended for diffusers 
with rounded corners, and the lower portion for diffusers with sharp corners. 


Contractions 

Contracting sections are subject to separation in the same manner as 
diffusers; however, the penalties are usually much less severe in the con- 
tracting sections. Separation of the flow can occur if the contraction is too 
short for the amount of area reduction. Figure 2(a) presents the general wall 
shapes suggested in reference 3 and figure 2(b) shows the design boundary for 
these shapes in comparison with the designs of several selected wind tunnel 
facilities. 

From this comparison it is evident that, while some facilities were 
designed more conservatively than others, no design severely exceeds the 
design boundary. Since none of the facilities considered has shown signifi- 
cant contraction-caused flow problems, the design boundary may be considered 
empirically reasonable. Further, reference 13 generally tends to support the 
positioning of the suggested design curve. However, the criteria of refer- 
ence 13 are more conservative due to consideration of viscous effects which 
were neglected in the study of reference 3. 


Corners 

The corner losses in a wind tunnel can be large. To minimize them, turn- 
ing vanes should be used for more efficient turning. Also, as with any other 
high-loss item corners should, where possible, be located in a large-area sec- 
tion where the flow speed is low. Corner vane losses can be minimized in two 
additional ways: (1) by selecting an efficient vane cross-sectional shape and 

adjusting it for proper alignment with the flow, and (2) by choosing the best 
chord-to-gap spacing. 

With reference to item (1) , turning vane shapes can vary from bent plates 
to highly-cambered airfoils. Some sources favor airfoil vanes as being more 
efficient (ref. 4, p. 63) while others claim that thin vanes can have lower 
losses (ref. 5, p. 93). But airfoil vanes with blunted leading edges may be 
more forgiving of misalignment with the flow. The thicker vanes may, 
therefore, hold some advantage. 

When considering item (2) , the best chord-to-gap ratio depends on the 
vane type. For thick vanes, a ratio of about 2.5-to-l is recommended (ref. 4, 
p. 62) and for thin vanes a ratio of about 4-to-l is suggested (ref. 5, p. 92). 
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Non-Return Wind Tunnels 


Non-return wind tunnels have presented some interesting problems in 
tunnel design. This type of wind tunnel has the advantages of less structure 
(and therefore lower construction costs) and of no exhaust-gas-purging or air- 
exchange requirement. Careful design can make the non-return circuit oper- 
ating power competitive with that of closed-return wind tunnels (the corner 
losses can be traded for inlet and exit losses). However, an area of concern 
for the non-return tunnel is its potential sensitivity to external winds which 
could affect both the required power and the test section flow quality. 

A recent series of NASA studies, which dealt with wind s.nsitivity prob- 
lems, showed that a non-return wind tunnel should have three matures ; (1) a 

vertical exit system, (2) a horizontal inlet, and (3) an enclosed area of pro- 
tection, with a solid roof, at the inlet. References 14 and 15 detail the 
development work for the end treatment considered in those studies. 

Reference 16 describes an inlet geometry that was developed to reduce the 
effects of wind. (This reference also presents a set of test-section flow- 
quality requirements by which the characteristics of any inlet treatment may 
be evaluated.) Although the end treatment designs shown in references 14 
through 16 could be revised or refined for additional wind protection, any 
additional inlet treatment would increase the structural cost and could 
increase the power requirement. 


PERFORMANCE ESTIMATION 


Although the performance analysis presented in this report was systema- 
tized and automated for rapid calculation of numerous cases or iterations (by 
the computer program described in the following sections) , the equations pre- 
sented are equally amenable to manual calculation methods. 


General Approach 

The equations were derived in forms that use the most common and conveni- 
ent defining parameters. The equations are listed and explained below and may 
be used for component after component, each in turn. 

The total pressure losses (proportional to power losses) of each compo- 
nent are calculated and summed to give the total circuit loss and operating 
power required. The computation technique is applicable to either closed- or 
non-return circuit types made up of any combination of standard wind tunnel 
components in any order. The flow conditions in the test section (velocity, 
and stagnation temperature and pressure) and the external atmospheric pressure 
are variable as required. 
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Problem Restrictions 


Three restrictions were found to be necessary in order to allow rapid 
solution of most cases with a minimum amount of effort. First, the cross- 
sectional geometries were limited to the most common shapes: circular, rec- 

tangular, and flat-oval (semi-circular side walls with flat floor and ceiling) . 
Second, air exchangers were omitted from this analysis due to lack of uniform- 
ity of configuration and a lack of definition as to the proper method of com- 
puting the losses. Finally, the drive system was assumed to be located in one 
or more parallel, annular ducts. 

o 

'i Computation Formulas 

The equations used in this performance analysis were synthesized from 
various sources. Some were used in their original (source) form and others 
were modified to make them more convenient for use in this analysis. The equa- 
tions used are presented below. 

Ftcw-state parconetet‘8- The basic flow-state parameters were determined 
from input information about the reference control station and the test sec- 
tion. These parameters were derived from standard relationships for 
compressible flow. 



(ref. 17, p. 8) 



(ref. 17, p. 51) 


(ref. 18, p. 19) 


(ref. 18, p. 4) 


(ref. 18, p. 4) 


(appendix A) 


Local conditions- The local flow conditions were determined for each end 
of each section. 
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1. Mach number: The local Mach number was found from a Newton ’s-met hod 

solution of the relationship 


- 




2(y-1) 

Y+1 


Y - 1 


= 0 


(appendix A) 


2. Reynolds number; The Reynolds number based on the characteristic 
length a, usually the local hydraulic diameter, was determined from 


RN = 


PqVq ^ /A, 

yrp 




0.76 


(appendix A) 


3. Friction coefficient: A Newton’ s-method solution was used to deter- 

mine the friction coefficient for smooth walls from the expression 


[log^^(XRN^) - 0.8] ^ - X = 0 (appendix A) 

Section pressure losses- The loss in total pressure caused by each sec- 
tion was calculated in a form non-dimensionalized by local dynamic pressure: 

K = Apj/q. (In this study the smallest-area end of each section was used as 
the local reference position.) The individual losses were based on the nature 
of the section, local flow conditions, and input geometry and parameter infor- 
mation. The most appropriate loss forms for typical wind tunnel sections are 
catalogued on the following pages. The nonstandard formulas, those which are 
not directly attributable to the literature, are developed in appendix A. The 
precise equations, which were developed from various curve-fitting and inter- 
polation techniques based on the plots presented in certain figures, are given 
in appendix B. 

1. Constant-area ducts: For closed, constant-area sections the pressure 

loss due to friction is given by 

K = ^ (ref. 7, p. 53) 


2. Open-throat duct: The losses from an open-throat test section may be 

found from the expression 

K = 0.0845 ~ - 0.0053 (ref. 7, p. 150) 

°h \®h/ 

3. Contractions: In contracting sections, where the major part of the 

losses is due to friction, the local loss may be approximated as 


K = 0.32 


All 

Dh 


(ref. 6, p. 528) 
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4. Corners with no net area change (^constant area”): A duct can change 

direction with or without the aid of flow guide vanes. For a constant-area 
turn employing turning vanes for efficiency, with a "normal” number of vanes 
(ref. 7, p. 241), and with chord-to-gap ratios between 2-to-l and 4-to-l, the 
losses resulting from friction and rotation caused by the vanes are 

The Reynolds number used for the turning vane loss should be based on vane 
chord. The turning vane loss parameter Kr^y is plotted as a function of turn- 
ing angle in figure 3(a), with the assumption that Kr^y = 0.15 is a reasonable 
value for a 90^ corner. Corners without turning vanes are less efficient and 
the loss function may be approximated by a sixth-order polynomial as shown in 
figure 3(b): 


K = 4. 313761x10"^ - 6.021515x10 ** (|) + 1.693778x10 (})^ - 2.755078x10 ® (l)^ 

+2.323170x10“'^ (l)*^ - 3.775568x10“^ + 1. 796817x10“! ^ (f)® 


(appendix B) 

This function assumes a loss value of about K = 1.8 for a 90° turn. The 
foregoing losses are those associated with the turning of the flow only. The 
losses for a corner system (with or without vanes) , with the walls of the duct 
to be considered as well, requires an additional term for the frictional loss 
of the constant-area duct based on the centerline length. 

5. Corners (diffusing): Corners with diffusion may well employ longer 

vanes in order to improve the efficiency of the diffusion process. For this 
reason they were treated as vaned diffusers with the addition of the rota- 
tional loss term of the turning vane function: 


K = |0.3+ [0.006(26 - 21.5°)u(26 - 21.5°)]}^ -^ - + -| 

(appendix A) 

where u(26 - 21.5°) is the unit step function and the turning vane loss 
parameter is defined as for a constant-area corner. This loss function 
includes the effects of friction. 

6. Diffusers: Diffusion produces both expansion and friction losses in 

the duct given by 

" ■ ["exp + (s sL 6 (appendlK A) 


7 


where the expansion parameter values, K^xp* plotted against equivalent cone 
angle in figure 4 and the technique used for estimating the KgxP values is 
described in appendix A. 

(It should be noted that there are more sophisticated techniques for 
estimating diffuser performance than the one presented here. However, they 
require boundary- layer calculations; for example, see reference 19. Experience 
with both the simple technique described herein and more complex techniques 
indicates that the two produce comparable results. Generally, little is gained 
by the significant additional effort required to use the more complex 
approaches.) 

7. Exit: The total pressure loss at the exit of a non-return wind tunnel, 

or of any expelled flow, is due to the loss of the kinetic energy of the exit- 
ing flow. This is given by 


+ - 1 

8. Fan (power) section: Fan drive sections are commonly made up of con- 

tractions, constant-area annular ducts, and diffusers. Analysis should be 
handled by dividing the fan section into these three component parts. 

9. Flow straighteners - honeycomb (thin walls): The loss through thin 

f low-straightener or honeycomb systems may be expressed as 

where the hydraulic diameter is that of the honeycomb cell. The friction 
coefficient is determined from a Reynolds number based on the surface roughness 
of the honeycomb: 

for RN^ 275, X = 0 . 375 RN"° ’ ^ ^ 

and for RN > 275, X = 0.214 

10. Flow straighteners - airfoil members (thick walls): Flow through 

adjacent airfoils will first contract and then diffuse. It was assumed that 
the point of minimum distance between parallel members would be at 30 percent 
of the straightener length back from the leading edge. The forward 30 percent 
was treated as a contraction and the aft 70 percent as a vaned diffuser. Thus, 

K = 0.096 ^ + {0.3+ [0.006(20 - 21.5°)u(20 - 21 . 5°)][ ( 

(appendix A) 



(appendix A) 
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where the hydraulic diameter is that of each cell of the flow straightener , 
the friction coefficient is determined from a Reynolds number based on that 
hydraulic diameter, the area ratio and equivalent cone angle are based on the 
exit and minimum flow areas, and u(20 - 21.5°) is the unit step function. 

11. Internal flow obstruction - drag item; The loss due to the drag of 
internal structure such as struts or models has the form 


K = Cd 


S 

^FLOW 


(appendix A) 


12, Perforated plate: Perforated plate with sharp-edged orifices, used 

as protection screen or as screen around the inlet of a non-return tunnel, 
produces losses given by 


K 


f F/l A ^FLOW\ / ^LOW^ l A 
" 11/ 2 V" ■ A ; ■ A yjAFLOwJ 


(ref. 7, p. 321) 


A / V ^ 

13. Mesh screen: The losses produced by a mesh screen may be expressed 


as 


K = KrnKmesh 




- ■) 


(ref. 7, p. 308) 


where the Reynolds number influence factor, is plotted against Reynolds 

number (based on mesh diameter) in figure 5, and the mesh constant, is 

1.3 for average circular metal wire, 1.0 for new metal wire, and 2.1 for silk 
thread. 


14. Sudden expansion: For a sudden expansion with ducting downstream 

(to allow reattachment of the flow and maximize the pressure recovery) the 
loss is 


K = (ref. 7, p. 128) 

15. Vaned diffusers: The pressure loss of a vaned diffuser, one in 

which splitter vanes are used to improve the performance of a short diffuser 
by decreasing the effective equivalent cone angle of each chamber, may be 
determined from 

2 

K = ■j0.3+ [0.006(20 - 21.5°)u(20 - 21.5°) 

(appendix A) 


where u(26 - 21.5°) is the unit step function. (See fig. 6.) 

16. Fixed, known loss: For a fixed loss item where the pressure loss 

value is known, that value may be used directly by definition: 
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K = 


Apij 

q 


17. Multiple ducts: In a system of multiple ducts, where the local flow 

passes through two or more separate, identical passages at the same time, the 
losses have the same value as those for the same type of single duct. Some of 
the pertinent parameters, such as hydraulic diameter and equivalent cone angle, 
should be based on the geometry of one of the individual ducts. The loss for 
the system of ducts may then be determined from the loss for a single duct: 


^^TsySTEM ^^TpucT 


18. Loss value transferred to reference location: Each local loss 

parameter is calculated base! on local conditions at the smallest~area end of 
each section and may then be referenced to the test section conditions by the 
formula 


“o / 1 +(^h 2 ) 


(appendix A) 


19. Overall and summary performance: The energy ratio of the wind 

tunnel under consideration is given by 


ER = (ref. 4, p. 69) 

Z Ko. 

i=l 


The pressure difference across the wind tunnel walls, determining the minimum 
required structural strength for each section, is given by 




- Plj. 


ATM 


PT 


SC 


- 6° S ■^^)l 




Y 


(appendix A) 


The power required to be input into the flow in order to drive the flow 
through the wind tunnel at a specified test section speed is expressed as 

(l) Koi)po"*o''o* 

'■input = 2^^ (appendix A) 

The actual drive power required is dependent on the efficiency of the 
fan/motor system: 
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COMPUTER PROGRAM DESCRIPTION 


The computer program was written in F0RTRAN IV language. It consists of 
a main program which calls five subroutines and/or six library routines, as 
required. Two of the subroutines are optional and may be abbreviated and 
simulated in order to save execution time and/or memory storage space. 


Method of Solution 

The general technique used is outlined in the computer program functional 
flow chart in figure 7. The program was developed in six functional units; a 
main program and five specialized subroutines. The main program retains 
general control over the computational flow and calls the subprograms as 
required. 

In the main portion (designated PERF0RM) , at first entry into the program, 
various section-shape geometry relationships and certain semi-empirical dif- 
fuser, turning vane, and honeycomb loss functions are defined. The case title 
card is read and checked for validity by specified code. The tunnel master 
data control card is then read, checked for validity, and checked for content 
of pertinent data by the data-checking subroutine. If any errors are found in 
either of these two preliminary cards, error messages will be printed. 

Although detected errors will not abort the computer run (unless a card of 
improper format is encountered where not expected) , the case under current 
consideration will not be computed - only the checking of input errors will 
then be performed on each section card. Prior to reading the section cards, 
the units of measure (International System or U.S. Customary System) to be 
used for the particular case are read. These units of measure are used as the 
basis for the development of the appropriate flow parameters and test section 
conditions . 

The section cards are read and operated on one at a time. They are 
checked for validity and input errors by the data-checking subroutine (called 
DATACK) and the input information, if sufficiently complete, is then used in 
the computation of the section upstream- and downstream-end geometries. 
Adjustments to these geometry calculations are made for any multiple-ducted 
sections. For diffusing sections where the expansion loss parameter was not 
input by the user, that parameter is generated from predefined functions. 
Branching of the computational flow then transfers control to the appropriate 
block of instructions for the remainder of the calculations which are peculiar 
to the particular section under consideration. 

After all section cards have been read in and operated on, each in turn, 
a case termination card is encountered. The termination card specifies the 
optional summary operations to be performed. The encounter of this card 
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signals the end of a case and triggers the final calculations. The codes 
contained on this card determine the printing of velocity-optimizing and cir- 
cuit summary information, the plotting of the summary information, and the 
return to the beginning for another case. 

The data-checking subroutine evaluates the master and section input cards 
for completeness of data (based on the requirements for the type of section) . 
Then, if any error was detected during computation of a case or if the 
appropriate termination code was specified by the user, the complete set of 
input data is tabulated. Messages about errors, omissions, or superfluous 
information are Included. 

The subroutine SPEED computes the local Mach number based on local cross- 
sectional area and determines the local flow velocity. 

The subroutine FRICTN calculates the local Reynolds number, visually based 
on the local duct hydraulic diameter, and the local friction coefficient for 
smooth pipes. 

The subroutine 0UTPUT accepts the calculated section parameters along 
with the section type codes describing the types of information to be output 
and prints the section information according to the appropriate format. 

The subroutine PL0TIT plots the summary information (cummulative total 
pressure losses and/or wall pressure differential) versus circuit centerline 
length if requested. The plot is scaled for centimeter or inch plot paper, 
determined by whether International or U.S. Customary units are used for 
computation. 

The program is terminated after the last operations on a case for which a 
no-return instruction on the termination card was given by the user. 


Computing Equipment Required 

Hardware and machine oomponents- Although this program was written for 
use on an IBM 360/67 with TSS Monitor, batch mode, an attempt was made to keep 
it compatible with any system that uses F0RTRAN IV . No magnetic tapes were 
used. In this version, input is made by cards and the data to be plotted are 
stored on a disc file for plotting at a later time in an off-line mode. How- 
ever, it is possible to use a typewriter-type terminal for conversational or 
real-time computation, typing the data by card-image format, and plotting 
immediately after the computation has been completed for a case. 

The total core required for compilation on the IBM 360/67 was approxi- 
mately 82 000 (decimal) bytes. If necessary this figure can be reduced by 
eliminating two subroutines, DATACK and PL0TIT. The sizes of the main program 
and of each subroutine were approximately as follows: 
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PERF0RM 

38 

800 

DATACK 

25 

700 

SPEED 


800 

FRICTN 


900 

0UTPUT 

12 

900 

PL0TIT 

2 

900 


The program was executed on an IBM 360/67, writing plot data on a disc, 
logical unit 10. Later the data file was accessed from a 14.8 character-per- 
second binary-coded-decimal terminal and plotted on a Zeta plotter with 
0.005- in. step increments. The plot page size was programmed not to exceed 
25 by 38 cm (or 10 by 15 in.). 

Software- This program was written for use on any computer with suffici- 
ent core and with a standard F0RTRAN IV compiler. 

The Zeta plotter routines, with minor exceptions, are compatible with the 
Calcomp routines. The subroutines AXIS, FACT0R, LINE, PL0T, SCALE and SYMB0L 
are alike in both Calcomp and Zeta plotting. 

CALL AXIS - draws the axis line and annotates the divisions at every two 

centimeters or each inch (depending on the units of measure specified) • 

CALL FACT0R - enables the user to produce normal size drawings with plotters 
which have either 0.01- or 0.005-in. increment size. The variable FACT 
must be set to 1.0 for 0.01-in. increments and to 2.0 for 0.005-in. 
increment plotters. 

CALL LINE - plots centered squares connected by straight lines through the 
coordinate pairs of data values. 

CALL PL0T - is used to establish a new point of origin for the pen and paper 
movements. Before plotting commences, the pen must be positioned where 
desired along the X-axis. The program will position it along the 
Y-axis. The plot-page size is defined by the values of YLEN and XLEN 
which are equated to 25 and 38 cm or 10 and 15 in., as required. 

CALL PL0TF - is an alternate plotting initialization routine which is avail- 
able in the Zeta but not Calcomp plot package. It is used in place of 
PL0TS whenever deferred plotting is desired. The first argument in the 
call statement indicates the speed of the terminal with which the plotter 
is interfaced. The second argument is the logical-device number of the 
plot file. 

CALL SCALE - examines the data and determines the proper scaling for the given 
dimensions of the plotter paper, 25 by 38 cm or 10 by 15 in. 

CALL SM0DE - is available only in the Zeta plot package. It permits the user 
to choose from extensive capabilities which affect several of the plotter 
routines. In this program the options have been set equal to the usage 
found in the Calcomp routines, and therefore, if Zeta plotter routines 
are not available, the call to SM0DE should be eliminated. 
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CALL SYMB0L - prints the input case title at the top left of each plot page as 
it appears in columns 2 through 80 of the title card. For reference pur- 
poses, it also draws a small plus sign at the origin of the plot. 

The library routines used are standard F0RTRAN routines: 

ABS - Absolute value 

AL0G1O - Common logarithm, base ten 

ATAN - Arctangent (result in radians) 

EXP - Exponential of the natural number e 

IFIX - Convert from real number to integer 

SIN - Trigonometric sine (argument in radians) 

SQRT - Square root 


Programming Techniques 

It was intended that this program be usable on as many different computer 
systems as possible. Therefore, in order to make them applicable to some 
machines, certain statements were forced into particular forms which would be 
less efficient on other systems (e.g., Hollerith instead of literals in format 
statements) . 

C0MM0N and DATA statements were used as much as possible to simplify the 
definition of parameter values. In the main program, arithmetic statement 
functions were used for three purposes: (1) for the definition of section 

hydraulic diameter, area, and equivalent cone angle geometry functions; (2) for 
the conversion function from local to reference-section pressure losses; and 
(3) for the definition of the least-squares-polynomial-curve-f it functions. 

The last group of functions includes: (1) the corner turning-loss parameters 

as functions of turning angle (see fig. 3); (2) the diffuser expansion loss 
parameters for the different cross-section shapes as functions of equivalent 
cone angle (see fig. 4); and (3) the mesh screen Reynolds number sensitivity 
factor as a function of mesh-diameter-based Reynolds number (see fig. 5). 

Certain functions not easily solvable in closed form were solved itera- 
tively (some by Newton’s method) to 0.01 percent accuracy. These functions 
include test section Mach number, local section Mach number, and local section 
friction coefficient. 

Numeric codes were used for specifying such things as section type, sec- 
tion end-shape types, and system of computational units; for decisions on 
requirements for inputs to each section type; and for case-termination pro- 
cedures and outputs desired. The various important input codes are listed in 
tables 1 and 4. All sections of the multiple-ducted type were assigned high 
code numbers for simplicity in selecting them for special handling. The 
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various section types were grouped in code decades for reasonable association 
of section code and component function. Where possible, the second digit of 
the code (if that second digit is not zero) reflects the basic characteristic 
of the section: constant area, contracting, or diffusing. 

The information input fields on the master data and section cards were 
arranged in three basic groups: (1) qualitative information (type and shape); 

(2) quantitative geometry information (number of ducts, cross-sectional end 
dimensions, and length); and (3) loss-related parameters. The case termination 
card employs the same format as the section cards so that it may be encoun- 
tered at random intervals without causing a program crash. For the tabulation 
of the input data (for error-location and record-keeping purposes), object- 
time formatting was used to compile the combination input and annotation data 
set for a convenient output. 

Much of the output of the program was set up on a demand (i.e., optional) 
basis. A section-by-section performance analysis is automatically provided. 

A brief summary of the variation of selected parameters through the circuit, 
and plots of those parameters, may be selected if desired. An annotated list- 
ing of the inputs may be requested or, if errors are detected, the listing is 
internally forced in order to provide a simple means of error-detection and 
correction and/or simplified record-keeping of case data. 


Source Code 

A source code listing of the performance estimation program is provided 
in appendix C along with the associated notation definitions. The source code 
includes the use of comment cards throughout the program for identification of 
the operations carried out by each set of instructions. 


Operating Instructions 

The basic source program deck arrangement is shown in figure 8. 

Input- Sample coding forms for the four types of input cards required are 
presented in figure 9. The special symbols required in the first columns of 
the title and master data cards are included. 

1. Title card: For the title labeling card, with the exception of the 

first column which must contain an asterisk (*) , the entire card may be used 
as desired. This title was programmed to appear at the top of each page of 
the case to which the title refers, including the plots. Only one title card 
per case may be used. 

2. Master card: The tunnel master control data card provides sufficient 

information for defining conditions in the test section (which is the refer- 
ence section for all calculations) and conditions of the surrounding external 
atmosphere. Table 2 details the inputs included on the master data card. The 
first column must contain a minus sign in order to identify the card as a 
valid master card. The remainder of the inputs should be positive, with 
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columns 2 through 6 containing five fields of integers only (no decimal points). 
Columns 7 through 10 were not used on this card and should be left blank. 
Columns 11 through 50 should contain floating-point numbers. These columns 
were divided into eight parameters of five columns each, Including decimal 
point. 


3. Section cards: The inditridual section information cards were based 

on the same format as the master card, except that the section cards require 
no special identifying code. Table 3 details the inputs contained on the sec- 
tion cards. The first six columns, containing four data fields, require inte- 
ger inputs. The remaining 74 columns were divided into two real number fields 
of two columns each (with the assumed default decimal points to the right of 
the second columns), and 14 real number fields of five columns each (with the 
assumed default decimal point between the third and fourth columns of each 
field). 

Although the input parameter requirements vary from section to section, 
certain requirements are basic to all sections. These items include: (1) the 

section type code, (2) the section end shape codes, and (3) the section dimen- 
sions (end height(s) and width(s) and/or diameter(s) and usually length). A 
detailed list of the additional, specialized requirements for each section is 
presented in table 4. 

Although not mandatory in order to obtain a correct total power estima- 
tion, it is advisable to input the section data cards in the actual section 
order so that the summary calculations and plots have relevance to the actual 
circuit . 


4. Termination card: The case termination card, which signals the end 

of the section inputs for a particular case, is identified by the constraint 
of blanks in card columns 3 and 4. The numbers contained on this card are 
used strictly as task codes; table 5 shows the details of these codes. In the 
event of a request for plotted information, the code determines the type(s) of 
information to be presented. For all other tasks the codes dictate a simple 
yes/no decision. 

As many cases as desired may be input in a single job submission. The 
same system of units need not be used in all cases. Any parameters may be 
changed as desired from case-to-case since there are no forced carry-overs 
(except the specific heat ratio, y> which is fixed at the time of program 
compilation) . 

Output- Based on the foregoing input information the results may be 
calculated and tabulated in five different types of information groups. 

1. Section performance analysis: The section performance tabulation 

fully describes the performance-related parameters of the wind tunnel circuit. 
Atmospheric and test section flow reference conditions are stated at the top 
of the first page. The various parameters are tabulated for each section in 
the order of computation with the upstream end information on the first and 
the downstream information on the second of the two lines for each section. 

The section sequence number and type (a translation of the code) and the end 
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shapes are given first. The geometry and local velocity information are 
presented next, followed by the section length and calculated total pressure 
loss values. 

2. Overall performance: If no input errors are encountered during the 

analysis of a case, overall performance values are presented at the end of the 
section performance tabulation. This includes the total circuit length, the 
total pressure losses and energy ratio for the circuit, and the total 
operating power required. 

3. Summary characteristics: If requested on the termination card and 

barring any errors, a summary of the circuit characteristics is tabulated on a 
separate page. This tabulation includes section sequence numbers, Mach num- 
bers, cumulative pressure losses, and local wall pressures, all as functions 
of distance through the circuit. 

4. Plots: Under proper condition codes, the cumulative pressure losses 

and/or the wall pressure differentials will be plotted as functions of dis- 
tance through the circuit (centerline length) . The straight lines that appear 
on the plots connecting the points are for reference only and do not represent 
the actual distribution in a component. 

5. Input data tabulation: Finally, if an input error was encountered 

during the analysis of the circuit, or if such information was requested by 
the user, the input cards are tabulated with annotations regarding missing or 
superfluous Inputs. A careful look at this section should allow the user to 
discover why a given set of input data did not produce the expected type of 
results. 

All of the foregoing types of output are shown for the test case. 

Computer system restrictions- Certain restrictions and/or assumptions had 
to be imposed on the computer system and its methods and abilities in order to 
perform the performance analysis within reasonable time, effort, and money 
constraints. 

1- Hardware: This analysis was programmed for a moderate-sized system 

with common components. No special hardware is required with the exception of 
a plotter if the plotting option is used. The output printing device is 
assumed to have available a minimum capacity of 120 characters per line, but 
the number of lines per page may be set by means of the LINEMX parameter in 
the main program. (Barring any special requirements, 45 lines for an 8.5-in. 
page or 60 lines for an 11- in. page are recommended.) 

2. Software: Certain software restrictions were imposed simply as a 

starting point to the problem solution. The input card formats were fixed as 
shown in figure 9. The specific heat ratio (y) and the number of lines per 
output page were fixed for each compilation of the source deck, although 
changes can be made by altering the values of G and/or LINEMX, respectively, 
near the beginning of the main program. 
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For reasons of possible memory limitations on smaller systems, the number 
of wind tunnel components in each circuit case was limited to 30 sections. 

This limit may be changed by assigning a new value to LMTSEC in the main pro- 
gram and by re-dimensioning the following variables as denoted by "XX”: in 

the main program (PERF0RM) , DELP(XX+2), SEKO(XX), SEL(XX), SMACH(XX), 
SSUMEL(XX+2) and SSUMK0(XX+2) ; in the data-checking subroutine (DATACK) , 
ENDATA(XX,20), NCHECK(XX, 20) and NDATA(XX,4); and in the plotting subroutine 
(PL0TIT), DELP(XX+2), SSUMEL(XX+2) and SSUMK0(XX+2) . If memory limitations 
are a severe problem and/or if computer-controlled plotting facilities are not 
available to the user, the data-checking and/or plotting subroutines may be 
"removed" by inserting dummy, one-card subroutines with the same arguments 
which would have no effect on the calculations. This would decrease the util- 
ity and power of the program, but would retain the basic performance estima- 
tion capabilities without crippling them altogether. 

The plotting routines were written according to the requirements for a 
plotter with 0.005- in. increments. 

Optional inputs- Certain of the parameter inputs are designated as 
optional and have built-in assumed default values in the event that the user 
knows no better values than the ones provided in the sources referenced herein. 
These optional parameters are shown in tables 2 through 4. 

On the master card (see table 2), the units of measure should be speci- 
fied and an error message will be given if they are specified erroneously 
(other than as type 1 or 2) . However, the units code will default to 1 (the 
International System) and case execution will continue. The test section and 
atmospheric total pressures will default to one atmosphere if not specified. 

On the section cards (see table 3) , the number of items in the duct will 
default to unity if not specified. The expansion loss parameter for diffusers 
defaults to a value based on figure 4. (It is computed by determining the 
shape of each end, the extent to which the diffuser is two-dimensional in 
nature (i.e., changing cross-sectional size in height or width only), and the 
equivalent cone angle, and then interpolating between the curves of figure 4. 
See appendix A.) The mesh screen loss constant defaults to 1.3, the value for 
an average-condition metal mesh screen (ref. 6, p. 527), and the reference 
Reynolds number for turning vanes defaults to 0.5 million (ref. 6, p. 527). 

The surface roughness for honeycombs defaults to the appropriate equivalent of 
0.00001 m, the value for new, commercially smooth, non-steel pipe (ref. 7, 
p. 62). The factor for the additional influence of a blockage on downstream 
sections (Ac) defaults to zero . 

Diagnostic messages- There are a limited number of error diagnostic 
messages which were built in to handle many, but not all, of the potential 
user errors. The causes and appropriate corrections of these errors should be 
evident in each message. 

1. Title card: If a card is in the position of a title card and does 
not begin with an asterisk as required, the following message will appear: 
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TITLE (’...(invalid title)...') IS INCORRECT OR IMPROPER AS IT EXISTS. 
THE FIRST CARD COLUMN MUST CONTAIN AN ASTERISK (*) TO BE IDENTIFIED 
AS A VALID TITLE CARD. 


2. Master card: An invalid master card is denoted by: 

MASTER CONTROL DATA ('...(card image)...') IS INCORRECT OR IMPROPER 
AS IT EXISTS. THE FIRST TWO CARD COLUMNS MUST CONTAIN A NEGATIVE 
NUMBER (-1 TO -9) TO BE IDENTIFIED AS A VALID MASTER CARD. THIS 
CASE WILL BE SKIPPED. 


A general omission from the master card of required information produces: 


CRITICAL OMISSION(S) IN TUNNEL MASTER CONTROL DATA PREVENT EXECUTION 
OF THIS CASE. ANY SUCCEEDING CASES WILL NOT BE AFFECTED. 

Two master cards, back-to-back, for a given case are identified by: 

MORE THAN ONE MASTER CONTROL CARD EXISTS FOR THIS CASE OR INPUT CARDS 
ARE OUT OF ORDER. CHECK DECK SET-UP. THE LAST MASTER CARD ENCOUNTERED 
WILL BE ASSUMED AS THE CORRECT MASTER CARD FOR THE SECTION CARDS WHICH 
FOLLOW. 

Encountering a master card where not expected (generally indicating missing 
case termination and title cards) causes this message: 

MASTER CONTROL CARD HAS BEEN ENCOUNTERED BEFORE CASE TERMINATION AND 
TITLE CARDS. CHECK DECK SET-UP. ERROR-MESSAGE TITLE WILL BE GENER- 
ATED AND SUMMARY OUTPUT, NO-PLOT, INPUT DATA TABULATION, AND NEXT-CASE 
RETURN TERMINATION PARAMETERS WILL BE ASSUMED. 


If an invalid test section upstream end shape geometry is specified, one which 
the program cannot handle, an error results: 

**ERROR “ INVALID TEST SECTION UPSTREAM END SHAPE CODE WAS SPECIFIED 
AS (code used) (SHOULD BE 1, 2 OR 3) . THIS CASE CANNOT BE EXECUTED. 

If an invalid units code is specified the message is: 

THE UNITS OF MEASURE CODE IS IMPROPERLY SPECIFIED AS (code used), 

(SHOULD BE 1 OR 2) . CHECK MASTER CARD (COLUMN 4) . SEE THE DATA 
TABULATION AT THE END OF THIS CASE. THE INTERNATIONAL SYSTEM OF 
UNITS WILL BE ASSUMED FOR THIS CASE. 

If the termination code requests power-matching but the input power value is 
such that the calculation would be meaningless, a diagnostic of the following 
form is printed: 

** ALTHOUGH VELOCITY-OPTIMIZING WAS REQUESTED BY TERMINATION CODE, THE 
INPUT POWER VALUE IS ILLEGAL (LESS THAN OR EQUAL TO ZERO). THEREFORE, 

NO VELOCITY-OPTIMIZING IS POSSIBLE. RECHECK INPUT VALUE ON MASTER 
DATA CARD. 
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3. Section card: A general omission of required data from a section 

card will cause this message: 

**ERROR -- CRITICAL OMISSION (S) IN SECTION INPUT DATA. SEE DATA 
TABULATION AT END OF OUTPUT FOR THIS CASE.* ** 

If an invalid section shape code is specified it is not possible for the 
program to properly compute section end geometries; as a result an error 
occurs: 

**ERROR -- INVALID SECTION SHAPE CODE WAS SPECIFIED AS (input code) 
(SHOULD BE 1, 2 OR 3). THIS SECTION WILL BE SKIPPED. 

An error which arises during computation and causes a non-positive total 
pressure loss for a given section prevents completion of the case analysis and 
gives rise to an error message: 

**ERR0R “ SOME INCORRECT COMBINATION OF INPUTS OR UNANTICIPATED 
SITUATION HAS CAUSED AN INVALID (NON-POSITIVE) TOTAL LOSS LEVEL. 

RECHECK SECTION (section number) INPUT DATA. 

If the maximum allowable number of circuit sections written into the program 
is exceeded by placing too many section cards together in one case, or without 
termination, title, and master cards between cases, this diagnostic will 
appear: 

MAXIMUM LIMIT ON THE NUMBER OF SECTIONS ( . . . (maximum allowable number 
of section)...) HAS BEEN REACHED. EITHER A CASE TERMINATION CARD HAS 
BEEN OMITTED (ALONG WITH TITLE AND MASTER CARDS TO BEGIN A NEW CASE) 

OR THIS CASE IS TOO LONG FOR THE PROGRAMMED ALLOWABLE NUMBER OF 
SECTIONS. THE CASE HAS BEEN TERMINATED AT THIS POINT. 

In this Instance, the inputs from the group of sections for which the limit 
was exceeded will be tabulated and the remaining section inputs will be 
evaluated and tabulated. If the user fails to cause the test section blockage 
amounts specified on the master control card to coincide with that of the test 
section card, erroneous analysis may result since inconsistent flow areas 
would be calculated. The section card value will be used (since the discrep- 
ancy may be desired) and this notice is given: 

**NOTE — TEST SECTION BLOCKAGE FROM SECTION CARD INPUT (...(section 
input value) . . .PERCENT) DOES NOT EQUAL THAT OF THE MASTER CARD INPUT 
(...(master input value) .. .PERCENT) . CHECK DATA DECK. SECTION CARD 
VALUE WILL BE ASSUMED AS CORRECT AND EXECUTION WILL CONTINUE. 

An invalid section type code will cause a section to be skipped and a message 
to be printed: 

*ERR0R — INPUT SECTION TYPE CODE (CARD COLUMNS 3 AND 4) CALLS 

INVALID SECTION TYPE. DATA CARD IGNORED.** 
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Any input errors were deemed justifiable cause for judgment as an incomplete 
case. As a result, reliable overall and summary information cannot be calcu- 
lated. To assist the user in locating the error(s), the input values will be 
forced to be tabulated and the following explanation appears: 

***DUE TO ERROR(S) IN INPUT CARD(S), VALID SUMMARY INFORMATION IS NOT 

AVAILABLE. REFER TO THE TABULATION OF INPUT DATA ON THE FOLLOWING 

PAGES. CORRECT THE ERROR(S) AND RESUBMIT THIS CASE. SUBSEQUENT CASES 

WILL NOT BE AFFECTED. 

4. Possible errors lacking diagnostics: Certain potential problem areas 

remain unprotected by diagnostic and error-recovery systems. 

No special provision was made for two test sections in the same circuit 
case. As long as the blockage values for both test sections match the one 
from the master card, no message will be printed. In any event, the execution 
will not be terminated. The test section shapes and dimensions from the 
master card are not checked against those of the test section card. Although 
a mismatch of these values could cause a mass-flow error, including and 
enforcing such a check could inhibit any meaningful tandem-test section cases. 
These problems could be avoided, however, by naming only one working section 
as a test section and referring to the other by general type. 

Also, there was no provision for checking the specified tunnel type 
against the types of sections actually used (e.g., checking a non-return, or 
open-test-section tunnel for exit or open-throat test section input cards) . 
This check is not critical and was left to the user. 

One error-check was not included due to the program complications which 
would have resulted* If a case termination card is omitted at the end of a 
case and a computer-system control card or a title card is encountered, the 
error will be disastrous due to mismatched format types. Execution and 
calculations will be immediately aborted by the computer. 


Test Case 

The NASA-Ames Research Center 40- by 80-Foot Wind Tunnel was used as an 
example of a typical wind tunnel. This tunnel, illustrated in figure 10, is 
of the single-return, closed-test-section, continuous- running type. It has a 
flat-oval test section 12.2 m (40 ft) high by 24.4 m (80 ft) wide and is 
powered by six 12.2-m (40 ft) diameter, six-bladed fans. It has an eight-to- 
one overall contraction ratio and uses multiple-circular-arc type turning 
vanes in each of the four 90® corners. 

A complete list of the test case inputs and computed information outputs 
are presented in figure 11. The machine computing time for this test case 
(without plots) was about 7 sec on an IBM 360/67. 

Although this test case was not an exhaustive exercise of all possible 
tunnel components, it does include most of the basic section types: diffusing 

test section, single-duct contraction and diffuser, constant-area single duct. 
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constant-area corner with turning vanes, and multiple-duct fan sections 
(contraction, constant-area annulus, and diffuser). Examples of other types 
of components are shown in the sample cases which follow. 


DISCUSSION AND APPLICATIONS 


Wind tunnel energy ratio, required power, and operating velocity are 
interdependent. The energy ratio is affected by velocity through the effect 
of velocity on the Reynolds number. The required drive power, influenced 
directly by operating velocity and inversely by energy ratio, is also con- 
trolled by the fan system efficiency which is often only an estimated quantity. 
Any estimate of operating velocity for a given power level is, then, dependent 
on the basic efficiency of the circuit (energy ratio) and drive system effi- 
ciency, assuming the best power estimate available to be that delivered to the 
fans. This interdependency means that an error in the prediction of energy 
ratio (and/or in the estimation of fan efficiency) will cause corresponding 
errors in power and velocity estimates. 

These errors resulting from an erroneous prediction of the circuit energy 
ratio can be found from the relationship governing required power, test sec- 
tion velocity, and energy ratio, assuming given motor electrical and fan 
efficiencies. For a fixed test section velocity, 

^^REQUIRED ^ 1 _ ^ 

^REQUIRED 1 _ 

ER 

and, for constant power, an error in energy ratio yields the performance 
penalty 



The expected true power and velocity levels can thus be obtained from the 
performance estimate: 


^REQUIRED “ 



^REQUIRE 


for a given set of test section conditions, and from 

1/3 

'^°Estimate 

for a given level of required power. 
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This adjustment of the estimated performance values is pointless for a 
known, existing wind tunnel, but necessary for new, or proposed facilities. 
Before the adjustment can be made the probable error in the energy ratio esti- 
mate must be determined. It is desirable, therefore, to consider several 
existing facilities of different circuit types in order to gain a degree of 
confidence in the performance estimation routine. 


Results 

The input parameters and output performance values for the several sample 
cases, other than the test case shown in figure 11, are compiled in appendix D. 
The estimated energy ratios for the seven sample wind tunnels are presented in 
table 6. The corresponding sketches for all these sample tunnel circuits are 
shown in figures 10 and 12. 

The actual energy ratios for the first three wind tunnels presented in 
table 6 were estimated from the best available information on fan and electri- 
cal efficiencies from known input power levels. The actual performance of the 
other four facilities was taken from measured data. 

The test case and first sample case was the circuit of the NASA-Ames 
Research Center 40- by 80-Foot Wind Tunnel as described previously in the test 
case discussion. The predicted energy ratio for this rather conventional 
tunnel was only 1 percent higher than the actual value when new. 

The performance of the NASA-Ames 7- by 10-Foot Tunnel was predicted at a 
slightly optimistic level. However, this tunnel is one with several known 
problems which complicate the prediction process. With the air exchanger oper- 
ating, the primary diffuser is known to have some local flow separation, 
having been designed at a 6® equivalent cone angle, an angle too great for its 
cross-sectional shape and length (see fig. 1). Also, the drive fan is stalled 
from the centerbody out to about 45 percent of the fan radius, causing some 
back-flow along the nacelle centerbody. (The impact of the stall on the fan 
efficiency has only been estimated; it was assumed that the fan efficiency 
would suffer by about an additional 10 percent.) In spite of these things, 
the predicted energy ratio was only about 3 percent too high relative to the 
original value of approximately 7.85, both values taken in the zero-air- 
exchange configuration. This agreement may indicate that much of the above- 
mentioned off-design performance is triggered by the air exchanger operation 
and is not as significant with the air exchanger closed- Although insuffici- 
ent data are available to resolve this question, the fact remains that the 
prediction accuracy, for the stated conditions, was good. 

The Lockheed-Georgia Low-Speed Wind Tunnel employs a tandem test section 
design. For this analysis, the larger, V/STOL test section was used as the 
only reference station. Because of the two area restrictions to cross sec- 
tions smaller than that of the reference area (those of the smaller test sec- 
tion and of the fan), the tunnel efficiency would be expected to be low. 

(This in no way reflects on the tunnel’s usefulness as a research tool or on 
its design or capabilities. The ”low efficiency*’ value results only from the 
point of reference used in the calculations.) In other than these features 
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the facility Is basically of conventional design. The computerized 
performance prediction was in error by less than 2 percent from the true value 
of about 1.10. 

The Indian Institute of Science 14- by 9-Foot Tunnel at Bangalore stands 
out among non-return wind tunnels as a facility with an unusually high energy 
ratio. Although the determination of circuit dimensions for the program input 
was somewhat hampered by the limitations of small drawings, the estimated 
energy ratio was within 1 percent of the true facility value of 6.85. It is 
interesting to note that the fan performance data of reference 23 would indi- 
cate a fan design efficiency of about 69 percent. The power requirement calcu- 
lations based on energy ratio and test section maximum velocity, however, show 
that the fan efficiency must be higher than was expected; in fact, greater 
than 90 percent. 

The Hawker Siddeley 15-Foot V/STOL Tunnel at Hatfield, England was con- 
structed under economy constraints and is a compact, cost-effective facility. 
The estimated performance was about 1.6 energy ratios higher (i.e., more 
optimistic) than the actual value of 2.38. This is an error of about 67 per- 
cent. The primary performance difference was probably caused by the fan sys- 
tem. The losses of the ducting in this area are difficult to predict because 
the area changes are not gradual and are even difficult to define. 

The University of Washington 8- by 12-Foot Double-Return Tunnel has a 
surprisingly high measured energy ratio of 8.3. This would indicate a very 
carefully designed circuit powered by carefully designed fans. The perform- 
ance estimate produced by the computer program is lower than the actual energy 
ratio by about 13 percent, showing that the achieved performance level is 
higher than would normally be expected. 

The NASA-Langley Research Center 30- by 60-Foot Open-Throat Tunnel is 
unusual in configuration, having a double-return system with the twin fans 
located less than two fan-diameters downstream of the test section. The loca- 
tion of the data point for this tunnel on the diffuser design curves of fig- 
ure 1 would not indicate that any diffuser-related problems should be expected 
forward of the fans. The diffuser between the fans and the first corner, how- 
ever, does have a rather large equivalent cone angle (more than 8®). If the 
fans cause or contribute to diffuser flow problems (see the Cautionary Design 
Guidelines for Diffusers) and if those problems lead to corner flow ineffici- 
ency in a region critical to overall performance, then the circuit energy 
ratio may be well below the normal estimated level. Although it is not clear 
whether this is the case in the NASA-Langley 30- by 60-Foot Tunnel, the perform- 
ance estimate was about 27 percent higher than the actual value of about 3.71. 


Evaluation 

To summarize what may be learned from the sample cases: 

1. The Ames 40- by 80-Foot and 7- by 10-Foot Tunnels and the Lockheed- 
Georgia Low-Speed Tunnel, although at opposite ends of the energy ratio 
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spectrum, are all basically standard, single return, closed-test-section 
facilities; the computer program estimates of actual performance were good. 

2. The Indian Institute of Science Bangalore tunnel, being of the non- 
return variety, is a different and less common type of facility; the computer 
program closely estimated its actual performance. 

3. The University of Washington double-return tunnel is a third major 
circuit type; the program produced a reasonably accurate prediction of its 
performance. 

4. The Hawker Slddeley V/STOL and Langley 30- by 60-Foot Tunnels are 
examples of facilities which may have flow problems due to too-rapid area 
changes and, as a result, lower than optimum performance levels for their 
respective circuit types. For these tunnels, because of their flow quality 
and not because of their circuit types, the program provided a poor estimate 
of actual performance. 

Based on these results one thing is immediately clear: the performance 

of a wind tunnel of conventional, conservative design can be evaluated accu- 
rately. On the other hand, the performance of a tunnel whose design generates 
or contributes to flow problems (separation or grossly non-uniform) will be 
overestimated by the loss equations and computer program. 

Flow peculiarities and off-optimum designs, even though seemingly only 
slight, can cause operational performance to fall significantly below the pre- 
dicted levels. Such problems can be expensive whether considered in terms of 
modifying the facility or in such terms as reduced testing capability and 
increased power costs. Judicious, iterative use of the estimation techniques 
presented in this report, simplified by computerized automation, can lead to 
the improvement of an existing facility through guidance of design changes or 
to the optimization of a proposed new wind tunnel design. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, January 8, 1976 
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APPENDIX A 


NON-STANDARD FUNCTIONAL FORMS 


Due to the nature of this analysis, certain of the local flow-state, sec- 
tion loss, and summary parameter formulas were used in a form more convenient 
than that usually found in the literature. The relationships which were 
altered or derived are outlined on the following pages. 

Local Flow-State Parameters 

The calculation of several local parameters was based on the local Mach 
number, determined from the relationship between the local area and the area 
for choked flow: 


A 

_a 

A 


Y+1 

^ + 1 ^2 (y-1) 




Y+1 

2(y-1) 


(ref. 18, p. 6) 


Solving for the area for choked flow, knowing the test section area and Mach 


number, 



Mach number- Another form of the same area relationship, 



2 

Y+1 



Y+1 



(ref. 17, p. 126) 


can be rewritten to produce a polynomial equation in Mach number which may be 
solved by Newton’s method if the areas are known: 





2 

Y + 1 




Y + 1 \Y + 1 



m2_ 


Y + 1 

Y - 1 



= 0 


26 



Reynolds number- The local Reynolds number was calculated based on other, 
known, local conditions and from basic principles: 


RN = 


y 


PVA = PoAoVo 

(conservation of mass) 

/ T \0 • 76 



(ref. 18, p. 19) 


(ref. 18, p. 4) 


76 


Friction coefficient- The Reynolds number-friction coefficient function 
used was 

= 2 log^Q(KN/x) - 0.8 (ref. 6 , p. 70) 

A Newton's method solution was performed on a rewritten form of the equation: 

[log^^CXRN^) - 0.8]"2 - a = 0 


Section Pressure Losses 

The losses for some types of sections were derived in forms not found in 
the literature. For others, a curve-fit of data points or a simplification of 
analysis was performed. 

Comers (constant area)- The frictional and rotational losses through 
turning vanes are additive : K — ^frxCTION ^ROTATION * Also , 

1 2 

Friction = ■3 ^tv Dotation '3 %v p* 


Assuming that the frictional loss value has a form similar to that for a flat 
plate, then at 90°; 


^FRICTION - 3 %Vgo 


0.455B 

(log^Q 


(ref. 6 , p. 527) 
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Thus, the constant B is dependent on the turning vane loss constant and the 
reference Reynolds number at which that constant was determined: 


B = 


1 / \ 2 . 5 8 

0.455 


Therefore, 


*^FRICTI0N 



\ / 


Since the rotational term is assumed independent of Reynolds number, 

^ROTATION ” additional complication of loss parameter varia- 

tion with turning angle is presented in figure 3 for a loss parameter at 90® 
equal to 0.15. It was assumed that the relationship between the actual and 
reference loss constants is linear: 


Ktv ~ ^TVso 



where f((f)) is the functional relationship plotted in figure 3. 
turning vane loss function then becomes 


K = K-pv 



/loglO RNrefY ‘^^11 

V J) 


The complete 


Comers (diffusing)- Diffusing corners were treated as vaned diffusers 
with the addition of rotational losses dependent on the turning angle. The 
expansion and frictional losses used were those for a vaned diffuser: 


%ANED DIFFUSER “ "*■ [0.006(20 - 21.5®)u(20 - 21.5”)][ ^^^ 


The rotational loss is as for a constant-area corner where 


^ROTATION ■ 3 %V ~ "3 ^Vgo 



The diffusing corner loss function is then 

2 

K * {0.3+ [0.006(20 - 21.5°)u(20 - 21.5°)][ f K.pv 


where u(20 - 21.5°) is the unit step function. 
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Dtffusers^ The diffuser losses are due to Doth friction and expansion. 
The friction term may be derived theoretically from 


^FRICTION 


' FRICTION 


Xpy2 

K PlVi'Dh 


Making the simplifying assumptions that the density and the friction coeffi 
cient are approximately constant and applying conservation of mass, 

^S2 

^FRICTION “ J TTa 
which, for a conical diffuser, becomes 


_ IbAi^X 
^FRICTION “ ^ J ^ 


and transforming variables from surface to centerline distances, 

„ 16 Ai^X dx 

(D, + 2x tan 9)5 

Completing this integration the friction loss becomes 


^FRICTION - 8 sin 6 




The influence of the expansion term is given by 

^ “ Expansion ^friction 


Thus, it may be rewritten; 


Expansion Friction 

rrri? 


(‘ - 


Ar2 - 1 
AR2 






EXP L8 sin 0 




^EXP 


Expansion 

/ AR - 
V AR / 


^ ■ ^FRICTION 
EXP ■ /AR - I f 
V AR / 


The expansion loss parameter curves shown in figure 4 were determined using 
the approximation 


Xrxp 


X - KfRICTION (CONICAL) 
MR - 
V AR / 


XeXP 


^■[s sin e AR^)] 


/ AR - 1 
V AR 




and figure 5(a) of reference 9, which shows complete diffuser losses plotted 
as functions of equivalent cone angle and independent of area ratio for circu- 
lar, square and rectangular, and two-dimensional diffusers. (This Implies an 
assumption that the expansion part of the losses is dependent only on cross- 
sectional shape, the extent to which the diffusion takes place in only one 
direction, and the equivalent cone angle.) Thus, the complete loss for 
diffusers is given as 


using K£jg> from figure 4. 

Exit- The kinetic energy loss at an exit of a non-return wind tunnel was 
derived from basic compressibility relationships and with the assumptions that 
the exit flow static pressure is equal to the atmospheric pressure and that 
the exit velocity is uniform. 

^ . [l- M^)]- 

Rewriting, the local total pressure is 


(rewritten 

from ref. 17, p. 53) 


Pt ■ «')]''■* 
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Also, since Ap^ * ® P.p “ P» the total loss parameter is 


^ I YPM2 

since 

q = pV^ » ^ YPM^ (ref. 17, p. 55) 

Simplifying, the exit loss becomes 

Flow straightenera : airfoil memhera (thick)- Thick flow straightener 

losses were assumed to be made up of two parts: contraction and subsequent 

diffusion: 


^ “ ^CONTRACTION DIFFUSION 


The contraction was estimated as being about 30 percent of the length of the 
straighteners: 


0.32X(0.30L) 
^CONTRACTION = ?T 


0.096XL 

•^CONTRACTION “ Dij 

The diffusion portion was based on the aft 70 percent of the length and on the 
exit and minimum areas for the computation of the area ratio and equivalent 
cone angle. As for a vaned diffuser, 


•^DIFFUSION - 


{0.3+ [0.006(20 - 21.5°)u(20 - 21.5°) -) 


Hence the loss for thick flow straighteners becomes 


2 

K = 0.096 ^+ {0.3 + [0.006(20 - 21.5°)u(20 - 21.5°)]{( ^^ 
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Internal flow obstruction: drag item- The loss due to internal structure 

may be derived from the relationships governing power losses; 

PinpuTdeag 2 ;;^ 

and Pdrag “ ^Ve = (l/2)pV^SCjje, where e is the factor accounting for addi- 
tional effects on downstream sections. Since ^INPUTj)g;^Q ” ^DRAG» loss 
becomes 


and therefore 


K . A A c c ilL 
° % Ao ^ Po^o 


T P - 

A P P o^o^o 


K = 




S 

A 


e 


PF 

P 


Since in general the flow density at the fans is unknown at the time a given 
section loss is calculated, and since for incompressible flow the density 
ratio is unity, the ratio of the densities at the fan and the local station 
was assumed as unity for the analysis. (If the user prefers not to make such 
an assumption, an approximation of the ratio may be made by way of a change in 
the downstream influence factor e.) The loss due to a flow obstruction is 


K = 



e 


Varied diffusers- The expansion and friction losses for vaned diffusers 
were combined into one parameter which is reasonably independent of area ratio 
and is presented in figure 6. The loss curves shown were approximated by a 
two-segment , straight-line curve fit so that, for vaned diffusers 


K = 


V AR / 


2 


and 

K = -10.3+ [0.006(20 - 21.5°)u(20 - 21.5°)]}( ^^~ 

where u(29 - 21.5°) is the unit step function. 

Loss value transferred to reference location- The change of reference for 
loss values is defined as 
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Using the law of conservation of mass, this may be rewritten In terms of areas 
and Mach numbers: 


pVA - 



MdlTy pressure differential- The pressure across a section of wall was 
determined from the exterior atmospheric pressure, internal static pressure, 
and cumulative pressure losses through the circuit. Since the wall pressure 
differential for a given section is Apy^ = “ Pi 

p 

[i+(^^Mi2j|]Fr 

and, using the test section as the reference location, 

i 

PTi = PTsc ■ 2 Koj 


The wall pressure differential may be written as 


Input power required- The power input to the flow required for operation 
of a wind tunnel circuit having predetermined losses was calculated from the 
pressure rise required at the fans, with the simplifying assumption that the 
static and total pressure rises across the fan are equal. 

P INPUT “ 

^ PpPoVo 

PiNPUT = APTpApVF ppp^AoVo 


Considering conservation of mass, 


Also, 


Thus, 


PiNPUT “ 


^Plp^o^o 



N 

Apt = Ap = E 
^ i=l 
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APPENDIX B 


NUMERICAL FUNCTION-APPROXIMATIONS 

The formulas that follow resulted from curve-fitting and/or interpolation 
techniques applied to certain functions arising from the loss analysis. 

Corners 

The corner loss parameters for corners with and without turning vanes are 
shown in figure 3. For a corner with vanes, using least-squares polynomial 
curve-fitting techniques, the turning vane loss function of figure 3(a) 
becomes, for 0® ^ 0 :< 30®: 

Ktv = 1.395066x10"^ + 5.672649xl0"‘' 0 

+7.081591x 10"5 0^ + 1.394685x10"® 0® 

-4.885101x10"® 0** (Bl) 

and for 30° < 0 90°: 

K.j.y * -1. 605670x10" 1 + 1.446753x10"^ 0 

-2.570748x10"'* 0^ + 2.066207x10"® 0® 

-6.335764x10"® 0‘* (B2) 

For a corner without turning vanes the local loss function of figure 3(b) was 
found using a least-squares polynomial technique and is given by 

K = 4.313761x10"® - 6.021515x10"'* 0 

+1.693778x10"'* 0^ - 2.755078x10"® 0® 

+2.323170x10"^ 0'* - 3.775568x10"® 0® 

+1.796817xl0"ll 0® (B3) 

For all the above equations, 0 is the flow turning angle in degrees. 

Diffusers 

The determination of the diffuser loss parameter is a complex operation. 
It depends on the cross-sectional shape and equivalent cone angle of the sec- 
tion. For a conical diffuser the expansion functions are, for 3® ^ 20 ^ 10®: 
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Kfxp^. , = 1.70925x10-1 - 5.84932x10-2 (29) 

^■'^Circular 

+8.14936x10-3 (26)2 + i. 34777 x 10 -'* (26)3 


-5.67258x10-3 (20)'* - 4.15879x10-2 (20)3 
+2.10219x10-2 (26)3 


(B4) 


for 0“ < 20 < 3°: 


Kt 7 vp = 1.033395x10-1 - 1.19465x10 2 ( 26 ) 

Circular 

and for 26 > 10° ; 

Kpyp = -9.66135x10"2 + 2.336135xl0"2 (26) 

i^'*^Circular 

For a square cross-section diffuser the expressions are, for 3 i 29 ^ 10 


(B5) 


(B6) 


Kexp 


Square 


1.22156x10-1 - 

2.29480x10-2 (26) 



+5.50704x10-3 

(26)2 _ 4.08644x10"'* 

(20)3 


-3.84056x10-3 

(20)'* + 8.74969x10“3 

(26)3 


-3. 65217x10- 

(20)3 


(B7) 


for 0° < 29 < 3°: 
KeXP 

and for 26 > 10°: 
l^EXP 


Square 


= 9.62274x10-2 - 2.07582x10-3 (26) 


(B8) 


= -1.321685x10-1 + 2.93315x10 2 (20) (B9) 

Square 

For a two-dimensional diffuser with a square upstream-end cross section the 
expansion loss functions are, for 3® ^ 26 _< 9®: 


Kpyp = 3.23334x10-1 - 5.82939x10-2 (20) 
2HRectangular 

-4.97151x10-2 (26)2 + 1.99093x10-2 (26)3 
-1.98630x10-3 (29)** + 2.06857x10-^ (20)3 
+3.81387x10-3 (29)3 


(BIO) 


for 9° < 26 5 10°: 

Kexp 


= 5.72853 - 1.21832 (29) 

^Rectangular 


+7.08483x10-2 (26)2 


(BID 
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for 0° < 20 < 3°: 


KeXP^d 


Rectangular 


and for 20 > 10®: 
^EXP 


^^Rectangular 


1.0x10"^ - 5.333333 x 10'3 (20) 
-1.36146 + 1.986460x10-^ (20) 


(B12) 

(B13) 


Since the expansion function for a two-dimensional diffuser with circular 
sides was not given (and is not defined), it was assumed for computational 
purposes that this value would be the same fraction of that for a two- 
dimensional recta^igular diffuser as the loss of a conical is of that for a 
three-dimensional square diffuser: 


/ ^EXPcircular\ 

^^2DRg(,tangular \ ^EXPgqygye j 


For cross-section shapes somewhere between rectangular and circular, such as 
flat oval (flat ceiling and floor with semi-circular sidewalls), or for dif- 
fusers which have one end rectangular and one end either circular or flat oval, 
a loss value between that for circular and rectangular may be more appropriate; 
thus, 

•^^Rectangular “^^Circular 

2 


and 


,, _ ^^^Square ^^^Circular 

^EXP.r, 2 

^^Average 


The extent to which a diffuser is planar in nature was computed from the 
ratio of the changes in size of the two characteristic dimensions from end to 
end: 


6g = smaller of 


^2 ~ hi 
W2 - 


or 


W£ “ Wj 

h2 - hi 


or if the ratio is negative. 


Then, based on the geometries of each end, the basic loss constant, 
may be selected from ^EXPcircular* ^^EXPsquare addi- 
tional loss fact. ^EXP^jjjjitional’ selected from the corresponding 
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. Finally, the applicable 


Kexp.t) , . ^EXPod, ^EXPod„ 

■^^Circular ^^Average ‘^‘^Rectangular 

diffuser expansion loss coefficient is given by 

KeXP = ^EXPB^ 3 ic + (1 - ^s)(KEXPAdditional 


~ Eexp 


Basic 


.c) 


(B14) 


Mesh Screen 

The mesh screen Reynolds number sensitivity factor plotted in figure 5 
can be expressed in functional form as, for 0 i RN < 400: 


Ern ■ 




100 


+ 1.01 


(B15) 


and for RN >. 400: 


Krn = 


Vaned Diffusers 

The vaned diffuser loss coefficient functions plotted in figure 6 were 
approximated by two line segments; for 26 < 21.5°: 

Kv = 0.3 


and for 21.5° ^ 29 < 90°: 

= 0.3+ [0.006(26 - 21.5°)] 

Thus, over the entire range of equivalent cone angles of interest, 

Kv = 0.3+ [0.006(26 - 21.5°)u(29 - 21.5°)] (B16) 

where u(20 - 21.5°) is the unit step function. 


38 



APPENDIX C 


COMPUTER PROGRAM F0RTRAN CODES 

The following pages contain the F0RTRAN codes developed to implement the 
wind tunnel performance analysis techniques presented in this report. 

The Notation section explains the variable names used in the program. 
(Note that in the notation sections, as throughout this report, all letter O's 
occurring in F0RTRAN names are shown with slashes, as 0; all number zeros are 
shown unslashed.) This notation section is similar to that for engineering 
symbols presented in the main body of the report, but this section was changed 
in two respects. First, it was rearranged alphabetically by F0RTRAN variable 
name. Second, it was expanded to include many variable names which were not 
used elsewhere and which have significance only in the context of the computer 
program. The ’’titles” shown in parentheses in the first column of this nota-- 
tion section are column heading titles which appear on the program output 
pages. 

Immediately following the Notation are the listings of the six actual 
F0RTRAN program codes: the main program (PERF0RM) and the five subroutines 

(DATACK, SPEED, FRICTN, 0UTPUT, AND PL0TIT) . Each program routine page is 
titled and numbered for clarity. The last seven columns of each line on each 
page contain a two-letter program routine name abbreviation and a line 
sequence number (in ten-count increments). Thus, the user can know at a 
glance to which routine (and where within that routine) any given line or 
instruction belongs. Each instruction line in the program is uniquely 
identified. 
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NOTATION (F0RTRAN) 


F0RTRAN name Engineering 

and/or (title) symbol Description 


A 


A 


cross-sectional area of local section, 
m2 (ft2) 


All 


cross-sectional area of individual duct 
at upstream end, m^ (ft^) 


AI2 

AL ApLow 


cross-sectional area of individual duct 
at downstream end, m^ (ft^) 

cross-sectional area of local flow, m^ 
(ft2) 


AMACH M 


local Mach number 


AMACHl Mach number at section upstream end 

(MACHl) 

AMACH2 Mach number at section downstream end 

(MACH2) 


AR AR 

(AR,CR) 


ratio of cross-sectional areas at 

upstream and downstream ends of section 


ASL 


speed of sound in moving flow at local 
section, m/sec (ft/sec) 


ASTAR A^ 


cross-sectional area for sonic flow at 
specified flow conditions, m^ (ft^) 


ASO 


speed of sound in moving flow at upstream 
end of test section, m/sec (ft/sec) 


AT arp 

AVGPWR 

AO Aq 

A1 Ai 

(AREAl) 


speed of sound in still gas, computed at 
total (stagnation) conditions, m/sec 
(ft/sec) 

average power consumed by each drive fan 
at specified conditions: PWR0P/ENFAN, 

W (hp) 

cross-sectional flow area of test section 
at upstream end, (ft^) 

cross-sectional flow area of section at 
upstream end, m^ (ft^) 
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F0RTRAN name 
and/or (title) 


Engineering 

symbol 


Description 


A10AO 

(Al/AO) 

A2 

(AREA2) 

A20AO 

(A2/A0) 

BLKAGE 

(BLKGE) 

CD 

CH0RD 

D 

DATA 
DELP 
(D EPS) 

DFAN 

DH 

DHL 


A2 


C 


D 


c 


V 


D 


^PWi 


Ae 


ratio of local section upstream area to 
test section area, m^ (ft^) 

cross-sectional flow area of section at 
downstream end, m^ (ft^) 

ratio of local section downstream area to 
test section area, m^ (ft^) 

blockage to flow in local section (at 
upstream end for all applicable sec- 
tions except fan contraction, for which 
it is at downstream end) , fraction of 
local area 

blockage to flow in local section (at 
upstream end for all applicable sec- 
tions except fan contraction, for which 
it is at downstream end) , percent of 
local area 

drag coefficient of flow obstruction, 
drag 
qS 

chord of turning vanes, m (ft) 

diameter of circular duct, m (ft) 

data array of master, section, and 

termination card floating-point inputs 

local pressure difference across wind 
tunnel wall, N/m^ (Ib/ft^) 

increment of flow-obstruction downstream 
influence factor greater than unity: 

£ - 1, (greater than or equal to zero) 

drive fan diameter, m (ft) 


D 


h 


hydraulic diameter: 

4 X (cross-sectional area) 
perimeter 


m (ft) 


hydraulic diameter of single cell in flow 
straightener , m (ft) 
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FORTRAN name 
and/or (title) 

Engineering 

symbol 

Description 

DHUB 


diameter of drive fan hub and/or spinner, 
m (ft) 

DHO 


hydraulic diameter of test section, m 
(ft) 

DHl 


hydraulic diameter of upstream end of 
local section, m (ft) 

DH2 


hydraulic diameter of downstream end of 
local section, m (ft) 

DMESH 


diameter of mesh element in woven-mesh 
screen, m (ft) 

D1 


diameter of upstream end of circular 
section, m (ft) 

D2 


diameter of downstream end of circular 
section, m (ft) 

EK 

(DP/QL) 

K 

local total pressure loss of section: 
APt 
q 

EKADD 

^EXP^dd i t iona 1 

additional diffuser expansion loss factor 
due to more diffusion in one plane than 
in another (i.e., partially two- 
dimensional diffusion) 

EKBASE 

^EXPEasic 

basic diffuser expansion loss factor for 
purely three-dimensional diffusion 

EKC 

^EXPcircular 

expansion loss value for conical 
diffusers 

EKCNTR 

^CONTRACTION 

local total pressure loss from contract- 
ing portion of thick-airfoil flow 
straighteners 

EKCSAV 


estimated expansion loss coefficient for 
three-dimensional , combination circular 
and square cross-section diffuser 

EKD 

^DIFFUSION 

local total pressure loss from diffusing 
portion of multi-loss-type sections 

EKEXP 

(KEXP) 


net expansion loss coefficient for 
diffusers 
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F0RTRAN name 
and/or (title) 

Engineering 

symbols 

Description 

EKMESH 

(KMESH) 

Kmesh 

mesh screen-type loss constant 

EKS 

^EXPgqygJ-g 

expansion loss value for three- 

dimensional expansion in square cross- 
section diffusers 

EKSTRT 


local total pressure loss coefficient due 
to strut drag in fan section 

EKTE 


local total pressure loss parameter for 
corners without turning vanes 

EKTE90 
(KT 90) 


vaneless-corner loss parameter for given 
corner at a 90® turn 

EKTV 


turning vane loss coefficient 

EKTV90 
(KT 90) 


turning vane loss parameter for given 
vanes at a 90® turn 

EKV 

Kv 

local total pressure loss coefficient for 
vaned diffusers 

EKO 

(DP/QO) 

Ko 

section total pressure loss referred to 

Api 

test section conditions: 

‘lo 

EKl 


local total pressure loss coefficient due 
to diffusion and vanes in a diffusing 
corner 

EK2 


local total pressure loss coefficient due 
to rotational flow in a diffusiing 
corner 

EK2DC 

*^EXP2Dcircular 

estimated expansion loss coefficient for 
hypothetical , two-dimensional diffusion 
with circular sides: 



^^Rectangular \ ^^^Square / 

EK2DCS 

E^2DAverage 

estimated expansion loss coefficient for 
two-dimensional diffuser with cross- 


section shape of some square/circular 
hybrid 
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F0RTRAN name 
and/or (title) 

Engineering 

symbol 

Description 

EK2DR 

EXP2DRectangular 

expansion loss coefficient for two- 
dimensional rectangular cross-section 
diffusers 

EL 

(L) 

L 

centerline length of section, m (ft) 

ELC 


length of contracting portion of thick- 
airfoil flow straighteners, m (ft) 

ELD 


length of diffusing portion of thick- 
airfoil flow straighteners, m (ft) 

EL0DH 

(L/DH) 


length-to-hydraulic-diameter ratio of 
flow stralghtener cell 

EMDATA 


data array containing master-card 
floating-point inputs 

EMF 


Mach number at the fan section 

EMU 

y 

flow viscosity, N sec/m^ (lb sec/ft^) 

EMUSTD 

*^8td 

standard-day value of viscosity, 
N sec/m^ (lb sec/ft^) 

EMUT 

Vt 

reference viscosity at a known tempera- 
ture, computed for a still gas (stagna- 
tion conditions), N sec/m^ (lb sec/ft^) 

EMWRIT 


master card output array containing data 
and/or annotatlon(s) 

EMO 

Mo 

Mach number at upstream end of test 
section 

ENDATA 


data array containing section-card 
floating-point input 

ENDUCT 


number of ducts in multiple-duct sections 

ENFAN 


number of fans in fan drive section 

ENITEM 


number of drag or blockage items in each 
local duct 

ENU 

V 

kinematic viscosity of gas, m^/sec 
(ft^/sec) 
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F0RTRAN name 
and/or (title) 

ENWRIT 


Engineering 

system 


Description 


EPS 

ER 

ETAFAN 

(ETA) 

ETWRIT 

FAC 

FACT 

FAF0 

FAR 

FDHC 

FDHF0 

FDHR 

FEKC 

FEKCH 


section-card output array containing data 
and/or annotation (s) 

e flow-obstruction downstream Influence 

factor (greater than or eqvial to unity) 

ER energy ratio: ratio of energy of flow at 

test section to the output energy of 
the fans 

rip fan aerodynamic efficiency, percent 


case termination-card output array con- 
taining termination request de-codings 

function defining the area of sections 
with circular cross sections 

scaling factor for plot size 

function defining the area of sections 
with flat-oval cross sections (flat 
floor and celling, semi-circular walls) 

function defining the area of sections 
with rectangular cross sections 

function defining the hydraulic diameter 
of sections with circular cross 
sections 

function defining the hydraulic diameter 
of sections with flat-oval cross 
sections 

function defining the hydraulic diameter 
of sections with rectangular cross 
sections 

function defining the diffuser expansion 
loss for three-dimensional, circular 
cross-section diffusers 

function defining the diffuser expansion 
loss for three-dimensional, circular 
cross-section diffusers at high diffu- 
sion angles (TH2 > 10®) 


A-5944 
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F0RTRAN name 
and/or (title) 


Engineering 

system 


FEKCS 


FEKS 


FEKSH 


FEKSS 


FEKO 


FEK2DL 


FEK2DU 


FKTE 


FKTVl 


FKTV2 


FTH 


A6 






Description 

function defining the diffuser expansion 
loss for three-dimensional, circular 
cross-section diffusers at small diffu- 
sion angles (TH2 < 3°) 

function defining the diffuser expansion 
loss for three-dimensional, square 
cross-section diffusers 

function defining the diffuser expansion 
loss for three-dimensional, cross- 
section diffusers at high diffusion 
angles (TH2 > 10°) 

function defining the diffuser expansion 
loss for three-dimensional, square 
cross-section diffusers at small diffu- 
sion angles (TH2 < 3°) 

function defining the change-of-reference 
station for total pressure losses from 
local section to test section 

function defining "two-dimensional" 
(rectangular) diffuser expansion loss 
for low diffuser angle range (TH2 < 9°) 

function defining "two-dimensional" 
(rectangular) diffuser expansion loss 
for high diffuser angle range (TH2^9°) 

function defining corner turning loss 
parameter EKTE for corners without 
turning vanes (based on a value of 
EKTE = 1.80 at PHI = 90°) 

function defining turning vane loss param- 
eter EKTV (based on a value of 
EKTV = 0.15 at PHI = 90°) for lower 
turning angle range (PHI i 30°) 

function defining turning vane loss 
parameter EKTV (based on a value of 
EKTV = 0.15 at PHI = 90°) for upper 
turning angle range (30° < PHI i 90°) 

function converting diffuser equivalent 
cone angle, TH2, in degrees to half- 
angle, TH, in radians 


F0RTRAN 

and/or 

name 

(title) 

Engineering 

symbol 

Description 

FTH2 



function defining diffuser equivalent 
cone angle, TH2 

G 


Y 

specific heat ratio of gas 

HI 


hi 

height at the upstream end of a non- 
circular section 

H2 


h2 

height at the downstream end of a non- 
circular section 

IFLAG 



parameter indicating the sequence number 
assigned to the fan section 

IPL0T 



decision parameter for selecting which 
(if any) plots are to be plotted 

IPRINT 



decision parameter for requesting or 
omitting output of summary character- 
istics page 

ISEC 



section type-description code 

ISEQ 



input section sequence number 

ISHAPl 



section upstream-end cross-sectional 
shape code 

ISHAP2 



section downstream-end cross-sectional 
shape code 

ITITLA 



assumed case-title array in the event the 
title card is omitted 

ITITLE 



input case-title array 

ITUNNL 



wind tunnel circuit-type code 

ITYPE 



code for type of output format required 
for printing section information 

lU 



units-of-measure type code 

LINEMX 



maximum number of output lines per page 

LMTSEC 



limit for maximum number of sections in 


any given case 
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F0RTRAN name 
and/or (title) 

Engineering 

system 

Description 

MCHECK 


master-card input-requirement checking 
code array 

MDATA 


master-card integer input data array 

MF0RMT 


master-card output format array 

MWRITE 


master-card output array containing data 
and/or annotation (s) 

N 

N 

section assigned sequence number for 
order of occurrence in circuit 

NCHECK 


section-card input-requirement checking 
code array 

NDATA 


section-card integer input data array 

NF0RMT 


section-card output format array 

NN 


section type number for printing proper 
section title 

NWRITE 


section-card output array containing data 
and/or annotation (s) 

P 


input tunnel total (stagnation) pressure, 
standard atmospheres 

PA 


input atmospheric (barometric) pressure, 
standard atmospheres 

PATH 
(P ATM) 

^■^ATM 

atmospheric (barometric) pressure, 
N/m^ (Ib/ft^) 

PHI 


corner flow turning angle, deg 

PI 

■n 

ratio of the area of a circle to the 
square of its radius 

PRSTY 


porosity of certain non-solid flow 
obstructions: AL/A 

PT 

Pt 

tunnel total (stagnation) pressure, N/m^ 
(lb/ft2) 

PWRI 


decision parameter for requesting or 

omitting the matching of power consump- 
tion with given input value 
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F0RTRAN name 
and/or (title) 

PWRIP 


Engineering 

system 


Description 


power required to be input to flow in 
order to drive wind tunnel at specified 
speed, W (hp) 


PWRMCH total power value for which the maximum 

test section velocity is to be deter- 
mined (if requested), W (hp) 


PWR0P 

^REQUIRED 

total fan motor output power required to 
drive wind tunnel at specified speed, 

W (hp) 

QO 

‘lo 

test section upstream-end dynamic 
p V ^ 

pressure: » N/m^ (Ib/ft^) 

R 

R 

gas constant, m^/sec^ ®K (ft^/sec^ ®R) 

RH0S 

P 

local static density, N sec^/m** 
(lb sec^/ft**) 

RH0SF 

Pf 

static density at the fans, N sec^/m‘* 
(lb sec^/ft^) 

RH0SO 

Po 

static density at upstream end of test 
section, N sec^/m^ (lb sec^/ft**) 

RH0T 

Pj 

density computed for total (stagnation) 
conditions, N sec^/m‘* (lb sec^/ft**) 

RN 

RN 

Reynolds number : 

RNREF 

RNref 

reference Reynolds number at which turn- 
ing vane 90°-loss constant, EKTV90, was 
determined 

RNV 


Reynolds number for turning vanes based 

pVcv 

on vane chord: ^ 

RUFNES 

A 

surface roughness in honeycomb cells, 
m (ft) 

(RUFNES) 


surface roughness in honeycomb cells, 
10'^ m (10-® ft) 

SEKO 


section total pressure loss array (refer- 
enced to test section conditions) used 
in summary calculations 
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F0RTRAN name 
and/or (title) 

Engineering 

symbol 

Description 

SEL 


section centerline length array used in 
summary calculations, m (ft) 

SERR0R 


section input error occurrence code 

SLAMDA 


friction coefficient for smooth pipes 

SLMDAE 


calculated friction coefficient in test 
section at the requested power-matching 
condition 

SLMDAl 


friction coefficient at section upstream 
end 

SLMDA2 


friction coefficient at section down- 
stream end 

SLR 

6 

s 

diffuser side length ratio: ratio of 

change in height to change in width 
from upstream to downstream end, or its 
inverse, whichever is less than or 
equal to unity 

SMACH 


section downstream-end Mach number array 
used in summary calculations 

S0A 

(S/AL) 


ratio of flow-obstruction drag area to 
local flow area 

SSUMEL 


summation array of total centerline 

length from start of circuit to end of 
local section 

SSUMKO 


summation array of total pressure losses 
from start of circuit to end of local 
section 

SUMEKO 

N 

2 Kq. 

i=l 

summation of all section total pressure 
losses referenced to test section 
conditions 

SUMEL 

N 

E Li 

i=l 

summation of all section centerline 
lengths (total circuit flow length), 
m (ft) 

T 


tunnel total (stagnation) temperature, 
°C (°F) 

TH 

0 

diffuser half-angle, rad 
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F0RTRAN name Engineering 


and/or (title) 


symbol 

Description 

TH2 

(2 THETA) 

20 


diffuser equivalent cone angle, deg 

TLIST 



case-fatal error occurrence code 

TLISTI 



decision parameter for requesting or 
omitting tabulation of input data 

TRETRN 



decision parameter for requesting return 
for additional case or final 
termination 

TSBLKG 



test section blockage used for computa- 
tion of basic test section conditions, 
percent of test section cross-sectional 
area 

TT 

Tt 


tunnel total (stagnation) temperature, 

Cr) 

V 

V 


local flow velocity, m/sec (ft/sec) 

VOC 



calculated test section velocity at 
adjusted power level, m/sec (ft/sec) 

VOK 



test section flow velocity at input 
conditions, knots 

VI 



section upstream-end flow velocity, m/sec 
(ft/sec) 

V2 



section downstream-end flow velocity, 
m/sec (ft/sec) 

W1 

wi 


width of upstream end of non-circular 
section, m (ft) 

W2 

W2 


width of downstream end of non-circular 


section, m (ft) 
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* 1 0 Q. ,H £ A P tS-ii-ULfl 0 ) HI I.LE 

1PA6E s ) 

H-.M—O. 

SUMfKO • 0,0 



TU8T « 0,0 

11*^' tITlE Card V*|.10ITV check 

IF (ITITLfd) ,E0, ITIo) CO WToF 

Tt»T ■-»<>■. 

HHlTe(6,P500) 

wpTTEf6.95oa< ITTTLf 


!!w!^5* TLlNi^gL’^l8HB CONTROL DATA CARD ObERATIQN I 

*!oI BEADf5,700n ISEUtlTUNM i «IUtI i.H A pI,I8MAb 1« (PATAf 
C 

^ HaSTEb CiaD validity CHE CK 

C.. .. . 


^.S-1 A-ft AIIU.1K I » V A aHA *-W fc-liJBl 

ir (TUI8T ,6T. -5,5) WRITE( 6 . 9500 ) 

.IF (T U s T ,S t» « ? t5 l TU I aT ■ «3 »-. 

h«IT£( 0,9505) 18EQ,iTUNNC,IU,I8HAPl,ISMAP2,(0ATA(I),I»3,l0) 


c !!!!ftAsUB C*bD EbcOb.CHECK aHQ Bah IC PApAM EtEb definiti 

*Io2 CA| , | DATACKM) 

IPA6E ■ t 

^ ■ ft 

8UHEK0 ■ 0* 0 

lUMEL ■ 0,0 

IF (TLI8T ,LT, -2.5) WRlTE(fc,900«) 

IF fTlIfT At. ■2«5 .AHQ, lu ,Eq . I) HrITE f6i9Qa&l. 

IF (TLI8T ,I.T. .2,5 ,AN0, lU ,EO, 2) HRITE(0,9o06) 

IF tTLIST At. ■2.5) GO tO 2Q0 

U.I8T • 0,0 

P ■ OATAtB) 

IF (P ,LT, 1,E»A) P ■ I, 

T • Data 


pA 9 OATA(IO) 


f t 9 • • 

AfilC I 




leRPQRM 


£At£ fl. 


PM 1510 

Imte r national system of U hiiTi ( gn eii.as 2 o- 

c.. PM 1510 

E?lLtAI.D.-l l,7flQ7E>S WU15A0 - 

PT ■ P*101325, PM 1550 

PA.1_M t PA* 101185, PM_lSW- 

R ■ 2»P',79 PM 1570 

IT_P T»271,15 PM 1580 

EMot a EMgsTD*(TT/288,0)**,76 PM iSPO 

PXH^CH t aATAf7lPl,l-6 pm lapp 

SO TO lo« PM 1*10 

PM 1*20. 

C., U,8, customary units pm 1*10 

PM 1**0 

lOl EMyiTD ■ 1,719E«7 PM 1*50 

PT P Pa2ll*,217 PM -1**0 

pATM a pA«211*,2l7 PM 1*70 

PM l * Stt 

TT • T+aSP,* PM 1**0 

Sh ut a 6MUSTD»{TT/518.<Ha*.7* 

PWHMCH a OATA(T)al,El PM 1710 

PM 1720- 

C,,,,C£'^ERAL«POrH DIHINSIONAI. PARAMETERS PM 1710 

C.I. PM 17 *0 

le« AT a 80RT(6*P*TTJ PM 17|0 

RMpT * P T / -X apTT ) PM 17*0- 

PM 1770 

C....TESN8ECTIOM MACH NUMBER mEmTOMS METMoO ITERATION PM 17 *0 

C,,, PM 17*0 

VLQ-.l ■ &A^AX*J PM- IBOO 

lo5 IP (lU «EQ« 1) VOX • V0*l«*«38 PM 1810 

IP-I.I.U iEia, 21 yOtt ■ V0* »5* l ft 8 PM_-1810_ 

EMO a VO/AT Pm 1810 

1-0* A.8Q a AT/SoRTM.Ar6»l.l/2.aEMOa*2) PiLUlo,. 

EM a VO/ASO PM 1850 

IP (ABStEMQ.EM./EMQ ,LT, 1,E.«) 60 TO 107 PM l8*0 

EMo a EM PM 1870 

GO TO 1ft* PM 1880 

1q 7 EMo a EM PM 1890 

C«t.im... PM 1*00 

C,,..MACH«NuHBCp.OEPENDENT PArAMEtEr* PM 1*10 

-C-M-, P M -1920 

EHOSQ a 1.4(6*l.)/2.*EM0*a2 PM 1*30 

ftMOaO a pHOT/EMOgO»A(l,/(6altl) PM 1980 

00 a RM080 *Vo*a2/2, PM 1*50 

C.»..,ti,.. PM. 1**0- 

C.l.TEsT section tHrOaT size and tHrOaT*ArEA«DEPENDENT parameters pm 1*70 

-X.,,, PM 1*B0- 

TS8UK6 a DATA(5) PM 1*90 

AC ■ PAR(DATA(H,DATA(a))a{l,«TSBLKG/100Q PM .2000- 
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g I ! o o «-i li 


£f££MM. 


RAfiC. —S 


IF (ISHiPl .eO, 3) aO I FAFO(DATA(3)iDATA(«))*a»»T8BLKO/lO0,) 
I F tlBH *p I ,FM, n AO 1 FAerDATA(Ani>^lt»T8 iLK0/i00, > 

DHo ■ F0M»(DATA(3),0*TA(a),A0) 

I ^ ( 1 8 H A p i..^LQ«..3l.atULX-£C HftoXDi-T 4.tXl4 QAUUUAQ 

IF (1SHAP2 ,EQ, 1) DHO • FOHC(DATA(A),AO) 

A8i*R..i 

HNoC ■ RHQSO*VO*AO/eMUT 


C»t.ii I* **»... 

C.i.iTUNNEL.TYPI NaHE ASSlONHeNT 


I * • • 


lp» 

lo** 


110 


(ITUNNI.'' 
108 I_I_ 


T ,ANO, ITUNNll .1,8* i) ®0 fo 108 


iCE| 

1,-3.. 
LT1T|.A(9) 


IF 

DQ 

LTlTLf(l) 

CQMIwUl 

GO TU n2 

U T I TLLlU-. ■ LIltLRLlI. 

UTITLCCA) I 
LTi.Tlf.m ! 

LTITIE(8) I 

LTlTLE.t3J ! 

IF (Itunnl 
I.T lTL.Eai.1 

CTITLEfS^ I LTIT|.A(2) 

GO To. 112 

IF fITUNNL ,NC, 2 ,ANO, ITUNNU ,N|, 5) 00 TO 111 


LT1 TLA( 8 i 

.4^Tl.IiiLtlil 

LTlUAdO) 

..wiii^Aasi 

I ,AND, itunnl ,NE, A) 
LllILAtU 


00 TO lie 


LnTLFtlL_t-LJ-liaxaJ 

UTITLE(2} « UTlTLAftt, 

GO TQ il2 

111 UTlTLE(l) B LTITUA(5) 

LTITLIiil • J.T lKA_jji.» 

112 IF (ITUNNL ,6T, A) 00 TO ill 

lTlTLft 8)..»_tnT4,AL_3J 

LTITLE(6) b LTITLA(10i 

. GO lO ll<» 

lU LTlTLEfS) B ITITLA(II) 

.. LTITLEIAI i^LlllLAlii) 



. I . , OUTpUL of Case bA G E t title , HEaOINCs and. OEFINIM JLAjJlMETlRt_ 


' 1. 1 « H R ITS ( t,, 9_05 OJ ..1 Tilif aUPAGIxLIUU. - 

IF (lU ,E0, 1) HRITE(6,9002) PA,PATH,P,PT,T,TT,V0,V0 k,O0 

if (lU .EQ. 2) WRjTE(8i30 03) PA.pAIMjP,PTj Tj11^YO^_¥ OK .^ 

HRITE(ft,900A) 

lFUlU.EQ.JLl_J^Rll£iAj-3JlP^ 

IF flu ,EQ, 2 ) wRITE(6,9006) 

LlNtCI.l_l3: 


IlR^TURN'pVTHrEREylQUS-pORTlON ONLY FOR BEGINNING CALCULATIONS FOR 
,,nEw wind tunnel circuit 


PH 2010 

PM 1020 

PH 2030 
...PH 10 AO 
PH 2090 
. _PM 20A0 
PH 2070 
PH 2080 
Pm 2090 
.PM 2100 
PH 2110 
Pm.2120 
PH 2180 
Pm IIAO 
PH 2190 
Pm .2190 
PM 2170 
PH 2180 
PH 2130 
PM 2200 
PH 22i0 
_PM 2220 
PH 2230 
. 8M-224Q 
PH 2290 
P.M 2290 
PM 2270 
,J»M 1280 
PH 2290 

l»M .1300 

PH 2310 
._PH 2320 

PM 2830 

_ELH.13A0 

PH asso 

PH_1390 
PH 2370 
. PM_2380 

Ph 2390 
.. PM 2A00 

PM lAiO 

PM 2A20 

PM 2430 

PM 2440 

PM 2490 
. - PM 2490 
PM 2470 
A PM 2A60 
PH 2490 
Ph .2900 
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PrKFOHM- 


PAtte ft 


Q************************************************* ********************** 

c..- «aXI«-Um number 0F_.8£CTX0N8 per CftSE-ULMIT.CilCCK ,AN0 riEaSAttE 

C I * 

8<)0 IP th ,LT» UMSEC) QO-IO~i0l- 

IP (LInEcT ,CE. (LINEMx-7)) IPaSE a IPAGEftl 

IP (LINECT ,6E, (UINEMF-71). KRITE.(ft,9a4U-tmUt*lPAfiE-. 

*«PlTE(ft,P505) LHT8EC 

“<NITE(6,9501) 

call DAtaCK(5) 

TLI8T i ,<|, 

'>: • 0 

IPA6E ■ IPaCE*.1 

N»ITE(ft,900t) ITITLE,IPAOE 

wSilTEtftrPOORJ 

IF (lu tFf;, n w«1tE(9»’008) 

IP (lU ,EO, 2j >"«ITEtft,9006.> 

LINECT a 0 


c:;ieefifir6:fx*55?srfioN8 

* 1 0 1 t(E A D ( 5 * 1 0 6 2 ) 1 SE fjT rate ri 8 R A P f» 1 8H-rpIV(DAT A ( I ) » I aj ,1ft) 

IP (ISEO ,(}E, 0) 60. TO 2o«_ - 

TLI8T a .2, 

IP (N Qj (JO TO jQ2 

wf>lTE{ft»9S0ft) 

LlMtT.«.LINEC,tft5 

60 TO 102 

idi »<«1TE(6,95o9} 

LINECT a LINECTtS 


C 

C 

c 

c 


op, J.0 J .1J 

ITITLE(I) 

. ^Oi.CPNUNyL . 


1, 21- 

I ITITLA(I) 


• • 


PIFJN] . 

termination Card 


l _POR _ QMITTED 


IP«INt a 1 

IPLOT-a.i! 

TlISTI a 1. 

PNRI.J Q 0 

TRETRn a •!, 


60 . T O, £002 — 

29*1 IP (18EC Ni, 0 .AND, TLI8T ,LT, .a,?) 60 TO 2o5 

IF (lSEC^Nt*_SL-..JLNCA.J[mTL^.a^^ 6Q T0 -2fllft 

c!!’** 0 |FINItI0n.J?,P_.iEbM|Nat 1P.N iAa-K»He8ui87_£A&AliIT|R 

c • * 

IPRInT a I8MAP2 


Pm 2910 
Pn 2820 
Pm 28J0 
PM 2S<10 
PM 2550 
Pm 2Sfto 
PM 2570 
PM 2580 
PM 2590 
pH 2ftOQ 
PM 2ftl0 
Pm 2620 
Pm 2ft3o 
Pm 26«Q 
PH 2fcSC 
PM 2660 
pM 2670 
PM 2680 
Pm 2690 
PM 2700 
PM 2710 
PM 2720 
PH 2730 
PM 2760 
PM 2750 
PM. 2760 
pH 2770 
. En .2780 
PM 2790 
PM ,2flaQ 
PM 2810 
.PM... 2620 
PM 2830 
. PM. 2640 
PM 2850 
_PM_2660 
PM 2870 
.PM -2 680 
pM 2890 
PM 2900 
pM 29t0 
PM. 2920 
PM 2930 
PM 29 <(o 
Pm 2950 
...P.M, 2960 
PM 2970 
PM 2980 
PM 2990 
. PM 3000 



ooo r»r>r» j looo 


.StfifJtWJL 


JUfiE- ... I 


IPLOT m IFIX(ab 8 ( 0 aT*( 1 ))) 

_,.TLl8n t ASSfOmiil) 

PwPI ■ AB8(0ATA(J)) 

... TRE.IRN -* . AeaLaA,TiIJU 

50 TO 2(*o2 


PH iOlO 

.J»JK-la 2 C! 

PH S 030 


* * 


Case-skip section check 


*205 


• t- 

20b 


N ■ N +1 

CALL OATACKtil. 

UO To 2 00 


PH Soso 
_P.H_JDAa 
PH 1070 
-PH- 108 Q. 

PM Sooo 
-.EM.-J 100 
3110 


PH 


KiOHviU section check 


N’ B N ♦ 1 

CALL OATACKtal.. 

IF fSFPROR ,(,Tf 

EL ■ OATAtSx 

SELfN) m EL 
SUHfL B SUHEL*EL 


1 «) 60 TO 200 


IlIIsECn^N .UPSTPEaH 1 -Hii-(IN.LET 0 . 1 . 

BlKASE ■ .OAtAXIOJ.. 

IF ( 1 $EC . 6 E, bO ,AnO, 0 aTA( 2 ) ,GT« I,E- 6 ) 8 |,KA 6 E 
IP fISEC ,LE, b .ANQ. Aft 8 (T 8 BLKG«ftA-Uila)j_.«(Lt 

1 ''RTTtfb.P 5 C 7 ) BlkAGE,T 8 BlkG 
IF ft SEC , LE b_» A 10 ^AMXlSaLK 6 .CLATAUO.J )-_.Al^_UEaii 
1 (.InECT a i 1 NECT *3 
IF ClSHiPi .Jig* , 2 v iO-IiJ- EOT 
IF (ISHApl ,E 0 , 3 ) Go TO 208 

Pi a D*TA{Ai 

A 1 « FAC(Ol) 

IF ( 1 8 E.C.. ^ £Q >..2 .gp.-_ 18 .e£..^ 19 ^ - 50 -. 0 ft 4 

1 |SEC (Eq, Ob) A 1 a A 1 « ( 1 ,-B|,k AGE/I 00 , ) 
t>Hi a FOHCf Dl f Al) 

GO TO 2 09 

2 0 7 H 1 a_D A1.A4J1.1 

a 0 ATA(fl 5 

A 1 t_ E ABfHl ,wn 



PH siso 

j»i 4 _ 4 Uq. 
PH SISO 
.J»M_S 180 

PM 5170 

Sl80 
Sl90 
S 200 . 
S210 
5220 
S 2 S 0 
S 200 
S 2 S 0 
5200 
5270 
5280 
S 290 
SSOfl- 
S 3 I 0 
5520 

5530 
iJM. 


PM 


55 S 0 

JJ 80 


jP fjSEC ,EO, 2 ,OR, 


ISEC .Eg. tt .OR, i 8 EC 
Al*{l,-8LKA6f/10(L,) 


208 


I I 8 EC ,Eg, 96 ) A 1 a 
OHi ■ FOH(?(Hi,i,i.AI) 

5 tLJ.Q_.ia 9 

HJ ■ DATA( 3 ) 

“tl L JiAlLUL} — 

Al a FAEofHl.rtl) 

IF ns&£_ ^Q.^._i,.,.o,R,_iSEc .es, 9 
1 ISEC ,EO, 9 b) Al a A 1 *( 1 .*BLKA 
DHi 1 FDHF-QCH.l.».aji,-A.i). 


,EO, SO .OR. 


PM 5570 
P-M 1380 
PM 5390 

.i»Jl_S«O.Q 
5410 

_lMo. 


PM 

_PH 


PM 

_PJ 1 


PM 

PM 


5450 

.l4flQ_ 

54 S 0 

-J 8 AQ- 


09 . 


ISfcC-,Eg. 8 j»_ 
U) 


^Qfi. 


PM 

-IM- 


5470 

JUftO 


PM 5490 
._PM_J 50 (L 


58 


gfcRFSSM.. 


PAttE a 


. • CW_0Q«ji4tgfc4M. ^ X Ll_._0fi._lMO - 2 J - Qiofit tfi Y . C0MPU,UllDllsia 

!<> IF,._aSHiP2._^EB^ ,,3J_SD. .TCL.211 

IF (IShAP? ,EO, 2) 50 TO 2I0 

02-. f. OAtAtTJ - 

A? • FAC(02) 

0^2 « 52 

SO TO iii 

.« DAUtaJ 

*<2 ■ OATAfTi 

- *2 t .FA.!jf H£4.1sU 

0H2 • FOHH(H2,w2,A25 

6i3- T0. It2 

11 H2 • DATA(6) 

'il. 1_ 0 A J.A tlJ 

A2 ■ FAF0(H2 ^2, 

0H2_ ■ _FpHFu(H2tMji,A.iJ 


1 . A t T.Efi.A T I_0*i3. 10. g| o>^E.IRV_9.gFJNlTloNS..5u.E..T& ■i?UUigL£ fi.ycTlM. 


iA_ Mp, _ i8e C 


12 IF fISEjC .L T. 60 .OR. ISE C .EQ. 8! 
EnoUCT a OATAd) 

iF dSEC .GE. 91 .Ai^O. tSEC .L E . ^ 
I fcNOOCT a 1, 

All ■ _A1 

Al 2 a A2 

*i 8 .*ll*|NP_UCl 

IF (ISEC .NE. 88) A2 « Al2*ENDUCT 

iAEc__>6e4-aj 

Fan s ection GEomEtPY 


GO TO 


P-^yi- ■ OATAfO) 

IF (ISEC ,NE, 91) GO TO fill 

OF A__M_8__0AU_tiLl 

OHi a OFAN.OHIJB 

QH2 a DHi 

Al a Pl*(OFAN*a2.0HiJB**2)/«,*(l,.flUFAGE/lOO,)*lNOUc'T 

_A2 a PI*f02a a2»OHlJ8*« 2)/4.* fl.«BLK A6e /IOO.)*EN DUC T 

00 TO 2)5 

1 3. IL lIS.le._,N E .^..9.2) GO TQ -.2JJ 

DFAN a OATAf?) 

0.M2 a_OFANaDHuH 

A2 a Pl«r6F'ANaa2*DHUBaa2)/«,aEN'6uCTa(I,*aUKAGE/lOO,) 

_ -*U-,S-A1/1.NPU£1 

AI2 • A2/CN0UCT 

...-6Q. T0_^l5 

l« OFAN a OATA(A) 

QHi a DFAN_DH Ua 


I 


# 


■ PI*(OFAN**2"OHUB**2)/fl«*(l«»8LKAGC^100,)*gNDUCT 
» 11 ■ Alye^^DUCT 


A20AO « A2/A0 

_ AR 1 Ag /AJ 

IF (AR ,LT, U) AR • Ai/A 8 

CAUL SPEEDtAliAMACa^ V -l 

CAUL FRICTN(0HI,AI,AMACH,8UAM0A) 

_-Ah A CHi .-1 A lhACtj 

VI • V 

8 LMOA 1 ■ 8 LAMQ A 

CALL 8 PECr)(A 2 ,AHACH,V) 

CALI FR!CTn/DH2.a2.AHACH.S^ . A mDA1 

AMaCHj a AMACH 

IM ACHfN^ .j Af iACtta 

V 2 a V 

— 8 L>!iflA 2 ■ 8U AM0 A 

IP (I 8 EC ,NE, 5 AND, ISEC ,NE, 4 ,ANO, ISEC ,NE, «0 ,ANO, 

1 I 8 EC .mE, 64 .AmD. fiatC.^siE. 91 -.OR. fAa/Algi ^.1 >Ul.^-L^iH- 

2 ,ANO, I 8 EC ,NE, 94, GO TO 224 

* 0 IHnItION op OIPPUSER.ONUT Pa«A^TER 8 “ 

TH2 a PTH2fA2,Al,EU) " 

IP nsec .GE. 60 .AND. iSEC^ NE^aSl I^H2_a p.TH2LAi2j.AU t£U.l 

TH a PTH(TH2) 

EkEXP a PATAflEl 

IP (EKEXP , 6 T, ®0 TO 224 

0ETau]^"^U8E’D NAT Ion OP ^IPPUSER EXPAN 8 lOIN U088 

parameter 

EKC a P E K C(tH2) 

IP fTM 2 ,UT, 3,) EKC a PEKC8(TM2) 

e 1^8 a FBK 8 (Th 2 ) " 

IP (TM2 ,LT. 3.1 EkS a PEK88fTH 21 

IP (TH 2 6 T, 10,) IKS a PEK 8 H(TH 2 ) 

Ek 2DR ■ PEk2DLtTH2J 

IP (TH2 ,UT, 3.) EK20R ■ FEK208(TH2) 

- tP (TM2 .gJE.^^4.1 Ek 20R ■ ..PgltiPUtrH^J 

IP (TH 2 ,GT, 10,) EK 2 DR a PEK2DM(TH2) 

EkC SAV a fEKC»EK81/ 2. 

IK20C a EK20R«eKC/EKS 

EKgPCa a fEK?0R*EK2 DC)/2, 

IP (ISHaPI ,N|, 1 ,OR, I8MAP2 ,Ng, 1) 60 TO 216 


60 




. — PAa£--io 


C -» t 


both ENoB circular 


SLR R I, 

EKHARg ■ Bur 


EKAOD p CKgDC 
SO TO 221 


2|6 IF (I9HAP1 .NE, i .OR, I8HaP2 ,E0. 1) SO TO 21« 


Ci. 

C- t ■> ■ 


UPBTrEaH end only circular 


IF ((H2«0l) ,6 t* 0,0 .AND* (m2«DI) ,6t* 0,0) Go TO 2l7 

9J.B- 1 -04 0 

GO TO 2)8 

■liZ- SLR .P _tH2aO-l)/.(W2»01) 

IF (SLR *6T. 1.) SlR • 1./SlR 

8l8 If f I8 H AP2. .EB* 2)-£KaA.aE-il -EKCIaY 


Pw «5)0 
Pm 0S20 

PM RSJO 

PM asdo 
PM 4550 
-PM a 560 
PM 4570 
P.M 4580 
PM 4590 
.RM .4600 
PM 4610 
-PM. 4620- 
PM 4650 
. J»M.J1640. 

PM 4650 
-PM-4660- 


IF (IsHAd2 ,E0, 2) EKAOO p EK2DCs 
IF (IS HoPa ,E0, 5) EKBASE p EKC 
IF (18HAP2 ,EO, 5) EKADD p EK20C 
GO TO PPl 


PM 4670 
.EIL488Q 


2)9 IF (ISHAP2 ,Ni, 1) 60 TO 220 


PM 4690 
-Pm_47AQ- 
Pm 4710 

_Rm_4J!21)_ 


C t • 

X.A- 


OoNnSTREAm only circular 


iT 


8lR 

I F fl SHAP ) ,EO, P) EKBaSE p EKCSaV 


I F iisHAp! !pq! i! bkbasf^/f kc^^ 


c • • 

C t a 


IF (IsHAPl ,Eo. 5) EKAOD p EK2DC 

.a5--TQ 223 


C. 


BOTH ENpg NON,CIfiCuLAp 


ho IF ((H 2» H1) .FQ. rM2»wm SLR p 1. 
IF ((H2»HU ,Eo. (n2»wI)) so to 222 

IF (tH2«Ml^ ,GT, 0,0 ,aND, fw2,W l 

SLR P O’.O 
JtOJLO 222 - 


-TO. 82 1- _ 


221 SLR ■ (H2-H1)/(N2>M1) 

J£-IjJ.R aXT* I .* } SL R ■ I , / SL R. 


222 EKBAgE p EKC.Ay 
ERAOP p EKgOrS 


PH 4750 

PM--,4740- 

PM 4750 

PH. -4 760 

PM 4770 

PH .4780. 

PM 4790 

£11-4800 

PM 4810 

£ii_4920 

pM 4850 

Pm .4840 

PM 4850 

eM_1860 

pM 4870 

PJ4-..4840- 

PM 4890 
PJi-Jt40(L 


IF (ISHaPI ,EQ, 2 .AND, ISHAP2 ,EQ, 2) EKBASE a EK8 
IF_(I8 H aPi ,EQ, 2 ,AND^ 18HAP2 1 0.^-JLKADp p EK2PR 
iF (jSRApI ,EQ, 5 ,AnD, t3HAp2 ,E0, 5) EkBAsE p EKc 

LF_dl i- - - 


Pm 4910 

£^lJ^S^o. 


>, TS 

SHaPI .EQ. 5 .and. !SHAP2 .6Q. 3) EK ADD ■ EK 2PC 


225 EKExP P EKSAsE4(1,«8Lr)*(EKADD«EkBasE) 



paging-chick before section information output 


PM 4950 
-_£M_4?40. 

PM 4950 
._£ft_«S60 
PM 4970 
- PH...4«ao 
PM 4990 
-PM-iooo 


61 


pe r fo rm 


J»ifiE — U 


2E« IF (LIMECT .LT, (UINEMX-2)) 60 TO 225 
IPA6E ■ 1Pa6E4^ 

PM 5010 
PW Sflifi 

WRIT|{fc,900l) 1T1TLE,IPA6E 

WRITE ffc.OftnUl - 

PM 5030 
.JilJlflOO 

IF (lU .1(3, 1) WRITI(6,9005) 
IF /lU .PU. WRITE 

PM 5050 
P-K -5060 

LIniECT ■ 9 
r 


PM 5070 
_ .P.M .5060 

c,.,.8fcCTJCi>^-TvPE Branching 

r . . . 

PM 5090 
PX -5100 

225 IF (ISEC ,E0, 1) 

Go To lOtO 

GO TO 1010 

Pm 5110 
PM 5120 

IF flSEC IeoI 3) 
IF rISEr .EO. 

GO TO 1030 
GO TO 1030 

PM 5130 
P..M 5140 

jF (jBEC ^EO^ 5j 
IF rISEr .EO. 

GO To 1050 
GO TP 1050 

RM 5150 
PM. 5160 

IF (I8EC .eg, 10) Go To MOO 
IF rlBEr _EG 90\ GO TO l2oo 

PM 5170 
P.M 5ieo 

IF (isBC IeoI 30) 
IF (I8EC .EG. 32. 

Go TO 1300 
GO TO 13QQ. _ 

PM 5190 
PM 5200 

If ^18E^ *EQ* 

Go TO 1330 
GO To l3tto 

PM 5210 
-PH 5220 

iF dSEC ;Eq; OO) 
IF /ISEf. .EG. U5> 

Go To 1^00 

GO TO ltt50._ . .. - . .. 

PM 5230 
_PM 5240 

If (Is'EC ,EQ, OP) 
IF (I8EC .EG. 51. 

Go To i960 
GO TO l5)n 

PM 5250 
PM .5260 

iF (tSEC ,E0, 52) 
IF fjSEJ .E6. S3) 

Go To 1520 

GO TO 1530 

PM 5270 
PM 5280 

\F (j8iC ,EQi 50) 
IF fISEC ,EQ, 56. 

Go To 1540 

GO TO 1560 . . 

PM 5290 
PM 5300 

IF fistC ,E0, 57) 
IF (I8EC .EG. 61. 

GO To 1570 

GO TO 1610 __ . . 

Pm 5310 
pM 5320 

If ^IIec *EQ* 70^ 

Go To 1620 

GO TO 17QQ . . . . _ 

pM 5530 
- PM 5340 

tF (iSEC lE0r71) 
IF rlSEC .EQ. 72. 

Go To 1700 

GO T(J 1700 

PM 5350 
_..PH 5360 

IF (ISEC ,E0, 73) 
IF (ISEC .EG. 7<». 

GO TO 1730 

GO TO 1740 . .. ._ 

PM 5370 
PM. 5380 

iF (isec ;Eq', 75) 
IF (ISEC .EG. sai 

Go TP 1740 

GO TU 1640 ... „ _ _ 

PM 5390 

. .. PH .5400 

iF (jSEC ,EG. B5) 
IF fISEC .EG. 86. 

Go TO 1850 
GO TO i860 

pM 5410 
PM 5420 

IF (I8EC ,EQ, 87) 
IF (ISEC j.EO^ 91 . 

GO TO i860 

GO TO 1910 

PM 5430 
Pm 5440 

iF fiSEC ,EQ, 92) 
IF (ISEC .EG. 941 

GO TO 1920 

SQ_1I3 1940 - -- — - 

PH 5450 
_ PM.. 5460 

jF (iSEC ,E», 96) 

IF (I8EC ,,E0»..9J4. 

C.,.. ...... 

C.,..sINGLE»QuC7 -sEClLQ. 

GO Tn I960 

GO TQ l97C . 

PM 5470 
. PM. 5480 

NS-IOTAL-PrESSUrE WOssES - - 

PM 5490 

PM 5500 
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P twroRM 


- PaSE iZ 


C • f • 




PM 5510 

Ctfift 




Pm S 520 

C«t 

TEST sections 



PM 5530 

X 4 -. 



_ 

— PM .5540 


C»t Pm 5550 

£*4 CL 0 SED» COMETAhifABEA TEST. 8 ECUa»i PK 5560 


C 


C 

£. 

C 


I. 

C 

L 


1010 EK ■ <tLMDAl»EL/PHl 

EKq ■ PEK0(EK,A0,A1,EM0,aMACH1«CHO8Q,C) 

8EK0(Nl ■ tKO 

SUhEKo « SU'^EKo^EXO 

CALL OUTPUTMfl) 

LINECT I LINECTaJ 

t lStC «ES, 1) Cq T.Q. 2tfl- 

>10 nEL IN JHE TE 8T fiE£.I.10N - . 

ISeo i|OA ■ DAtA^BI 

CO ■ OATA(IS) 

EPb ■ l.»DATAf 161/100 ^ 

EK ■ C0*S0A*CP8 

EKo m FEKO^K.AO.Al.EPO.AHACHl.EHOSO.fii 

SEXo'eN—i SKMfNj'tEKO 

8U w£kO ■ S UWEkOaEkO 

CALL 0UTPUT(2»5) 

kliLlCt—i . UMC til 

GO TO 200 


CLOgED, DIFFUSING TEfT sECTlON 


PM 5570 

PM.5580 

PM 5590 

Em 5600 

PM 5610 

PH 5620 

PM 5610 

Pj«_S660 

P** 5650 

pJd-5660. 

pM 5670 

PJl-_5680 

PM 5690 

_JLM-S700 

pM 5710 

PE -5720 

PM $7j0 

PM- S7«0 

PM 5750 

PM .5760 

PM 5770 

PJl- 57 flQ 

pM 5790 


^ PM 5800 

1030 EK ■ (EKEKP«SLMDAl/f8,*8lN(TM))*(AP9l,)/(AR«l«))*((AP»l,)/AP)**2 PM 5810 


IF fE K .LTt fi|L MOAUEL/DHl)i EK ■ g LMDA UEL/O Ml p« 5820 

EKO ■ FEKO(EK,AO,*i,£HO,AMACHi,lM««, 6 ) Pm S 830 

- 8 £KOiN.i...i E.KQ p.M - .5860 

SUMEKo ■ 8 U'^EKo*EKo Pm 5850 

CALL 0 UTPUTfl , 2 ^ -PM_ 586 Q 

LINECT « UINECT*! Pm 5870 

I^ LIaEC .£«♦ 3 ) GO_jO_jLQO_ PM .5680 

60 TO lo2o PM 5890 


C t • 

_c. 


OPEN-tHrOaT TEST SECTION" 


foSO EK ■ ,0 b65aEL/RH1«,0053*(EU/dH1)**2 ‘ 

EKo 9 FEK0 ( EK.AQ.AlfEM0fAWACHl^eM 0tQf6> 

SEKOfN) ■ |K(| 

SU''<£KQ, a Su m eho^PKO 


TM? ■ FtH2(A2,AI,EL) 

CAtL-QUIP,il.Ti3j«) 

LInECT a LlMECT +1 
IF -tI 8 £.t- 


PM 5900 

Pm 5910 

PM. 5920 

PM 5930 
PM ,5940 

PM 5950 

P_M. 5960 

PM 5970 

PM 5960 

PM 5990 
Pm 6000 
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■PgHRflBtl 


...PjUiE -U 


GO TO 1020 

cir** »TRAISHT, NoN-OIFFU^lNOTum 

Co^laTA NT»AR EA P UC-T - 

C.. 

1100 NN > s 

EK • 8LHDA1*EL/0H1 

EKQ ■ PgK0fEK.40.AUEMQ.AHA CHl»EM08QAGl,- 

ITYPE ■ I 

GO TQ 2000 - - - - 

0 • • 

CoN T RA£J-Lq!i 

1^00 - - 

TMj m PtH2(*1,*2,EL) 

IK •_*32»iL?^DA2<^ELy DH2- - 

EKO ■ FEKOfEK,AO, A2,EM0 ,AMaCH2,EH08Q,G) 

- 

GO TO 2000 

C,, Cf.R^iERS AK-0 TURNS 

C,, 

c.. Cn^sTANfAREA CqRNEw •• TURI^InG YA^ES-QNLY 

i SoO NN > 7 - - - 

CHOi?0 » OAlAfO) 

RHI m AR8(DATAtll)) 

EKTVRo « OATAMZ) 

h^REP ■ l>AlA(ia)*ig,**b _ . _ 

IF (l)ATA(ia) ,LT, l.Ewb^ RNREP ■ »5*10,**6 

IF CEkIVRO ,LT* l.E-Ol LkTV’O 1 ,15 

IF fPHi ,LF., 10,1 * AKTVI (PHl)*EKTV90/,l5 

tF (PHjt ,Gt, 10, ) EkTV ■ FKTv2CPHi)*EKTy90/,15 _ 

RNV ■ RNOC*tHORO/Al*(l,+{C,l,)/2,*AHACHl**2)** 76 
fcK ■ EKtV*(2, + CA( OS10(RnrEF)/A, oGl0(HNV))**2,56)/J, 

IF fISEC ,EU, 30i <iO TO liol 

. - -- 

c,, CoNsTANt.arEA CdRnEr inItH TURNING vaNES anO walls 

c.. . -- - - 

NN a ft 

£K , EK + SL^0A1*EL/0H1 . 

Ucl EK(. a FEKOfEK,AO,*l,£NO,AMACHl,EM08Q,6) 

ITVp£,^ B_l .. . 

r,Ci TO 2000 

f. ,, - 

t., CO^STanT-AREa corner wITn walls AND WiThOUT TURMnG 


PM 6010 
... PM i 020 

PM 0010 

..PM ftOMO 

PM 6050 

. PM 6060 

PM 6070 
. . PM. 60 B 0 

PM 6090 

PM 6100 

Pm 6110 
Pm 6120 
PM 6110 

PM 6160 

Pm 6150 

PM 6160 

PM 6170 
PM 61 B 0 
Pm 6 l 90 
RM 6200 
PM 6210 

P.M 6220 
Pm 6250 
P.M 6260 

PM 6280 

Pm 6260 

Rm 6270 
PM 6280 
PM 6290 
PM 6 JC 0 
PM 6310 
PM 6 l 2 U 
Pm 6310 
PM 6160 
PM 6580 
PM 6360 
PM 6370 
PM 6360 
PM 6390 
PM 6600 

RM 6610 

PM 6620 
PM 6630 
PM 6"60 
Pm 6650 
PM 6660 
PM 6670 
Pm 6680 
VANES PM 6690 
PM 6500 
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PftPjJJJJM. 


9 A&E 


USO NN • « 

»H 1 ■ ABXDATA(im 

EKTE90 • OATA(tZ) 

lf-t£lLlI 9 Q EKTE 9 i) 

£KTE ■ FKTE(PHl)»|KTE90/i,«0 

CK ■ EKT E»8LH P AitEL/ DHl 

EKO ■ FEKO(EK,AO,AI,EMO,AMaCH1,EHOSQ,G) 
ITYFE ■ 1 

90 70 aooo 


C » t - 1 .1 — . ■ - — — , 

c.« diffusing Corner wIth turning vaneb and wauls 

iJ«0 NN « let ' ‘ ~ 

CilORfi_I_BilAi9.j - 

RMi ■ AB8(0 AtA{U)) 

EKTV90 ■ PATAfiai 

RNREF • OATA(1«)*10,**6 

IF tEK TVRQ , LT, EKTVRq • .15 

IF (OATAflA) ,IT, l.E-6) RnREF a .5*10. **6 

J*t4- JiLa -EITY. ■ . FK T V H PH I j *tK T V9 0/« 15 

IF (PHj ,GT, 50.) EkTV a FkT V 2 (PMi } *tKTv90/, 15 

!Hj- «-g.TH2iA8,_AbEu..._., 

RNV a RNOC*CHOPD/AU(l,*(S«l,)/2,*AHACHUa2)**,r6 

EkV a . 3 _ 

IF (TH2 ,0E, 21,5) EKV a eKV*,006*(TH2»2i ,5) 

EkI a EkV «((AR»1.)/Ar)*»2 

EKg a EKTV*2,/S, 

CiU. FHigTNiCHflRD4AljA,MACMl,SLAN0Aj 

EK a EKUEK2 

£'<9. -»...Fll9i|H.,.4i),Ai,E«0»i.N 

ITYPE a 5 

_ _ 60 TO 2 000 

— — J3lFf.UlJih; . ... 

c.. 

"c!' oTFFUSER 

t2'oo~ nn li~fi 

E*' ■. lEKEXPWSt«'lAl/(8,*9I.McTH))*(ARwl, )/iAR«l a) *£(*«•! »)/AR)*»2 
IF fEK {SUHOAl*EL/Owin £K a 8LHDAi*EL/0Hl 

EkQ a FeK0tej<a0j_*l|Ei»04AMACMl,EM08(S,Gj 

ITYPE a 5 

. 8U.. T0 iQQQ 


C * , 

C.. EXlT_KJNE.tL£_.EN6i!Si,FRQM FLOW DUHP._,. 


C • • 

1950 


NN a 12... 


EK a 2,*((|,+{G»l.)/2.*AMACHi**2)**(tf/(G«l.))-l,)/(G*AMACHl**2) 
I.A.O.a F£KO tlY,A.O# Al,eMo,AWACHl,EHQaCiS) . 


PM 9510 
Ph 9520 
P« 9 SJ 0 
PM. 95(10 
PH 9550 
PM 9590 
Pm 9570 
PM 9590 
pM 9590 
Pm 9600 
Pm 9910 
-PM 6620 
Pm 9950 
PM 99 a 0 
Pm 9950 
PM 9990 
PM 9670 
PH 9980 
Ph 9990 
PM 6700 

Pm 9710 
P.M 9720 

PH 9730 

PH 6740 
pM 6750 
pH.i 76 Q 
Ph 6770 
.PJlJj760 
pH 6790 
Ph 6800 

PH 6810 
PH 9820 
PH 9610 
PH 9040 
PH 9850 
Pm 9090 
pH 9070 
PH 9600 
PH 9890 
PM 6900 
Ph 6910 
PH 9920 
Pm 6930 
PM 6940 
pH 9950 
PM 9990 
pM 9970 
PH 9960 
PH 9990 
Ph 7000 
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_e£Ef.Q8!5.. 


J?AQE. -45 



ITYPc • 2 

GO jO PAOO - _ 

PM 7010 
7020 

c« • 

r . - 

iUUOEN EKPAAifitllN .. 

PM 7030 
PM -7040 

c.. 

i4fcn 

HN a 13 ... 

PM 7050 
PM 7060 


TH2 a 90, 

IK a t £AP« 1 .)/AR)**2 

PM 7070 
-PH 7080 


EkO a FEkO(Ck.AO; AI.EmO.AMAChI.EmOSQ.G) 

PM 7090 


TTVPE 8 t -- - -- 

. PM 7100 

e 

GO TO 2060 

PM 7110 

.-PM. 7120 

f • 1 

c.. 

f . - 

FLOW obstructions 

PM 7130 
PM 7140 

— »*-§_• 

C i i 



HONEYCOMi .TMDi._E.UlM 5 TP.a1.6HTEnERS 

pM 7150 

RM...716Q 


c « « 


• fc) 

..«_AIS!D,. 

,AMO. 


lu 

lU 


• 5.s«. 

.53. 


lx 

2) 


RuFjxES. 

RUFNE8 




.000032608 


r«oi T».j/iug. 

Cam SpEEPf A^j AMA.CHiV) 


■ «H0f/Tl,+‘(G«lo/2,*AHACH**2j**(l,/(C.l,j) 

• E«UT/ tUtX45i-i)->?2,**MACH**2)*R.7.6 

■ fcMIJyRHOS 


c • « 

C « t - 


NN ■ la „ __ 

ELODH ■ 0ATA(8) 

PRSTY, a PAIAilQ)_ 

uUANiEs I OATA(1«)*10,**( 

IF (PAlAtlRJ 
IF [OATAfM, ,LT, I.E.6 

F El/ElOP^-- - - 

AL » FRSTy.A j/lOO' 

E*<u 

E'^U _ - . 

Rn * V*RUFNES/DHL/tN,y - 

SIAHDA ■ ,375*f R uFNE8/DMU**,«/RN«<»,1 

IF (Rn .GT. 275,) 5, AmDA ■ ,2t«*(RgFNE8/PRi.)**.^ — 

EK • 3UAHOA*tEUOC.Ht3,5*(lOO,/PRSTVj**2*aOO,/PH8TV»l.)*Ai2 
EKi; a FEK0£EK,Ap,Al,tM.(LiAPACHl,eM08Q»8) 

ITYPE a t 

GO TU loao- - - — --- 

AlafPlL. TWICK.._F4.0 w .STBlififlTENEBS 


NN » 15 

ELOOH a OAT*f«) 

PR8TY a PAJAL1P4 

IF (EL6^H ,LE, 0,0) 


,3tEw- ■ — 


2. 


ELOOH a 

ElC 

ELD a EL.ELC 
*L • PH5TY/100,.*A1 
Afv e A2/4L 
I'f^L a EL/ELOOH 
CALL SPEEO(aL» aMaCH,V) 

Call FrICTs;(PHL,aL,AmACH,8LAmOA) 


PM 7180 
PM 71Ro 
-PM 7200 
PM 7210 
J»M 7220 
PM 7230 
PM 7240 
PM 7250 
PM .7260 
PM 7270 
.PM 7280 
PM 7290 
.PM 7300 
PM 7510 
PM 7320 
PM 7530 
. Pm 7340 
PM 73S0 
-Pm 7360 
PM 7370 
. PM 7380 
PH 7390 
PM 7400 
PM 7410 
PM 7420 
PM 7430 
PM 7440 
PM 7450 
PM 7460 
PM 7470 
PM .7480 
pM 7«90 
PM 7500 
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Pace. 16 


„ PfREURr 


EK^NTR s *32«e(.C/0Ht*SLAH0* 

tH? ■ pTHgfAg,4L.ELDl 

CKV ■ ,S 

IP fTWa .Sr. 21. rKv m fkv»^ ooo.* ^ 

EKO • CKV*((*6«1,)/Ar)**2 

EK ■ EKCNTR 4 EK D 

EKo « FEKOfEK, A0,A1,EM0»A^'ACM1,EM0SQ,G) 

illPI -I_3 

60 TO 2000 


ctl pii^RATED pLaU WlT« 8H*RP»VbGEb'' oRiriVlS 

^ 11 - ■ 

ISSO ■ 16 



iTVpC ■ 2 

J50- _T_P i OfiCL 


C * • 

Cxs wOiiLN Ji'laH.. st«i£N 

c.. 

js ao NN j ^7_ 

OhEs'h i DAfA(P) 

PRSTY ■ OATAdoi 
EKME8H* OATAC12) 

If (fKMESM ,LT^_i.E-6) EKMESM * j,3_ _ 

Ri^O* • RWQT/(i.*(C«i,)/2.*AMACMUi*2)**(l,/{G»l,)) 

■ S«UIiiltY.L6-lj)/2.*AMACHl*t»2.3**,76 

£MU I EM(j/bHOS 

« Vi*t>‘'E8H/ENg 

tK 9 EKH63M*(l,-PR8Tv/l00,) + fl00,/pR8Tvl,)**2 

If..(RN .LT,_600._.AND, 

1 Ck .LT, (EK*(78,5**{lf»RN/J5R,)/100,*l,01))) 

1__1H f,E,K*£7a.5**(t*-RN/35«,)/lO0^tl,eij 

EKo ■ FEK0(eK,A0,AJ,ER0,A«ACMl,EH0so#G) 

LtYRE » » _ _ _ 

GO H 2000 


PM 7SlO 
P« 7520 
PM 7530 
Pm 7540 
Pm 7550 
-Pm 7560 
PM 7570 
PM 7580 
PM 7590 
Pm 7600 
PM 7610 
PM 7620 
PM 7630 
Pm 76a0 
PM 7650 
Pm 7660 
Pm 7670 
Pm 7680 
Pm 7690 
PM 7700 
PM 7710 
PM 7720 
Pm 7730 
PM 7740 

PM 7750 
PM 77b0 
Pm 7770 
PM 7780 
PH 7790 
PM 7600 
PM 7610 
Pm 7620 
PM 7830 
PM 7840 
PM 7650 
Pm 7660 
Pm 7870 


c.,^ Pm 7 eeo 

c.t internal structure (Drag ITEM($)) at upstream end op a pm 7890 

C.A SECTlpM PM 7900 

C,, PH 7910 

1560 NN s 18 p,., 70^0 

ENITEM ■ 0ATA(2) PH 7930 

- SOA. p, data (Si. PM 7940 

CD » 0AtA(13) Pm 7950 

- .E'»S t-l.T0ATA£l6)^.L9J}^ PH 7960 

IP (ABs(ENItEM) ,Lt. 1,E-6) ENItEM • 1, PM 7970 

EK « CD*SDA*EP8*ENIT£M Pm 7980 

EKO • PEKOfEK, A0,A1,EH0 ,AMaCh1,6M0SQ,6) Pm 7^90 

- ITYP 6 * 2 PM 8000 


pfWroaf^ 


PACiE 17 


60 TO 8000 

ct1 P 1 XEO , KNOwN.UOC *L»L08^ I TE^Tf THI^U^TRE rn’ENO OP 'V 

JL.-, 8Er.TI0*j 

C i » 

_IS10_NN_1_13 - 

EK ■ DATA(t}} 

EKq ■ PEK0(EK.A0f AI.EMQ. a HaCHUEHQSQ.G? 

ITYPE • I 

CO TO 2000 . . 

..LOSSES 

C 1 1 • 

c,, • T R * J G hY 7 To N ^ 1 p p“u s’! no’ Ducts 

c!! OfiNsJANi-Aj-EA -PuC-ia - 

c • « 

LSlfl NN_r 2_Q _ 

£K ■ SL^OAUEL/DHI 

f'Q Gi 

ITYPt « 1 

GP T.Q 80.00 - 

C ( • 

, c.t - C0<^,TaACl.i9NS — 

c • • 

. .1620 • 21 - ■■ 

TH2 • FTH2f#Il,AI2,EL) 

Ek » - - - 

£K0 B fEK0(Ek»A0,A2,EM0«AMaCh2|Em0SU«g) 

iTYPE.i. ii,_ 

f,U TO 2000 

c|,* *' cornErs V^o Turks' 


c , , C.o N s T A N.T.- A.S E A.. C 0 r N E « s . . I U 8 N J N G. . V A E 8 0 N l. Y 

Oft 

1700 « 28 - - 



pM 80)0 
P^< fl020 
PN 8030 
PM 8000 
PM SOSO 
... PM S060 
PM 8070 
PM 8060 
PM 8090 
...... PM 810 0 

PM eito 

Pm 8120 

PM 8130 
PM 6100 
PM 81S0 

PM 6180 

pM 6170 
PM 6180 
PM 8190 
Pm 6200 
PM 8210 
. . PM 8220 
PM 8230 
PM 8200 
PM 8250 
PM 6280 
PM 8270 
Pm 8280 
PM 8290 
PM 9300 
PM 6310 
PM 6320 
PM 8330 
Pm 6i00 
Pm 8350 
PM 8380 
PM 8370 
PM 6360 
PMi 8390 

PM 8000 


CHORD ■ DATAcpj Pm eOin 

RHj ■ ABSCPAIAXIIJ] PM 8020 

EhTY 90 t OATAM 2 ) pn 8030 


RhreF B DATA,mj*10.,*p8 - PM 6000 

IF (6KTV90 ,UT, l.E-6) EKTV90 b ,15 Pm 8o50 

IP tOATAtlO) ,LT, l.t»8) RnREP ■ i5aIQ.*b 8._ pM 80oO 

IF fPHI .EE, 30,5 E*<TV a FKTv I ( PH I j pEKT 090/ , 1 5 Ph 8u70 

tPMj ,6Tt IO.J.EkTV ■ PkTv2(PMi)*£kTv90/,15 ..... P" 8080 

KNv B RMOC*ChOPD/AI*( 1 ,+(G«l ,)/2 ,bAM*CH1**25*p,78 Pm 8090 

i* ■ EKTV*(2, + tALOG10(HNH£F5/4|.OGlO(R^JV5)p*2,58)/3, Pm 8500 
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-ElftfJIfiB 


PAGfc l« 


IF (JlEC .EO, 70) 00 TO I70t 


PM GSIO 


- - — . PM 6520, 

C,, C0N8T*MT»A«Ea corners WITH TUrMhC v*NCS aND WALLS Pm 6530 

PM 6540 

NN 9 24 PM 6550 

EK ■ EK»6LHOAUEL/OH4 - - PM 6560 

IP (I8EC ,EQ, 7l) NN ■ 23 Pm 6570 

- - — - — - ■ - Pm 6 S 6 Q 

C»» CONsTaNT«ArEA corners WITH TURNING vaNES AND ONLY ONE Pm 65R0 

JLm aiDE.WALL Each pm 6600 

C(* P<^ 6610 

If. (ISEC ,Er, TU-EK_f_EKA3«/4, Pm 6620 

l70l EKO ■ PEKO(eKiA0»AlfEMO|AMACHl»EMOSO#Q) PM 6630 

IllPl_J_J - . PM 8640 

GO TO 2000 Pm 8650 


JLm PM 8660 

Ct» CQNsTAN'j.ArEa COrnEr wItH wALLs ANO wITmOUT TURNING VANES PM 8670 

PM 868 C 

J730 NN a 25 PM 8690 

PHI t_A 8 JiDATAQi _ PM 8700 

Ek'T£P0 a DATA (12) PM 8710 

IP (EKTE 90 .LT. l,E» 6 ) EKYE90 a 1,80 _ Pm 8720 

EKTt a PKtEfPHI)*EKTE90/l,80 PM 8750 

Ew a EmTEaSLmDaUEi^/Ohi pM 8740 

■■ ENO a PEK!)(EK.A 0 »Ar, 6 Mo,AMACHl,EM 08 (a»G) Pm 8750 

ITVPE a 1 PM 8760 

GO 70 2000 PH 8770 

— - - — — - . — ■ Pm 8780 

C,, OiPPUsiNS COrNErS with TUrNINO VANE 5 and WALLS Pm 8790 


X.._ - - PM 8800 

1740 NN a 2a PM 881G 


CHORD a 0 aTAI9_^ __ _ _ 

PmJ a ABS(0ATA(11)) 

EKTV9Q a pAT.A.(L2J 

RnrEF a OAT*(14)*10 ,*a6 

if (0atA( 14) ,Lti _Iif"6i RjjREP « •S*J0tJt*6 

IF (EMTV90 ,LT, 1,E«6) EKTV90 a ,15 
- iF fPHj ,LEi 30^_EK,TX_»_.FK.TVlIP_HTjAtKTy9J/,lS 

Ip (PhI ,GT, sc.) EWTV a FKTV2(PHI)*EWTV90/,15 

THg a PTH2(A2,AltlLl - 

R^V a RN0C,CH0RD/M*(l,*(G«l,w2.*AMACHU*2,a*,76 

Ekv •..3 

IF (TH2 ,Ge, 21,5) EKV a EKVt,006W(TM2-21,5) 

Ek.1 • EkY*(1AR»1,1/AR1,**2 

EK 2 a EKTV*2 /3, 

CalL.fr I ClN£,tH(iRDjA.l*,AMACHiASLAM . 

C,. 

C,« diffusing cOrNErS wItH turning vANEs and only one 

r,, SIDE-WALL EACH 


.PM 8820 
pM 8830 
. PM 8840 
Pm 8850 
_PH„886Q 

Pm 8870 

PM._8860. 

PM 8890 
. .-.PM...8900 
PM 89io 
PM 8920 
PM 89j0 

PM..8940 

Pm 8950 

PM. 8960 

Pm 8970 
PH 8980 
PM 8990 
PM 9000 


o o o r> jn o 




—PAfiE . 4» 


IF (I8EC tEO, 7«) £«! i EKI-EU*8 L*Moa/OHI/«, 

IF tlsEC .EQ, 75^ tiN « .27 

EK a EKiaEKE 

gKO a rEKO ^eK.AQ.tl.EMO.iMACHl.EMOBQ.fi^ 

ITYFE ■ J 


BIFFUSErS 


ihuo NN a 28 


) ) a t A P_a 1 , ) ^All 4 E.t 4 - 

) Ek • 8Lh0M*El/0H1 


iF (tK ,LT, (8 uMDA1*Ei./DHI)) Ek • 8U>»0Al*Ei./0Hl 

EKO a FEKQfEK.AQ.Al^EMQ.A MAC^liE>*QSiLiiLl 

ITYPE • 3 

2000 

C * * * * * uaNEO DlFFuaEB^ 

.5fs..»_2a 

EKv a ,3 

T,H2..a,_FlHitA4^.Al^iuV -■• 

IF fTH2 ,GE, 21,5) EKV a t'KVt,008*{TH2«2l ,5) 

a..t:KVtiiAR.v,,,^A_R)a.a.2. - 

E^o a FEK0rEK«A0,Al,EM0,AMACHl,EMO8Q,6) 

ITYFE-.?-- 3 

GU tO 2000 

C,, .5MQaE.^£l?-JLN.8jON f Pq^1 MULTIPLE. O.UcT-I .TQ .8ING.Lt UUcT 

C • t 

1860 ■ 30 - - ... 

Tm 2 a 90 

a ((*R-Ul/AiiJ**.2 

£><0 a FEK0f£K,A0,Al,EHo,AMACHl,EM08G,G) 

IF . f IS.EC ,EQ^__ai)._NN..a.ii._ 

ITyPE a 3 

60 TO 2000 . 


_ __OPIV|-FaN 3Y,8.T.£M. .. . .. - PM 9900 

PM 9910 

■ _ PM 9920 

. . Fan annular DUCT( 8 ) WITH MOTOR. SUPPORT STHUTfS) PM 9930 

. . . PM 9990 

t«l 0 N ^ « J 2 '" PM 9950 

ENITE.M 1 . 0 ATAt 2 J_ . - . - FN 9960 

jHA a DAtA(«) F^ ^PTO 

to a OAfAfiS) - Fm 9980 

ETAFAN a OATAflS) FM 9990 

£PS a l,tOATA( 16 )/lO 0 , PM 9500 
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o r> o o rin«-»n ’ oooo 




.. PAiiE .20 


r«2( 


IF feNlTE*' .UT, J,E»6i ENITE” » 1* 

.. IP . CFTAF*N_,LT»-i >fc»fcl-£T-*-£i-N • 100. 

E»'F a At'iCMl 

. ,£‘JP,AN i..£NQU£JL 

IFL.A& a N 

«H 0 SF 1 . RHflJ/tl,*t.G-.l>.)/2»*.A.MACHU*2)*atl./tGaUi4 

tKO a SI.>'0AU£U/DH1 

-IF (tA2/Al.l.4_.6T.. J.E.fej. EKD a UKEXP*5WM5AlZ.CfiL.*SmiilJJ_t_ 

1 (AR+l,)/(AP*l,))*fCAR»l,)/*R)**2 

IF f£»<D..LT._.t8J.MPAl*EL/DHni EKO .« 81.HuAl*EL/QKt 

EKSTrT * C0*50A*EPs*£^lTt*^ 

£•' s eK!?+6K§TR.T ...... . 

£!<C a FEK0fEK,A0,Al,EM0,AMACHl,EMO3Q,G) 

. ITYPE 1.1 

IF ((a^/ai,1,) ,GT, 1,6-6) ITYPE ■ 3 

-,.GQ TO 2Q(I0.„,. - 

Fan CONTiiA.CiION(s) tO ANNuLAb PuC.T.t$J...M.TH....PC)TaR-JUPP_QP.T, 

STphT(S) 

I A'N 8 33 ' 

THg a FTh2f AIl,Al2t£U) _ _ _ _ 

6*^ « ,32»sUH6AUtL/OM2 

6«0 3 FP.kO(Ek,AO, A2^EM0.AMACH2,EM0Sa,0) _ 

ITYPE a u 

GO TO 2000 . - .. . . 


t • 

• • 

• • 

• • 

l9ao 


Fan PI^FUSER(91.FR0M jygNULAR 0MCTC8Jj> 
CokE^SHapED ce^TEpBo^y 


.1 4i £ms^ 




3« 


6 • 1 6 5 1 X P * SLM 0 A| / 19 , . S X N t T H J ) a ( A R A 1 , ) /4 4 R O * U.A .R.-LUi^Mt.*i 
iF {6 k ,tT. (St'^0‘l*E|./0hl ) ) tp< a SUMOA 1 *E(^/DHI 
t KO ■ ■ * • ‘ ..... 


. s F. E 5. iiJLi.A-0j, A 1 1 E 9 i. H 1 1 E 0 8 « ». Q 1 

ItyPE * 3 

GO. TO. ..2X10.0 


PN 9510 
P.M 9520 
PR 9530 
PM. 9500 
PM 9550 
Prt 9560 
PM 9570 
PM 9580 
Pm 9590 
PM 9600 
PM 9610 
PH 9620 
PM 9630 
.PM 9660 
PM 9650 
PM 9660 
Pm 9670 
PM 9660 
pM 9690 
-PM 9700 
PM 9710 
.PM 9720 
PM 9730 
PM 97 ( 10 . 
PM 9750 
. Prt 9760 
pM 9770 
PM .9780 
PM 9790 
.Pm _9600 
PM 9810 
_£iS. J92fi 
PM 9850 

.firt 98.60 

PM 9650 
.JJ 4 _ 9 fl 60 


« * 

, • 

C., 

c.. 

1963. 


.1 L Ow 0 e&l 4 iU.C J 1 QMS.. 


iNTERNAt structure (DRaG iTEh(S)) AT UPSTREAM END OF EACH 


8 3 § 

ENITEM a 0 aTA(-2) 

SOA a DATA.X3J 

CO a OAT*n3) 

EPs ■ 1 Y-QAtA( 16)/1.Q,Q.^ 

IF fENITEM ,UT, l,E-6) ENITEm a 
£k. ■ C0*SOA*£PS*EN-il£rt 


1. 


PM 9870 

-Prt 9880 - 
PM 9890 
..PM .. 9900 
pM 9910 
-Prt -9920 
PM 9930 
_Prt... 9 . 96 Q . 

Pm 9950 
_£M. . 9 . 9.60 
PM 9970 
. Prt -9980 
PM 9996 
-PJAlXlOOO 
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o o oo <->0000 


P.EB1MH ._..P4fiE— 21 


EKO > FEK0(EK>A0»AliE^0,A^ACHt»EMO8O,G) PMiOOlO 

_PM10020 

GO TO 2000 P^IIOOlO 

c, , PMiOOPO 

Cii FIXED, KA4OrtN«LOCAL-L0«8 ITEm AT UPSTREAM END OF EACH OuCT PMlOOSO 

C ,. PMlOOPO 

jPTO NN » 36 PM10070 

_EK_ i_oAx*aii _ pMiooeo 

EKo . FEK0(EK,A0,*l,EH0,AMACHl,EM08a,6) PMlOOPO 

__ I TYPE ■ f PMlOlOO 

C,, PMlOUO 

Cri*-*COMEaN£lilJm*attTlJ3il-Mf.gAlNCEI).-.ta5LS..-ai!Kl4mXlJi - PMIOIEO 

c.,, PHlOIlO 

2000 SEKOf Kil ■ tK o PMlOUO 

SUMCKO a 8UHEK0«EK0 PmIOISO 

PMlOlbO 

C*.,.SECTIoH PERFqRHaNCE information output pMlonO 

pMioiao 

IF tE«0 ,ST, 0,0) GO TO 2001 PwlOlPO 

»«RlT Ef6,9 SflO ) M PM10200 

TLiST ■ .2, Pm10210 

. I.I»^ECLI„.i„LlNEtl*.3. PM10220 . 

«0 TO 2o0 PM10230 

.2001 ,C Akl_.Q.UIPy.T.t.mPi,.!!tN J J,_. PM J 02«0 

lInECT ■ lINECT*3 PM10250 

GO TJ 2(lii _ LL PM10260 

PM10270 

************^jj *jftjk**«*tJLt.****t*AAf>** A ******* A ***A.A* ft jJA*** ********** A.* fl -PM.102 b 0 

* * * • . V. . , . * PM 1 0290 

»,C.A8£.lEjjHJLNAllQN._TAa!l8--- - - - PMlolOO 

,,, PHlOSlO 

* t At t , 1 i — - .... — ....PmIoJEO 

,,, .summary *nO overall performance Cai,CUl*T10N8 PM10330 

,,, PMlOiOO 

,,,,, PMIO3S0 

C.. .iSRijR-CAUSED.SK I P..A aiP HEssaGE . . PMIOJ 6O 

C,, PMlclTO 

.J.0.0 2. ..I F .. ( ( I.k 1 ST-.Ll * --• U . Qa* A B a t s U ^ E K 9 ) , L T . _ I . E • p 1_ ^A N 5, P H 1 0 36 0 

1 UNECT GE, (LlNEMx.ijj IPaCE > IPA 6E*1 Pmio390 

IF ((T IIST .lT , ,1. .or. A68CSUMEK0) AT^..UEi,6y.^Nj?. PMlOROO 

I l^nECT ,CE, (lInEmX-D) mRITE{6,9001) ITITiE»IPa 6E PMlORiO 

I_F I ( T IJ_S T „ ,J.T ^ ,1 OR A 8 S ( SUMS K 0 ) . , LU- L, E«6A._.,ANe., PM 1 0 « 2 0 

1 LINECT ,GE, (lINEmx-D) WRITE(6,900R) PM10«30 

IF. ({TWLST .LU_.l...^OR^..A88(3UMEK0j-,i.T^.,E.9) .aNQ.^ .PMiORRO 

1 LjNECT ,GE, (LiNEMX-n .AnO, tU , 6Q, 1) MRjTEfP.POOS) PM10R50 

IF ...QP,- ABStSUMEKOl. *LT.,-1,E.9) .♦ANB., PMlOttSO. 

I lInECT ,GE, (lInEmX- 1 ) ,anO, lU ,Eu, 2) wRlTE(6,9006) PMloa70 

IF aTLI8T_,LT,..l..,0R, ABS(SUMEKO) ,WT,.l,E*bl .ANO, PM10480 

1 LINECT ,GE, fLINEMX-D) LINECT ■ lO PM10«90 

IF tlLI.8T. .IT, «l, .,0R, a8S(3UME«Q) ,LT, 1.E.6) WR1TE(6,9501) _ PmioSOO 
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PiRFDHr. 


Page 22 


IF (tU8t .LTt •!. .OP, aBSOu'^EKO) .Lt. UE»8) LINECT « U^tCT + S PHI0510 

IF_ (TLl5I.-»LU_-a^--»0fi^-.lMi&U'ieKQl.^l.T,_i^.A.)_.aQ_ia JU515 PkIo520 

C P«10530 

C,, EJjERjgJf. HaT-IQ .. . PnloSflO 

C,, P^IOSSO 

l9..»- U/SUM£Kil . PMiOSfeO 

pmiosto 

C,, BfiESSUSk-filfltaEAjUL ACRt)88 SECTIOA»_Jf!AWU_X£flMPUTEJ) AT THE PM10580 

C,, OOt^NSTPEAH E^O OP EaCH SECTION) P-'tOSPO 

C,*.. . - - - . Ph10600 

Do 2005 I a 1,N PM10610 

IF (I -♦?<£,. 1) .60.10 ^QOl . . PM10620 

SSUWfLd) ■ SPLei) PM10630 

SSU^FOUl 1 .gE-KHiJ. - PMlObOO 

SO TO 2oo« PM10P50 

.2'JOi S§'J.HEUL1. l_8S0ilELcI?lJtSELcn — PMlOfcfeO 

SSOMKodj B SSU'^KOf I.1)+S|K0(I) Pr'IOPTO 

2000 iF (iFl»G ,.|Q, tJ 880 .HkO(t) ■ S8UMkO(.)«8UM£kO PHlOOeO 

2005 OEtP(I) « PATH«r(PT«Q0*S8UMH0(I))/(l,A(G*l,)/2,A8MACH(I)**2)** PM10690 

I (6/ (.6-1,11) PH10700 

IF (LINEcT ,LT, (LINEMX-D) 60 TO 2006 PmIOTIO 

IPACfc a 1PAG£,1 . _ Pf', 10720 

wPlTt(B,900i) niTLE.lPAGE P*^\O7J0 

WInECT * 6 - . . - - P‘<107<I0 

2006 IF (III ,e 6, 1) wPITf!(6,9300) SUMpL PM10750 

IF (lu. , 10,-21 -■^PlTFlP,’5.Ql) 8UHE)U „ PH0760 

LI'JECT a lI^ECT^a P«10770 

C , , , ,., — ... - - _ Pm 1 0780 

c,, paging.checf pgiOg to Output of sumpabv pepfophance pmiotpc 

c,, _ INFQRMaTIQN PmIOOOO 

C,, pMlOdlO 

IF (Ll?sllLT.,uJ,_aIM.HX-«U-60 TO 2007 ... PH10820 

IPaGE 8 IPAGt+l PM108SO 

- “PITE £9,9001) -IimEi-IPAGE Pm108P0 

UNECT a g PMioeso 

2Q0T IF (IU ,e6,-2) 6&-.I9 2008 - PM10860 

C,..,, PMloOTO 

C,, PO'*ER.CAl,CUtATIo.HS FOR SI Units _ __ PMlo880 

C,, PM10890 

PrtRIP « A0*W0*A36SUHiK0aRH080aa272,/«MO8F . PM10900 

P»<ROP a P»iPIP,100,/ETAFAN P«10910 

AVGPNP s PwPoP/ENFaN _ PH10P20 

-PITe ( 6,9302) SUMEKO,EB,PWPIP,PWPOP,AVGPaR,ETAFAN,ENFAM PmIOPSO 

LInECT a I.INECT+5 Ph10990 

GO TO 2009 PM10950 

C...,, . - PM10960 

c,, POiKiER Calculations pOp u,g. Customary untts pmio97o 

C-, - . Pm1o980 

2008 PwpIP a AO*VO«*3ASUMEKO*nHOSO**2/ilOO,/iiHOSF PM10990 

P/iPOP « PNRIPaIOO./CTAFAN . — PMllOOO 
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f’A&t- 2l 


AVqPi^R ■ PWROP/ENFAN 

RllE_tfex9ig3l_AjJjieK(>iERi-gj!>i 

LINECT • UINECT+5 

«O0‘> 


PMUOlO 
f»«ll020 
P«4U030 
.PMllOiO 

,,,,,,, PMilOSO 

....vtlpC lTy ADJJia-TMEj^l.-iAkCiik4.t Iaw5_f -INFUT. PMi 1060- 

,,,,PO>^ER VtLUE — DETERMINES APPROXIMATE MAXIMUM TEST SECTION VEI.OCTV PWllOTO 

• t.t^aR SPECIFIED POWER \ E VEf PMltOSO 

PMllOOO 
_.JlllU0-0- 


FAlCTN (OH0 ^O4-EMQ.^t. M(ii£.i_ 
Vft 


CALL 
WOC 

VO ■ - V 0 * £ P w R M£ -H z PJ»IMP > * *.£ -L, ^ 

£MF a fcMFavO/VOC 

EmO a.6;.«o#y.0./yo.C 

RHOSF a RH0T/(1,«(G-1.)/2,*EMF**2)**(1«/(6«1,)) 


R MoS 0 a RHqIZ 1 1 * t £.S?i »J./2 ».a£m 0 at21*-«( 1 » ) ) 

RnoC a RhOSO*VO*AO/£MUT 

Call FHlCTN(,DM0,AO»lMa,8LMOAE) - 

IP (lU .ECl, 1 ) PwROP a Ao*Vo**3*SUMEKo«SLMOAEaRHOSOa«2/8LMOAC*SO,/ 

I RHOSP/ETAF-Ati - RM-LL^O... 

IP (lU ,ED, 2 ) PwROP a AOaVO«*3*SUMEKOtSLMOAEaRHOSO*a2/8LMDAC/l 1 « PMitSlO 


PMUUO 

PR1U20 

PMU130 

PMlllSO 

PMtUSO 

PMUUO- 

PMIUTO 

P-MliUO 

PM1U90 


2010 

2011 


I /RH05P/IT*P*.»: 

IF (aB3((PwRmcm. 
ip (ly .aQ, ii- yos-i-voaURRift 
IF (I'J ,F&, 2) VOK a V0*,S92 a 8 
UO a RH380*V0**2/2, 

IF (uINECT ,LT, (LINEHX.iO)) 

IPA(?E. a IPX6E*1 

WRlTEf6,900h lTlTLe,I***6E 
iF (iU_,FQ,_ 1). WKiTEL6#930«1 
IF (lU .E9, 2) wRITE(6»9305) 

IF (Ip.fiUT _Q1 6 o., TJ3 20l« 


PwROPj/PwRmcHj ,6T, l,E«6) 60 TO 200^ 


RHU220. 

PM11230 
JH 11490 


60 TO 2010 


PMlUSO 

_PMUlfctt_ 

PmU270 


PWRoPjy.QjViaKAlilOjfiO- 

PwROP»Vo»VOK,EMo,00 


\ ! I Ic r«i.Oi,l_SUMMAR.T_j:H ARiJCJERlStlCS Pa6E. 0U.TPUI... 
LliiECL t 


l.N 


2012 

20l3 

20l9 


00 2013 I a 

IF ClInECT. .i.,T ,.-.lI.nEmX 1_.6Q TO. 20l2_ 

IPACE a IPAGEaI 

- R J T 6 ( 6 j 9 0 fill .. II I.T LF.j lP.A5E- 

w«IT£(0,9aoU 

iF tiU. .EQ. 1 )..wR-IIE(6j,9402) 

IP (lU ,EQ, 2) “RITE(6»9«03) 

L1.^ECT -1 .15 

write ( 6, 9400, I,SSUMELfI,,8MACH(I,,S8UMK0(I,,0ELP(l) 
LInECT ».-LL^LCTal --- 

CONTINUE 

IF (IPLOT ,EQ^.0v6O to 20l5 . 


PMU290 

P M 1 1 3 60- 

PM113I0 

P.M1142Q 

PM11330 

PM11340 . 

PMU350 

PM1136H- 

PMIUto 

PM11380 

PMU590 

e.M.U9oo 

PMU910 

!IMtl920 

PHU430 
PMXUOO 


C,.,.CIbCuIT SUMHAbT characteristics -PLOT £s). 


PMU950 

.,fPll960- 

PmU«70 

PM11480 

PMU490 

-PMllSOO 
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o o o r> o non o o o 


PERfaH^; .. - . PaUE 


C»*» 

Cali. pLuTlTtr.iD£LPiSSUM£.LjS3U'^KU*lu«lPLnT,niTLE4TR£TRN,P.L0T0N) PMU520 

P»U!>30 

AN'jpTA.TtD .TaSULA.TIO!'( M.lNPyT DAlA,.CARc5. FUK. CURfiE(LT CaSE, P.mUS«o 

• • • P^^ 1 1 550 

20i5 If CTLISTI »GT,»..l*£-P-40a*-TLl8T ♦LT, •♦ft J, XiiU.- tUIACJsa) . PMUSftO 

P*-Ui>T 0 

,,,tfcN0-0F-C iSEs OR SILT UR CHECK. . PMU5B0 

PMUS90 

tft.9007) ITITLE pMlibOO 

IF CTfifeTW''! ,GT, l«E-ft) Go TD 100 Pi'iUftlO 

IF fT»£T8N ,lT._ ,,5v ..GQ m.ici2 pMllbZO 

STOP 1*^66 PhllftSO 

- - PWUfcAO 


C 

c 

f: 

c 

c 


c 


PMllbftO 

..I^PUT REaO Fir.'aTS PMilft70 

... - -P.MliftBO 

7000 FURt'AT Ml,19Aa,A5) PMjiftftO 

7u0l FGR.*aT a2,ali,4A,ftF5 2 j __ _ PMH700 

7002 FuwmaT (2i2,2M,2F2.0,l4f5,25 ” PmUTIO 

.ti... , PMIJ720 

, .OUTPUT FoRmaTS PMH730 

... P '111790 

.......... PmU750 

«...P*^RF'OrmaNCE I.NiFOaOAYiUN La8|LL1'^G A^U OutpUL.FOr^aT.S. . P^'ll7ft0 

... PMIJ770 

ROOO fORwAT flWt//20)(^Al,19A9,A3,l3x,9HpAGt,13//' P«117ft0 

1 2ftA,944,29M »iIn0.tuNNEL PEftFOftMANCE/) Pi^ll790 

9ool format I'1H1//S* Al,j9A<(^A5,ftX,16M,,,C0^'TlNUEp ,,,ftf,4HPAGE#n//) PMilBQO 

R002 format f2«H aTmOSPoERIC PRESSURE « ,Fft,i,15H ATMOSPHERES ■ ,F9,l, PHUSIO 

A ftM v/SiJ M / _ _ _ PM11820 

M i7M TEST Section coroitiqns .« / pmusso 

C. 2li" TC7 TaL.PRE 8SUR£..« ...Fft, 3. ISN. atmospheres «.,..F9,1. PRliePO 

n SH N/sa m;/ pmusso 

E 2«H , total temperature a ,Fft 2 9H DEG C f ,F7.2.7H OEG •<,/ PHllftftO 

F ISH velocity b ,FT,2,9H m/s£C t ,F7,2,28H KNOTS. DYNAMIC PRE PMU870 

GSSURE *. |Fij, jigH m/SQ m,/) . . _ , PMllSSO 

9003 Format (2«m aT*,0SPHErIC pressure » iF6,3,15h atmospheres « #FH.2| PMne90 

4 lOM Le/Sa FT^/ . _ _ PM1J900 

^ ?7h tEsT section COnoItIONj .• / Pm119Jo 

c 2lH. total pressure . 1 ,Fg,3^lSH .ATMOSPHERES 8.,fS,2, Pm11920 

0 lOh LR/Sa FT,/ Pm11930 

F 2«« total TtMPER.ATyRI.« ,Fs,2 9H qEG F...»_ f7.2.7H..DEe «,/ PMU940 

F i5H VELOCITY * ,F7,2,10H FT/SEC • ,F7,2,2»h KNOTS, DYNAMIC P PMU950 

.GhESSURE *.,.F7.2ilO.H lh/SG ft./) PMll9ftO 

9fM)9 FOfjMAT (ISOM NO, jECTION TYPE HI “1, 01 AREAl PMU970 

A ai/a?i aR.CR 2 TmETa . Vi HaCMj CENffTH nP/QL OP/QO PMll9eO 
M / SOK, 33H H2 *2,02 AREA2 a2/AO ,17y, PmU990 

C 22H v2 maCHS .... i . PM12000 
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EfiJEMi! ... EASE. 25 

900W0RMAT(2ax,26H METgRS MET|R8 ~Tq ^ 17X, 32H0E6REES M/S 6C RM12010 

t meters /120H f» .RH.12Q20 

9 ^mmm^ ^mwm^ ^mmmf ^ummmif ^***4 ^'■■*** RM12030 

f* ) ,...PM12Q«0 

900<> Format (50X» 26H feet feet SQ ft , l5X,3lM0EttREES FT/ PH120S0 
AS .EC FE£I /12 QH » .» + - m ♦ * P M 1 2 0 6 0 

A******4 t pMl2070 

e»» »» *•••* ^ PMiioeo 

9007 format (/9H ** i,Ai,l9Aa,A3,39Hi CAsE" C0HP1,EtE6 OR TERminaTEO, ** PH12090 

.A J. — PMiSlOO 

PMi2tl0 

C^.^. summa ry PERFQRmaMCE |,A 8 E^l I>jG AnD OUTPUT FMmaTS PM12120 

C,,, PM12130 

_?300_ f ormat f9«X.8H»» — «»»/6aX f 2bM TqT A| _ CEnT ER|_^1m E LENGTH • , F8^2^... PhUUO 
A TH METERS) PM12150 

9301 for mat / 68X.»8H TO TAL C fL MTERLlNE LEh lSTH. .« -P H 1 2 i 6 Q 

A 6 H FEEt ) PH12170 

93c2 format (21HoPER F0RMaNCE SUMHa RT „/ „. pMi2i60 

A IX, SOM TOTAL Pressure loss (DP/qO) ■,FS,5,6X* l9hE^ER0y RATIO ■, PM12190 

— ,e. -f7.3/ .PH12200 

C «X,14mT0TAL power ••/ PM12210 

J2 .Uia INP UT TO PLOW cy J PU I_.R^aiUP,£a AY18A.6£_ P t R. F A N p H 1 2 2 2 Q 

E Fam efficiency TOTAi nUmBEr qF FanS / Pm12230 

_ F . _F 1 5 ,0 , 6 H watts, F l3,0,6H W ATT8 ,FU , Q , 6 H WAT TS ^ F li.* 1+9 M._...PiR CE.N T , P M 1 2 2 « 0 

G Fi6,0) PM122S0 

93o3 for mat f23MQPERFQRWANCE S UMW^Ry „ / _ . pM)2260 

A 3X.30M total pressure LOSS (OP/ 8o') •,F8,5,6X, j'«MENERGY RATIO ■, PMi2270 

B FT, 3/ PM12260 

C RX,1«Mt0tal POwEr ••/ Pm)2290 

_ _Q__iioH INPUT TO PLOW OUTPUT R EUUlR En AVERAGE. PE.R FaN PM12300 

E Faa* efficiency total NijHBER OF FaNS / Pm12310 

F FlS,0,4M MP,iriS,Q,9H H P pi6 ,0,qM HP P 13,2,9M . PERCENT , F U . 0 ) PM12320 

93oR format (///S 5X,5 oMMA XI mum VELOCITY FOR a SPECIFIED POWER cONSUMPTI PMJ2330 
40N //aSM TmE Maximum test section flow .achievable wITH.FIQ.Q, _ pM 12390 
B SEN -ATTS nP Po'^ER AVAILABLE 1$ APpRoXi mATELV AS FOLLOWS ••/ Pm12350 

C 15X 13H vELOrlTY .. FB.a . B H M/8E C ■, F6,2.bH KN OTS/ . PM12360 

“ 0 15X.16M macm Number -• ,F5,2/ Pm!2370 

E 1SX,21 H DYNAMIC PRESSURE ■» ,F9 ,2 ,.7 M _ N/S.O _M ) - . P m 1 2 3 0 0 

93o5 format (/// 35 X,SoMMAXIMUM VELOCITY FOR a SPECIFIED POwER CONSUMPTI PMi2390 

AON the maximum test SECTION FLOW aCHIEvaBLE „> 1 T M , F 0 , 0 P.M 1 2900 

B S2 M HQqsEPOwEr available Is APPROXIMATELY as FOLLOWS «/ Pmi2<*10 

C iSX.)%H velocity ,Ffl.3.9H F T/SEC « ,F8.g ,fcM KNOTS/ Pm12«20 

0 ISx.lEH MACH NUMBER •• ,F5,2/ PM12930 

E 15X.21H DYNAMIC PRESSURE ■■ ,F7,2, 9M LB/8 Q FT) Pm12«90 

C, PM12950 

C.... CIRCUIT summary information P aGE FURMaTS PM12960. 

C,,. PM12970 

990 0 FO oMAT r25 v.l5,Fl7,2. F19,3,FlS^5.Fl7.i) ___ . PM12960 

9901 format {38X,aaH WINO-TUNNEL CIRCUIT characteristics summary/ PM12R90 

A ji0 X,RQH takEm at d ownstream EnD oF EaCM_ sECTlO.M/y - PM12500 
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.....P1RF0-Rm Page 26 


p 2SX,71 h section cumulative MaCH CUHUUaTIVE PM12510 

. C i^ALk_.P^£aSJJfi£y, PMi2i20 

n 25*,71MAS8 IGned CIRCUIT NUMBER PRESSURE PMt2530 

.. „ -£_ D I.FF1R£N11aL^ P.M1256 0 

P 25X,7lhSEOUENCE LENGTH LOSS PM12550 

S.-£JLl!LQS&HLRIC_y _Pji42560 

H 6RX,27m(!5p/Q01 • internal)) PmUSTO 

, 9Ao2 forma t f2SX,7Q H Hg TERS PiUlSflO 

A N/S8 M / PM12590 

9 2SX t J AWff<»««IPJ».AA ^mmmmA ATy****^^* Bill26i).0 

C PMI26I0 

. 9R0?_ Form at f25i<^ 7QH FEET Pltl2620. 

A L8/SfJ FT / PM126S0 

_ §..... 2 5 JL». 7 1 f ♦»■■■■■■ a* PLM-1 E6il G. 

C PM12660 

r,,.,, , . ,, , Eilti6AIU 

C..t»ERRfiR*t^J*®MpSTTC MESSAGES *N0 Formats PM12670 

C ... : : — EHllMa. 

9500 format HHi) PMI2690 

9501 format f/li8H ■>*. OU e TO ERRORtS, IN INPUT CAWPfS), VaLIQ SU MMARY _PF1jL27aQ- 
ainFCRmaTjOn iS not AVAtlaBlE^^ ^FER To THE TABULATION OF INPUT / PH12710 

„ _ fl__ i 2M DA TA ON THE F 0 L LO R ING P A GE8'. CORRECT THE E R R QR l.S J_..ANC_ ?. M 1.2 120... 

c ris'ubmit this Case; subsequent Cases rill not be affected,) pmi27So 

952_2 format title ( i . A 1 . ^ 9A q. A3. R7HI ) 18 INCORRECT OR IMPRO PER / PM127R0 

i i20H AS IT EXISTS, THE FIRST CaRO COLUMN MUST CONTAIN AN ASTEr PM12750 

risk if\ T O Hi I DENTI FIED Aj_A VaLID TITLE CaRD J._ PM 1 2760 

C 28h tHjS case WILL BE SKIPPED,/) PmU770 

RSoi . forma t Master c ontr ol data 7fs.2.6ih i ) is inc . (yu2,7 6 Q.. 

aorre‘ct or improper as tT Exists, the f*rst twq card/ pmi279o 

s I££H columns must contain a NEGaTIVE NUHbER to »9) TO BE ID PMi2SOO 

CEnTIFIEO as a valid master Card', this case will be skipped,/) PH12810 

95o« fo rma t f//)17H MORE T haN ONE MaSTER CONTROL CaRR EXISTS £QP._IH18_ PM.l281.fl 

A Case or input CaROS a’re out 6f order, check deck set-up, / phusso 

R 1 2f)H TMg last master Car d EN C.QIMJEMD KI L L BE -laS UM EB-AS. JiiE-Cl) -P Ml 2 8 4fl. 

cnrect Master caRo for the section caRos which follow, /) pmusso 

-- 9 5.o5 FDR Mat f/atH m aximum LIH Lt_QIO HE number OF SECTIONS (. 1S.69H) HA . f W12fl 6 0 
AS been reached. EITHER A CASE TERMINATION CaRD HaS BEEN OMITTED/ PMUBTO 

C 12f>H (ALONG WITH TITLE a N 0 HU T E R_£aR DS TO BEGI N -A_N£.M_ CaSEJ- QR-PJtjiB B-Q 

0 ThjS Case j8 too long for The programmed AllOHaBLE number / PM12890 

g 59H OF sect io ns . I.H£„CLa8,L HaS BEE N TERMINATED aT THIS f,OJiLU/4-^M12900- 

95f)6 format f/ll8H master CONTROL CARD MAS BEEN ENCOUNTERED BEFORE CaS PM129J0 

aE T ERwi.si AT iOn and TtT| E CAfiOS, CHE Ck OECk SET-UP. I P»4l291fl.., 

R 120H ERROR-mESSaSE TITlE wIlL BE GENERATED AND SUMMARY OUTPUT, PM129J0 

CNO-PLQT INPUT Data Tabu lati o n an 0 . next- caSe rpturn / Piil 2 9 R o 

0 <hH termination Parameters will be assumed,//) pmurso 

9S n7 format f/ 60 H »« NOTE — TEST SiC-TlPN- BLO CK AGE FRO M SE CT IM iARft_l_E!U2?6il 

a»,PUT (,FS,J,69H pERCEnT) OqES NOT EQUAl THAT OF THE MASTER CARO/ PM12970 

, 8 8H I NPUT f, F5 ,l,lQ9 H PERCENT), CHECK DATA DECK. SEC TION CaRO V. PMA29BQ 

CalUE will be assumed as CORRECT anD EXECUTION WILL CONTINUE, ) PM12990 
95f>6 for mat- f/>i 7H altho ugh _.V£LQC.m.P PT IHI2ATI QN - H AB R E.flU28TE0_fl-Y ...P-MlSOaO 


77 


££fi£jasjl 


£a6E SH 


ATErMINaTION code, the INpUT pOwER VALUE IS ILLEGAL (LESS THAN OR/ RMllOlO 

8 180H EOUAk-10_l£Ml, _ IttEBEEQRf A_NQ_l£LQillT Y«QPTlHI7l NB 18 PQ_fM.UL020 

C6Sl9|,e. RECHECK InPUT ValUE qn hASTER DATA CARO, ) PHlSOSO 

15.ai -£.6 P-H At - [ /I tSU AA ERROR some INCO RRECT COWaINaTTON OF (HPUTS OR ■ PmAOS 0 

aunantiCipated Situation has caused an invalid (non.positive) / phisoso 
B J4H total l oss LEVEL^_J»XCH£CK__S£cmN-,-ll,-llJtLJma.J)llA^--PJUlDAD 

ENo PMIJOTO 


^ . 
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DAT4£K 


subroutine oaTackinlist) ok 

C * • t * * * r>K 

C******************** A************ ******* A****************************** OK 

C THIS SOUTI NE A subroutine OF THE _Ha IN. PRQQRaW PERfQRM, CHECKS KOR OK 

c errors in inRuts or Master control card and section cards, and, ir ok 

t . RESUE8TED, I T rONTRQLS THE aS SEHSLY AND PRINTINC OR THE ANNOTATED OK 

c Tabulation or the input inroRmation, ok 

£ A*,* * * ** *** * ******** * ************************************************ *** OK 
C**«* ******* A***************************************^******************* OK 

COHHON/ 8LQCK4/i8EO,ISmaPi,I8haP2,N OK 

COHHON/BLOCK8/OATA(U),ITXTLEf21),IPA6E,lPLOT,IPRINT,I8EC»ITUNNL, OK 

L- - IUALIMCT,LlNEHKtR R .R-USERR Q M>.I tT i TL l tT l, TR E TRN OK 

OIhCnSIOn EmDATA(13),EmWRIT(26)iCn0ATA(30,20),ENWRIT(R0), OK 

L. ETW RIT(l2i^MCHECKfl3),MOATAjS>,MRORKT(3Q?,M*RITEM2l. qK 

2 NCHECK(Jo,2o)#NOAtA(3e»«)»NRORMT(R2),l^HRitE(8) OK 


^t^toB^K'TlHE RORHaTTINC aRRaYS 


M. 

OK 


10 

-IS 

30 

-00 

SO 

-fcfl- 

70 

-BO- 

RO 

loo 

110 

12iL 

130 

IRQ 

ISO 

iKfl- 

170 


c ♦ I 


DATA ML£PT,NLERT,IaRL02,1ARLOR,1IRI.OS,1RLOO,1RRLDO,IRRLOI,1RRL02, 

_L IRRLD 3,lMLDR ,i COHHA, i 8PACC, l 8PACl'lRIGHT/ 

2 4Mc1Ry,«H( ,3HA2,73HAa,,2HlS,4MR5,0,RHR6,0,ttHP6,l,RHR6,2, 

3. .4 HR S.3,4HP6.4,1H,,4H,1K,,3h1x,,2M/)/ 


C* !**IN^UrT A8ULATlON annotation HESSaCE aRRaY (INTEGER VALUES^ 


C • f • 

- 2 4 HHPTY^ !___ 

c!!|!!NS urTA8ULATlQN annotation ME88aQE AR R AY {RLQaTING POINT VaLUESi 

data RBLN K2,RBLNK4,RH861,RH 862,R H863, RH8S 4.RH8 gS, RH8GB,RH8S7, 

1 R*^SS8/2H ,4H ,2H E.4HXTRA.2H E.RmRRORsH 0,4hPTiN,2H E# 

4HHRTY/ 


c| !!!H S hination Card parameter translation array 

C i i 

*_ * da T a TH 8Sl,Tw8G2,TH8G3.TH8S4,TM8B5,TN8G S A,TMSShB,TH8G7,TH8G8, 

1 TM8G9,TH8Bie,TH86n*TM8Sl2,7M8613/ 

2 ■ 4HYE 8 ,4H NO .gHN0NE.4HPRE8.4M8. L.4HQ8S .aH08S..4H HAL,4HL PR, 

3 4M£8S.,4H(CHo,4M8En)»«H(Ror,4HcEO)/ 

— t 

C... .transfer to applicable section or subroutine •• 

part 1 (STATEMENT 1000) ROR MaSTER CONTROL CARO CHECK.OUT 

C., Part 2 (STaTEHENT 2000) ROR section card check-out 

_Jtx4i PART 1 fSTATEHENT 40001 RQR O U TPUT TABULATION OR INPUT CaRDS 

1!_* g o to rlOOO. 2000, 4000), NLlST 

C 

C*** .*-*A *. ********.*************************** * **.**********■ ***************** 


-- DK U-0- 

DK IRO 

_QK 20Q_ 

OK 210 
OK _.22Q 
DK 230 

DK 2S0 

_JU< ifcO._ 

DK 270 

-DK 2Jfl_, 

DK 290 

OK 310 

—QK 320 

DK 330 

DK IM- 

DK 350 

—UK 360 

OK 370 

OK 390 
DK 400 
OK 410 

__QM 

OK 430 

—OK 

OK ttSO 

_DJS ftfcO- 

DK 470 

—aK— aaq- 

DK 490 
*._0K ioo_ 
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XUUCK 


a 


c|.D £hNlfioroF*MltfER CONTROL INPUT WgOUIREH£.NTl. — 

C., CODE lOi indicates INPUT NHICH 18 NOT REQUIRED POP ANY 


DK 910 
JUS ilO- 


C.. 


PURPOSE fN n WHTCH MAY BY QMlT TfB PROM IMPUT CAHB 


OK 930 
DK IRQ. 


L 


CODE •)' indicates mandatory INPUT OK 990 

CODE iPt indicates optional INPUT WITH DEFAULT PROVISION OK — 9*0- 

CODE lit INDICATES N0N*RE0 UiRE0, CONVENIENCE INPUT OR DK 970 

optional input NMICH May be correct as EERO OJS—iAO- 


0 • • • 


c ,,, .integer data 


OK 9r0 

JUS — *00- 


1000 MDATA( l) m ISEQ^^ 


DK 010 
-OJS — *i0- 
DK *30 
JUS **0- 


MDATAfS) 

MDATAfSM 


lu 

ISHAPl 


DK *90 
-OK — **0. 


hOATA( 9) • iSHApj 
C*ICt^t^*nN6 POINT Data 


OK *70 
JUS *10- 


DO 1001 I ■ *.13 
EMDATAfI) w DATAfI,31 


DK *90 
JUS 700- 


1001 CONTINUE 


OK 710 
JUS 710- 


l-t A AAJ 


C....1NPUT requirement DEFINITIONS 
.0 . . ... 


00 too* I ■ 1,13 

MC HECKfll ■ < 


1001 con t^nue^^ 


mCMECK( 5) ■ 3 

IF fl SHiPi ,EQ, is MCHECK< *> ■ Q 


MCHECK( 0) « 3 
JlCHiCUXioJ-S-JL 


OK 730 
JUL-lOO-. 
OK 750 

. DK J!*0- 

DK 770 

-OK 7SO - 

DK 790 
OR— -SOO- 
OK 810 
JUS- -820 


mCHECk(II) • 2 
-N£H|-£Kj12I_9-X 


OK 630 
._OK — 8*0. 


MCMECK(13) • 2 
cttTT^Nfi&lR input error.check 


C* 


00 1003 t « 1.9 

IF ( McHECKtl) ,NE. 0 ,aND, H DA TA.I) 
1 TIIST • ..9 

IF (Mr.HECKfI> ,E0 , 1 ,aND, M QATA(I) 


1003 CONTINUE 

C.. 

c • » 
c.. 



UNItS»oF«MEASURE Error OEtECTIoN 


EMERR ■ 0,0 

IF- ( m«2j-l- ,Efl»-fl v-fifl. T Q . t00«_ 



850 

dk. 

-_8fcO- 

OK 

S70 


-8*0 - 

DK 

690 

rsK 

900- 

OK 

910 

-93 OK. 

920 

OK 

930 

(US . 

. .990 

OK 

950 

OK 

9*0 

OK 

9T0 

DK 

980 

DK 

990 

OK 

1000 
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JULA£JL 




MRlTE(6,a006) ITITlE»II>*GE 
MLTJUA«.gQ02) lU 


LINECT 
JU l-l. 


■ 6 


EMERR ■ S, 

1000 IF 

IF (EMERR ,er, 2,Sj GO TO 1005 

HRlTEf6,e oo6) ITITLE, I RiGE 

LINECT ■ 4 


lOOS i»RJTE(6|80QI} I8HAP1 


TLlST ■ ml, 

100 8 IF f Cfl8HiP2»n*f I8HAF2»2)»rISHAF2«I 

c !!**^ l62Hng point input err or.cmeck 

c 

P0_ 1007 I ■ 6, 13 


IF fMCHBCK(I) ,NE, 0 ,*N0, A88 (EmDaTA( I) ) ,|.T. 1.E.6 ,ANO, 

1 T LI8T ,6T. » , 5i TU &I j_m5 

IF (MCHECK(l) ,Efi, 1 ,AN0, AB8(EM0 aTA(I)) ,CT, l,E»6) TUI8T ■ •». 
1007 CON TINUE 


IF (TLI8T ,0T, *2,5 ,0R, TLI8T ,l.T, .3, 5) RETUR^ 
wHlT Ef6y8000) 


IF fCMERR ,i.T. 2,5) MRITEr6»8006) ITITlE»IPAGE 
IF fE^ERR ,GT. 2.5> LINECT « LINECT*4 
ITTenerr Itt* Ilsi LiNECf 4 

AaM*L 


OK 1010 

OK 1020 

OK 1030 
— OA-4060 
OK 1050 
1060 
OK 1070 
QK 1080 
OK 1090 
-OK_1100 
OK 1110 
^_0K--U20 
OK 1130 
J)K 1140 
OK 1150 
._.HK_ll60. 
OK 1170 
_DK 1160 
OK 1190 

QK .1200 

OK 1210 
...llK.a 220 
OK 1230 
12.40. 
OK 1250 
__DK_1260. 

OK 1270 
....J3K.1260 


OK 1290 

>*»»*»**«»*<>»» r>K 13 0 0 . 


f, .DEFINI TION OF SgrTlON INPUT ReOUIREHENTI 

C«* CODE lOi indicates INPUT NHICH 18 NOT REQUIRED FOR ANY 

Cmm PURPOSE aN D mHI CB. .MA i ,.6E_imiTTEP-F.R0M IfjgljT. CaBD 

C-, CODE «i» indicates HANOaTORV INPUT 

C... tQM_.i2i_lNDl£.Ali.8. ■ " 

c«« CODE (S' Indicates nOn«reguireo, convenience input or 

. optional input wh ich Hay be corr ect aS_.Z1RQ 

2M1L.8ER RQ R R.jLJ) 

C, 


OK 1310 
.DJ<_1120 


Data 


DK 1330 

DK-1340 

OK 1350 

OB .1160 . 

OK 1370 
_ JlK_13ao 
OK 1390 
_HK -Uoo 
OK 1410 
H(L_1-420 




n 0ATA(N,2) • ISEC 
Jt OAT -A iN Al) r llHAPl. 


n0ATA(n»4) • I8HAP2 

c • • I»float1ng point Data 

D >»» 


DK 1410 
_QK_U40 
OK 1450 
_0JL1460_ 
OK 1470 
fiK_14eO 
OK 1490 
.OK_1100- 
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20Dt CONTINUE 

C,,,,CENER*L IN^UT REQUIREMENT OEFINITIONE 

C I « « 

DO 2002 I ■ 1,20 

NCMECK(N,I> ■ 0 

2002 CONTINUE 

NCMECK^N, n ■ a 

NCHECkINi 2) a 1 

NCNECKfN, 11 ■ 1 

NCHCCK(N, 4) ■ 1 

IR fISMARl ,NE, i) NCHECK(N. T) ■ i 

NCHECK(N, 8) • 1 

NCHiCKfN. 41 ■ t 

IF (ISHAR2 ,nE, n nCHECk(N> 10) ■ 1 
NCHECK(N,in 9 i ^ 

c. 

C, ,,.SECTION,TVPE BRANCHIN G 

C... 

IF flSEe ,EQ. n GO TO S OOO 

IF (IsEC .Eq. 2) 00 TO 2020 

IF flSEC .EQ. l\ 60 TO 20S0 

IF (ISEC ,CD, 4) 80 TO 2040 

IF flBEC .EQ, 60 TO loop 

IF (ISEC ,eu, 4) 60 TO 2060 

LF ( I 8 EC »E8, IQl .Sfl-IU .iMO 

IF (IsEC ,Ea, 20) 60 TO 5000 

IF (ISEC ,EQ . 10) 60 TO 2500 

IF (ISEC .EQ, 52) 60 TO 2500 

I F (ISEC ,EQ, 35) 60 TO S550 

IF (ISEC ,E0, 5«) 60 TO 2540 

IF MSEr .EQ. 40< 6 0 TO 2400 

IF (IsEC ,Eo, 45) 60 TO 2«5o 

IF ( ISEC ,EQ. 46) 60 TO 2 460 

IF (ISEC .EQ, 51) 60 TO 2510 

IF (ISEC ,E0. 52) 60 TO Z52C 

IF (ISEC .EQ, 55) 80 To 2530 

IF (ISEP .EQ. 54) 60 TO 254 0 

IF dsEC ,EQ, 56) 60 tO 2560 

IF fIS EC . EQ, 57) 6Q_TQ ^5l0 

NCMECKfN, 5) ■ 1 

IF" |lsfc~TlQ! 62) 60 TO 5000 

1 F -MSE C .E6. 70 1 GO TO 270 0 

IF (IsEC ,EQ, 71) 60 tO 2700 

IF (I8EC_*EQ^ 72) .fiJ3-JtO.-27-OQ 

IF (ISEC ,EQ. 73) 60 TO 2?50 
- IF (IsEC. *EO,-74i-_fiO_Xa -27-40 


OK 1550 

OK tl40- 

OK 1550 

OK 1 56 6- 

OK 1570 

OK 1560 

OK 1590 

PK UftO- 

OK 1610 

OK 1620- 

DK 1650 
ftK 1640 

OK 1650 

OK 1660- 

DK 1670 

QK_X6.aa 

OK 1690 

DK--i700 

OK 1710 

0HLJ1720 

OK 1750 

J5K..174Q.. 

OK 1750 

0-K. 17AQ 

OK 1770 

OA-ITBO 

OK 1790 

- - Q.K .1900 

l)K 1810 

- OK. 1920 

OK 1650 

OK 1940. 

OK 1850 

OK 1860 

OK 1870 

OK..1980_ 

DK 1890 

..OK 1900. 

OK 1910 

DK 1920 

OK 1950 
OK 1940 
OK 1950 
.OK 1.960 
OK 1970 
OK 1980 
OK 1990 
OK 2000 
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-JJlIJltJl 


PAfiE._.-5 -~ 


(I8EC .E8. 75) 60 TO a7«0 
6 0 I fl- a8tlL 


IF (I8EC ,E0, 05) 60 TO 5000 
_1E_ i laEr ,gg, bo to UkO 


IF ^lilc *Iq* B?i CQ TO f??n 


nK 2010 
_0(L-I02fl. 
DK 2030 
-IU1-20JO 


IF (ISEC .EO. 92) 60 TO 2920 
IF-tl&IC_4E0^-J^9) 60 TP 29j Lfl_ 


IF (Z8EC ,C0, 96) 60 TO 2960 


OK 2050 
-OK_206i»— 
OK 2070 
-OK.-iOBft_. 
OK 2090 



fiatc ^7) .BQ TQ 

OK aioo 
DK 2110 
DK 212 a 

-1 

INViLlD SECTION TVPE.CODE MESSAGE 

C • i 

IF rllNEcT tLT, (L1NEMX.2)) 60 TO 2003 

DK 2130 
nK 2)40 


WRITE 6 ^»ll^*tTTLE TP t 

DK 2150 

QlL-1160 

OK 2170 
UK 2180 . 

LINECt ■ 5 

2003 KRITEfO. 80041 M 


LINECt ■ L1 N|CTt3 
SERROR ■ . 2 , 

OK 2190 

OJt_1206— - 

OK 2210 

. pM 222 A 

C t 

GO TO 3003 

C t • » • 

UnO^e.ouct 8 reciai.ized input requIrEmEnts definitions 

DK 2230 
pK 224Q 


DK 2250 
.flK_226a- - 


TEST sections 

c.. 


OK 2270 
pK 22 j(L 

c.. 

£.«. . 

con 8 T*nt»ar£a test section with MODEl 

OK 2290 

DK-.2300 

2020 

NCH£CK(N,12) ■ t 
_NCH|£!<^N,1<U_ 9 3 

DK 2310 
D,K,2320 


NCHECK{N,i 7 ) b \ 
NCHEcLK(N,20) B 3 

DK 2330 

. .. . DK 2340 


60 TO 3000 

DK 2350 

OK 23<ifl 

C 1 • 
L.. 

0 IFEU 8 IN 6 TEST SicTlON •• EMPTY 

DK 2370 

nK 238Q 

2030 

nCMECK(N,i 6 ) b 2 
GO TO JOOO 

DK 2390 

nK 2400 

C t f 
C.t 

OIFPUSIN 6 TEST SEcTlON WITH MOnEL 

OK 2410 
OK 2420 

c » » 

2040 

NCHICK(N,12) B \ 

OK 2430 

OK 2440 


NCHECK(N,14) b 3 
_kCMECitlN.16) B 2 

DK 2450 

OH 24<^0 


NCHECK{N,17) a 1 
NCHECK{N,20> b 3 

DK 2470 
nK 2480 

C-»^f — 

GO TO 3000 

OK 2490 
_08_2400_ 
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D4TACK 


JUfie 6. 


C • i 
c.. 

open-throat Test section with mooeu 

OK 2510 
riK 2^2^ 

2060 NCMECK(N,i2) m 1 

OK 2530 

NCHECKfNflttl • 1 

pa; 

NCMECK(N,17) « 1 

OK 2550 

^CHECKfN.^01 ■ % 

OK 2<^6D 

GO 

TO 3000 

OK 25T0 

Cff t^l 


_ OK 2SS4- 

c.. 

CORNERS and TURNS 

OK 2590 

C.t 


D“ 

Cat 


OK 2610 


.CQNSTaNT-ARE* turn NlTM VaNES 

UK. -2 620 

C t * 


OK 2630 

2300 ^CH£CK(N,^ij i 1 

.. M_2660. 

NCHECK^NliS) ■ 1 

OK 2650 

NCHKCKrN. 16) ■ 2 

. .nK .2660. 

NCHECKfN.lB) ■ 2 

OK 2670 


TO 3000 

_ UK -26fl0 



c.. 


OK 2690 

C^a 

constant. AREA TURN WITHOUT VaNES 

. „ .OK .2 TOO.. 

c., 


DK 2710 

2330 nChECK(N,15) *1 

. OK 272j1_ 

NCHECKfN,l6) ■ 2 

DK 2730 

GO 

TO 3000 

DK 2740 

C.aatt 


OK 2750 



. _ . QK.2760 

C.. 

DIFFUSING corner 

OK 2770 



OK 27S0 

23«0 ‘^CHECK(N,n) ■ 1 

DK 2790 

nchECK(N,i 5^ a \ 

OK. 2SOO 

NCHECK(N* 16) ■ 2 

OK 2610 

NCHECKi-n. 181 J_2 _ _ _ _ _ __ 

DK 2620 

GO 

TO 3000 

DK 2630 



nK 26«o 

Cat 

nlFFUSlON 

OK 2650 



OK 2660 

Cat 


DK 2870 

^ a a 

diffuser . . 

OK 2860 

c!! 


DK 2690 

6J ■2 

DK 2900 

GO 

TO 3000 

OK 2910 


EXIT-FLOW KINETIC ENERGY . 

DK 2920 



Cat 


OK 2930 

2450 '“CHECKfN. 91 6 0 __ _ 

DK 2940 

GO 

To sobo 

OK 2950 



QK 2960 



c.. 

suoden expansion 

DK 2970 



OK 2980 

2460 MCkECK(N, 9) ■ 0 

»K 2990 

CO 

-10^ 3400 - - - - - - 

- - OK 3000 
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HONIVCOMft THIN FLOW 8TR*I0HTENER8 


C.t.t* 

£•,. tv 

c • • 



C » • 



2510 NCHECK(N,12) 

j^ HECK(N ^nn 

NCHECK(N,16) 

GO Tft lOftn 


• 1 
4_1 

■ 2 


Cf • 

airfoil thick flow STRaICMTEMRS 

Hio NC HECKfN^tZ^ m 2 

NCHECK(N,H) 9 I 

fid TO SOOQ 


c • « 

. • perfo r ated PL.4TE 

c • » 

25SQ wch ECK(N^ 9 0 

nCHECK(N,14) 9 I 

BQ T .P gftiLfi 

C • « 

-JL,., HE8h SrREEN 

C t • 

2S aO NCHECK^N, 9i 9 0 

NCHECK(N,1J) 9 I 

NCH ECKfN,14l) ■ 1 

NCHECK(N,16) « 2 

60 TO 3QQ0 

Cat 

- Cij Internal structure 

HfcO NCHECKfN. 6) 9 t 

NCHECK(N, R) 9 0 

NCHECK(N.12> 9 1 

nCHECK(N,U) 9 J 

^HE CKfN.l T) ■ 1 

NCHECK(N,20) ■ J 

SQ I fl 

c » § 

P lX tOa KHQWN LgJAL-L.088 

2570 nCHECK(ni R) ■ 0 

NChECK(N,17) 9 I 

GO TO 5000 

■CI!!!nulH^LE.DUCT 8PECULI2ED INPUT REQUIR E MINTI PEtlNITIQNS 

C t * ( 


dk 8010 

OS 8020 

OK 3030 

OS_3000 

OK 3050 

OK -30b0 

DK 3070 

. OK 3080 

OK 30R0 

-OR .3100. 

OK 3ll0 

OK 3120 

OK 3130 

.-OK 3l«0 

OK 3150 

OR 3160- 

DK 3170 

DK 3180 

OK 3190 

OS 1200 

OK 3210 

OS 3220 

OK 3230 

DK. 3240 

OK 3250 

DK 3260. 

DK 3270 

OS. 3280 

DK 3290 

-OK. 3300 

OK 3310 

DK 3320 

DK 3330 

OK 3340 

OK 3350 

UK 3360 

OK 3370 

QK..Jifti) 

DK 3390 

OS 3400 

OK 3410 

OS-J420 

OK 3430 

DS_14«Q 

OK 3450 

OS. .3460. 

OK 3470 

O.K. J48Q 

OK 3490 
OS -.35 00 




corners and turns 


IHSICI •! 


NCHECKCN,!^) • 1 
NCHECKfNplh) « 2 
NCMeCK(N, 16 ) « 2 


OK S 510 
-_OK_ 4 Saa 
OK }S 30 

(MC tKilft 


OK iSSO 


C I « 

_i 7 io N CMecK(N ,iai--i-i- 
nCHECK(N, 16 ) ■ 2 
GO TQ loop 

0 1 • 

QlfiUll^ 

2^ 0 NC HECK(N.H) 9 1 
NChECK(N,i5) ■ \ 


NCMECKcN, 18 ) ■ 2 
GO TO 3000 


t DiRRm iQN 

C.. 

. -L.4-* 

C,, OIRFU 8 E* 

^S'ao NCHE'c^rNji^r^ 
M- 10-300.0 

C i • 

_i:« SUDD EN t 


OIRFU8ERS 


ISfcO NCHE CK/N, 9) ■ 0 
GO TO 3000 

" c . oRfyETFAN syTtYm” 

. A 

C » t 

..C„ .Fa-^-AN-NI, 

c.t 

2Rl0 NCMlihLU*_5j_.i_2_ 

NCHECKfN, 6 ) ■ 2 

NCHECK tN.12) ■ 1 

nCwECK(N, 13) ■ t 

NC^’6CK(N,14) 1 3 

ncHECK(N,17) > 1 
. - nCHECK4.N,19J ■ 2 


.FjL^-ANjyyL aR_J 3 UC 14 83 
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..JJAliCK. 




nCHECk(N,20) 

a 3 

DK 4010 

60 LQ 3000. 


HK 4020 

c*. Fan 

CONTRACTIONfS) 

OK 4030 

Cf . 


DK aoao 

2920 NCHECKfN, 5) 

■ 2 

OK 4050 

NCHECKfN, fc) 

• ? 

DK 4060 

NCHECKfN, 13) 

1 1 

OK 4070 

NCHECKf.M,14L 

a 5 

- -D-K 4060 

G'J TO 3000 


DK 4090 

Cff 

DIFFUSERfS) _ .. 

. _ OK 41-00 

c * • 


OK 4110 

2940 NCHECKfN, 5) 

■ 2 

.. OK 4120 

NCHECKfN, 6) 

a 2 

OK 4130 

NCHECKfN, 11) 

a 1 

. pK 4140 

NCHECKfN, 14) 

P 3 

OK 4130 

NCHECKFN, iAl 

a 2 

_ . -OK. .4190. 

GO TO 3000 


OK 4170 

^ H 91 D f f 


DK 4160 

C.. flow 

OnSTRUCTIONS 

OK 4190 

_£.aa 


—OK 4200 



iNTgBNAL STBUCTUSe 


OK 


4210 

J220 


NCH ECKfN^ ■ g 


NCOECK(N, 9) ■ 0 
h.CHECKf^..lg) ■ 1 
NCHECK(N,1<o ■ 3 
NCHEC KfN.1T) « 1 


C » t 

-Caa- 


NCHtCKfN,20) 
GO TO 3000 


■ } 


C * • 

2’je. 


FiXEOp KnQHm L0 88 


NCHE CKfN. 9) ■ Q 


• 1 


NCHECKfN, 17) 

cIlIIiNtlcER input error. chick AND ieTfl’NGllF ERROR FlkS 



jooo 


OK 

_.0K 

OK 

--0K 

OK 

OK 

...QI4 

OK 

..QK. 

OK 

-OK 

OK 

-_QK- 


DO 3001 I m\ U 
If fNCHECKfN.n 


1 8ERR0R ,CT. ,.5 
IF fNCHECKfN,!) 




5 ) SERROR a . 
,Ka. 1 .ANP, 


I® 


,E Q, 0) SERRQR a ,2. 


4230 

J240 
4250 
4290. 
4270 
.4260 
4290 

4310 
4,12a_. 
4330 
_4340 . 
4350 
.4.390- 
OK 4370 
OK .A3iO- 
OK 4390 
i)K...44M. 
OK 4410 
DK 4 426- 


3001 CONTINUE 

cIlM^oIfiNG, POINT INPUT ERROR-CHECK ANF'srTTINC OF ERROR FLAG 

C » > I 


00 3002 1 a 5.20 
J.F_t.NeH g CI<lN'|) ,N E . 


OK 

J5K 

OK 

-OK. 


4430 

JLftilO. 


4450 

4490 


1 SERROR ,GT, *,5) 

- IF . (NCHE CKfN , I) , EQ . 


,LT, ItEaO .A N il, 


.fi._A-A.Nt* _ .AJ.8.{ LN.tAlA.{.N *.I )4 . 

SERRQR a «,5 

L-aMJ) *_A.B a t IMATAc&U-U-X-flT. -1*.E«9) 


OK 

-OK. 


OK 

-GK. 


4470 

_M6fi 

4490 

-4500- 


87 


DtTiCK 


JUiii — W. 


1 SIKROR • *2$ 

3002 CON T I NU E 

C..M* 

iNVtLtD ggC TlDN SMiP E CMtCK tND MES8iBE 

c • • 

300 3 If f(fI8MAPl»ntfl8HAPl»2l>fl8HAPUI)) Q> go T° 

Iff (UNEcT ,LT, (UNCMX-a)) 80 TO JOOK 

lMfiJL-l_leJLfiJUJ 

WHnE(6»ei|l) ITlTUEflPAGE 

wBiTe^hfana?) 

IE (lu ,CQ, n WRITE(6«8008) 

If- f 

LINECr • 9 

_lMi.Jflllt»» 80 0l) n ,»HA E| 

UNECT • LINCCT«3 

8EB P 0R...a. mlt 

3005 IF (((I8HAP2«t)*/I8M«P2.2)*/ISHAPa.3)) ,EQ, 0) 80 TO 3007 

_ . IF ( U INgcT AT , lUlNENj(^ n . 60_T-Q loafc 

IPAGI • IPAftEfI 

WAl TE(P f5Un ItItU i IPMI 

WRITE(6,8007) 

IE (lU .EQ. 1) FPITEtfc.SO 0 81 

ITfliTTlQ, 2) F«ITE(8»8009) 

UI NEC T ■ 9 

3009 MPlTEf9*800S) N»I8HAP2 

l.IJlE^T_t UNECt*3 

8EPP0P • ,2t 

__3.Q.0,7_ if fOtPPOP .6T. .lo GO TO 3009 

IF (LINEcT ,l.T, (LINEMX-S)) 60 TO 3008 

IPME-8 IFACE»1 

WRITE(9,8U1) 1TITUE»IPAGE 

WRllflOjeMli 

IF (lU ,E0, U WRITE(9»8008) 

2) wRITEf6i800iI 

UNECT ■ 9 

. 3OJ08. wR I TE (» ,80 05 ) N 

WINECT ■ LINECT+S 

3 009 IF (TLIfT agBPOtn TLigT 9_8 £r«Q8 

return 

_c 

c************************************************************** 


ctt*ftiftU6’fAifi£Afi6ft OF INPUT Data with annotaTiong aG rebuireo 


/ « 0 0 0 . 61 J_ JP15X* 1 

NRlTE(8«ei00) ititleupaoe 
- - -FRiiitA^imj 


OK «510 

OK 8520 

DK R530 

OK .GSR-q- 

OK 4550 

QK 4380- 

OK 4ST0 

pX 8 t 8 0 

OK 4590 

OK-i4aO- 

OK 4810 

OK 4830 

QK 4 A4fl- 

OK 4850 

04-4880- 

DK 4870 

04-4860 . 

OK 4890 

OK. -47.00- 

OK 4710 

nK 4720 

OK 4730 

J14 .4740 - 

OK 4750 

OK 4T.8(L. 

OK 4770 

QK. 4760- 

OK 4790 

DK 4800 

OK 4910 
OK 4620 
OK 4830 

OK 4640 

OK 4850 

OK 4680 

OK 4870 

— - OK 4980.. 

OK 4890 

- OK 4900 

)***«* OK 4910 

OK 492(L. 
OK 4930 

DK 4940 

OK 4950 

J}K 4980 

OK 4970 
DK 4980 
DK 4990 
. OK 5000 


88 


DAIAC«..._ 


Fa6£ U 


IP au .ea, n wRiTc(6,»iiOt) 

I P-Xlu _ »£0«-JD-XRlTR(6»8LILlj 

wHITE(8,810<*) 

cTTr^aeWnloNnor ha’stPTata inf'oW*tiYn"7no~fo»mat" arrays 



lNTm.H-lNFOMAllQN 


c « « 


__ia£MnlX_U--*-MtE.tl 

MRORMTf 2) ■ ICOMMA 

_ ■ngMMl(13)_l ISPJLfU 

MPORHT(Jo) M IRIOHT 



log ■ 3 

on 0007 1 ■ i,% 

IF (HCHECK(l) ,Ng, 0) GO TO flOOl 

_ .MJRIlE4J-QVJ_A.-iaiM2 

iWRITF(IOV+n • IBINKO 


.laJLGfiS 


’^WRITE(IOV) . 1M861 

M xR ITEIU^iH i l i t 8 G 2 

GO TO 400$ 

_JLO.0 1 IF <HCHECK(I) ,EG. 0 MDaTA ( IT 

HWRlTEflOV) B MOATA(I) 



tE Q , 0) GO TO 4002 


Ijjjfj- P IQVa I... 

MFORMT(IOF) B IIFI.08 

iiFOR HTfTOF»n ■ t8F*CC 

60 TO 4006 

4002 I F fMCHECK fl) ,NE, 1 ,0R . WPATAfl) ,NE, 0) CO TO 4001 
MwRITiflOV) a 1W86S 

HjRITE( IO YAl) t IH8 Q 4 

GO TO 400$ 

4001 IF fHCHECKf!) .NE» 2 .OR. WO AT Afl) .NE, 0) fiO TO 400 i„ 
MPRlTfflOV) a IHGG5 

*^4RITE(10V*lj a 1M8G6 ... 

60 TO 400$ 

4004 HWRITE(IOV) a 1(«8G? 

MwRITE(10V*t) a 1M806 

400$ IQV a tOV»2 

MFORKT(IOF) a IAFL02 

*^F ORMTfIOF»n a 1AFLD4 

4006 lOF a 10F*2 

40 07 CONTINUE 

lOF a lOFtI 

lOvT a TQV.l 

IQVR a lOVUl 


e,,,,, 

C«. FLOaTING-POINT INFORRiTlON 




OK 5010 

. OK $020 

OK Solo 

OK $040 

OK $050 

- - OK $060 

OK $070 
OK $080 
OK $040 

OK 5100 

OK SltO 

OK 5120 

OK SIlO 

OK-$U0 

OK $150 

Oil- .51 60 

OK 5170 

OK 5180 

OK 8140 

DK 5200 

OK $210 

OK 5220 

OK $230 

- DK 5240 

OK $250 

. OK 5260 

OK S270 

JBK 5260 

DK $240 

DK ISOO 

OK 5110 
DK $320 
OK 5330 
OK 5340 
OK 5150 

DK $360 

OK 5370 

DK $380 

OK 8340 

J3K 5400 

OK $410 

_ DK 5420 

OK $430 

DK 5440. 

OK 5450 

OK 5460. 

OK 5470 
DK 5480 
OK 5440 
. . OK -5500 


EMWRIT(IOV) « RBLNK* 



^MiXi « ■ • it& n 


IF (ABS(EMDATA(I>) .it* ItC-B) Oq TO 4012 

EMWBITMOV) ■ RWaSi 

CHXlRlT(IOVtl) ■ RH802 

QO TO «012 

ROQB if (MCHECK(I) ,EQ, 0 ,OR, AB8(EMDATA(I) ) ,tT, 1,E"8) 60 


TO 4009 




lOV ■ lOV^l 




OaTA.MAQNITUOE.CONTROUEO FORMaTTINO 
MprORNiT{I0F') ■ IFFUDO 

IF fEHDATAfn .LT. 1000. > HFORMTflQFl ■ IFFLDl 
IF (EMOAtA(I) ,lT. 100.) MRorMt(IoF) ■ 

IF fEHPA TAfl) ,LT, 10^) MFORMTflOF) ■ I F FLDl 
IF fEMOATA(l) ,UT, 1,) HFORMTflOF) « IFFLO<» 
MFORMTf I0F*1\ ■ ICOMMA 


BO TO «01S 

■400-9, IF (HCHE C Kfl) .NE. 1 , 
l^'wpIT(10V) • rMsSI 

E M WR ITjJ 0 V 

60 TO «nl2 

0010 IF fWC HE C K rii ,NE. 2 . 
Ef*WRlT(10V) ■ Rf-865 

£J! “!U1,( IQVtll ■ RW8BB. 

60 TO 9012 


El^wRITflOV+t) ■ RH8G0 

<» 0i2 to y _«_10V*2 

MFORMT(IOF) ■ IAF102 

gf.aFJa.UflEjt.il.*-IA£U1.4 

9013 lOF a IOFt* 

9Q19 co n tinue 

lOV ■ lOV.l 

WRITg(b.wFQRMT) tHWRITetllfl a_L^lQVI>. 

clIlloEFINITIO N OF t e rminati on CQNT.R5t^ _CilDl ^ 

C ( • « 

p ^ _ 

c,, suV^Try Information print 

** ETwBl'ft’T) a Tf<SC2 

Lt_£lPR l!iT-ait*-0.1--£J.WRll C-1.1 auJi?S6L 


-BUmmJlRY PL-OIB- 


. IF. 1 1.PL 0 1 ..*.NE-^ Q y 6 0 _TJJ .90 1 5 


DK SS30 


OK $550 
OK 5590 
OK 5570 
px 5590 
OK 5590 


OK 5610 
QK 5620 
OK 5650 
qK $990 

OK 5650 

aK-596fl- 

OK $670 
OK 5 660 
OK 5690 

OKJ7A0. 

OK 5710 

QK...5220.. 

OK 57J0 

QJL 5790.. 

OK 5750 

0 .K, 5.7J0 

OK 5770 

0K_i7Bil. 

OK $790 

OK,SflM- 

OK 5810 

nJS,-5820. 

OK 5830 

OK-589.0. 

OK 5850 

J)K_lft6(L 

OK 5870 

D.K-5880. 

OK 5890 

JIK .-5900 , 

OK 5910 

OK . 5920. 

DK 5930 

OK 5990- 

DK 5950 

0K_596O 

OK 5970 

-DK 5960 

OK 5990 
-OK .6000.. 



I 


-QATAtK. 


PAGE l3 


C • * • » ( 
JU* 

C » « 


ETwrITC 2) ■ RSLNSit 

IThRlTt Jl. ■ TMlGi 

ET*it»IT( «) ■ RBLNKit 

QP 

IF (IPLOT ,N|, 1 ,AM0, IPLOT ,N6, 3) GO TO 4ol6 

ETwBltX 

£TwrIT( 3) • T"SG5 

|TwR|Tf 4)_J_ T»*?G6A_ 

IF flPLOT ,GT, 2) gT».RlT( «) « TMSG6B 

..GO.-Ta^-^LM.? - 

£TwhIT( 2) ■ T^'SGT 

E T - R 1 1 { -J 3 „l._T «aGa. - 

ETivnlTt 4) ■ T^SGR 

IF X I PLO.L_^I.._ J V-Q.0_.iQ_,4ai? 

DO 4018 I ■ 10,12 

.i T wRITi l) « R8LMKU 

KO To R020 

.tT«RlltJ.Q3,.r T?^S41 

ETi>.'RIT(U) * T»^SG8 

.lTwHlTa20.-P TMiM 

L, -JLN.NOTA.TlQ_USUi.ATiQN. QF_ . IN.PuT LJAT A 


4022 
C4 < « «.i 

C.-^ 


ST » I T ( .5 j_ .•-.iMiQjL 

IF (TlI8T ,GT, GO TO 4021 

b) i; T»< 8Gi2 

ET>^RIT( 7) • T'«3G13 

...SO... ,10 .4022. 

ET*!RIT( b) 9 T^'SGlO 

-£ I*' Ril{_ 3 G U 

ETwrITC 8) « Tt/SG2 

' HI N 6 A NO Te L 0 C i t Y .“0 Pl'l^l Z A Te 0 U E 8 T 

,GT, 1.E.6) ETHRfT’(~eTT~T^GT 
...t Tw RJ_Lt <») ■ Th S q2 

’ NEXT,CA S e.R £ TgRN OR T£R H! NaT10N REQ UEST 

J-L-LTRlIR n ^ gT. tLT.W-PJ.Ii-.ll_»_Dl8&l 

wRITE(ft,8l05) 

JlRJ Tg^b.aiQ b\ ItTitiRIT 


Ladings fo r listing of section data inputs 


WR1TE(6,8100) ITITL6,IPaGE 

»«RITE (6f6ln7j 

IF TiU tEO, 1) WRITE (6»eioe) 
1 F_t_LU_^£S » .13 . WJII J 9.).- 


iU-lACK 


PiflE ML 



wrITEC 6»8U0) 
LINECT a ?7 


OK 6510 
OK 652A 


00 9039 I a 1,M 

IF fLlNEcT ,LT, <L1NEMX*15) GO 

TM anal 

OK 6530 

nK 659A 


IPAGE a IPAGE+i 
WRITE^6,81 11) ITITLEflPAGE 


OK 6550 

QK hSafl 


K«1TE(6,8107) 

IF flu .FfJt n WRITE f6, 81081 


OK 6570 

DK 65BA 


IF (lU ,E0, i) WRITE {6, 8109) 
WRITE(a,ei 101 


OK 6590 

QK 6660 


UINECf a 19 


OK 6610 
D& 6620- 

C • • t 9 

4023 

6E^MiTI0^ OF Section data information and format arrays 

OK 6630 

.DK....4440_. 

NFORmT( n a NLEFT 

NFOSMT (921 a IRIGHT .._ 


OK 6650 

. _ __ ..H&_.4644_ 


lOv a 1 
lOF a y 


OK 6670 
DS._466Q. 

C • t • ff 



OK 6690 


integer information 


DiS_.4liUL. 



00 9030 J a 1.9 


OK 6710 

DK...412CL_ 


IF (NCHECK(I,J) ,NE, 0) GO TO 
N^PITEflOV) a IBLNK2 

9029 

DK 6750 
DK 6790 


N‘'RlTE(IOV + l) a IBUNKa 
IF fNOATAfl.J) .Ed. 0) 60 TO 9028 

OK 6750 

. OK_iL74jOL 


M«SITE(I0V) • IM8GI 
N J.T E t tOV*ll a 1MSG2 


OK 6770 

. .. .DJi.-4.760... 


GO TO «t28 


-.iGU- 


aO?« IF fNCHECK fl, J } 

nwr:TE(10V) ■ NDATA(I»J) 

.IJ3V...1 _1qVjtJL 




NDATAflfJ ) ,EQ, Q) BO TO H Qii.. 


NPORMTcIOP) ■ 1IPL05 

5iE,0»!iTxlOPilJ_i_15PJC.t,, 

GO TO «fi29 

CHECK f l.Ji .HE. I 


.OR. ND ATA<I.J) .nE. Q) O.Q. TO. JoJ^ 


NwsiTeaov) . IM8S3 
'lial Tlll!^*n ■ IHsOJ - 

GO TO Afi29 

a 01b_lf fNCH E CKfl.J) .NE. E .OR. N OATA(ltJ) .NE. Qt GO T 0 « fl J 7 
MKHITeCIOV) * IM 8 G 5 

K^P 

GO TO «02ft 

"J.2L 


NwRITEclOV^l) ■ 1K8G8 

« 0 2 B_ 1Q.V . V.*JL 


NAOPmT(IOF) 
MfORMT(_IOF»U 
U029 lOF a lOP+2 
R030 CP^tINuI . 


1AFLD2 
I .lAltO<L_.. 


OK 6790 
DK .69Qi 

OK 6610 
QK 662il. 
OK 6830 
. JIK 66.8iL. 
OK 6850 
J3K-6660_. 
OK 6870 
. .OS. 4610 . 
OK 6890 
0.K..6 9 CO- 
OK 6^10 
OK 6920_ 
OK 6930 
UK 6990... 
OK 6950 
OS 49.6.0 
OK 6970 
DK 6980 
OK 6990 
.DK 7000- 


92 


-OAlAfllL 


. FaSE. iS 



lOVl ■ IQV.I 

lOvfl • 10V1*J . 

OK 7010 

.... . _ DK 7020. 

fill 


OK 7050 

* _ * 

fi QaTInB»PQInT INFOBmaTION 

OK 7090 

** DO 90X6 J 9 5.i0 

DK 7050 
.. . . _ OK 7060 


IF (MCHECK(I,j) ,N6, 0) 60 TO 9031 
EnwRxT(tOV) • «8i,Nli2 _ _ 

OK 7070 
. . .DK J.QfiO 


ENWRIT(iOV^l) ■ RBLNK9 

IF_f_AB8fEvOATAfI, Jn l»F-6) 60 TO 9036 

OK 7090 

... .- ..oK. ,70.00 , 


EN«RIT(I0V) t RM8G1 
ENwRlT(IOVfl) ■ RM8G2_ __ 

OK 7110 
.. DK 7120 

9051 

60 TO 9036 

IF fNcHECK^I.Ji .EQ. Q ,OR. aB 8 f ENQATA / I . Jn .LT. 1.6*6^ 

PK 7130 

OK .7190 

(56^T0^9r>XX’ 

DK 7150 

. .. , JJX..7160 


lOV • lOV+1 

IF fj ,NE, 5 ,AN0, J ,NE, 6i 60 TO 9o32 _ . 

OK 7170 
- .JlK 7160 


NrORMT(IOF) I IFUOO 

OK 7t90 
OK 7200 


BO TO 9037 ^ 

OK 7210 

J3K ...7220 . 

• • 

9032 

OATA.MAGNITUDE.CONTROt.LED FORMATTING 

DK 7230 
DK 7290 

NFORMT(IOF) a IFFLDO 

tF {E^DATA^J.J) .LT. JOOO.) nFORmT^tOF) a tFFlDI 

OK 7250 
.DK ,7260,, 


IF (ENDaTA(1,J) ,LT, 100.) NfORMT(IOF) a IFFLD2 
tF (ENOATAfTiJ) .LT. 10.) NFORMTaOF) a jFFldJL 

OK 7270 
DK 7260 


IF fENOATAfl.J) ,LT, 1.) NfORMT(IOF) a 1FFL09 
NF0RMT(10F*1) a ICOMMA 

DK 7290 

OK...730Q 

9033 

GO TO 9o37 

IF ,NCHECK#I.Ji .NE. 1 .OR. ABSfENDATA^I.Jii .GF. 1.6-6^ 

DK 7310 
DK 7310 

t GO TO 9039 

ENWBiTtTQVT ■ BH8GX 

OK 7350 
DK 7390 


6NwRiT(IOV+l) a RM8G9 

GO TO 9aX6 .. 

OK 7350 
... DK 7160 

UOlU 

IF (NCHECK(I.J) .NE. 2 .OR. aB$(EN0ATA(Ii J) ) .GE. iiE«6) 
L BO TO 9ft35 

DK 7370 
JJK. -2120 


ENMRIT(IOV) a RM8G5 
.ENwRiTdoVAl) a RM8b6 . 

OK 7390 
__ OK 7900 

903$ 

CO TO 9o36 

_£NwR_lTfIOV) a RM8CT _ 

OK 7910 
.0i_7920 

9036 

ENWRIT(IOVaI) a RM8G8 
_tOV f lOV+2 

OK 7930 
OK. 7990 


NFORMT(IOF) a 1AFL08 
-NFQRMTf lOF+ll a IAFL09 

OK 7950 
DK 7960. 

9037 
-M3 6 

lOF a IOF+2 
CONTINUE 

OK 7970 

. _.aS M60_ 


lOV 9 JOV-1 


OK 7«90 

m< 7<?A0 



LJNEcT ■ LINEcT+2 

<>019 CONTINUE 

RETURN 

r _ _ _ 


C..DATA error NEsSaSE FORMATS 0« 7550 

C,». Q iv 7560 

6000 format (/10X,62HCRITICAL OMI 8 SION(S) IN TUNNEL MASTER CONTROL DATA OK 7570 

A PREVENT EXECUTION OF THIS CA 8 E,yloX,tt>HANY SUCEEDINS CASeS N ILL N OK 7560 

SOT BE affected,/) DK 75r0 

6001 format f/tOK.TSH** ERROR T NV aLID TEST 8 ECTJQN UR SIRl AM END 8 MAP OK 7EAfl.. 

aE code ^as SrCCIFIED As»I2«23H ( 3 HOULD RE 1, 2 oR S),/ DK 7610 

S 13X 29 H THIS CASE CANNOT B E EXECUTED ',/ 1 p K ..- 7 -6 20 

6002 format (/ 5 «h THE UNITS OF MEASURE CODE IS IMPROPERLY SPECIFIED AS OK 76J0 

_Ai_IaAJ5.tt-H ISMOUI D aE 1-0R .-»T-, -.-CMEC-K M A5.TE8. £ARD_l£.QLUMii _4 Jiii._liiitO- 

P 120H sEE the. data tABuI-AtION aT T^E END OF tHIS CASE, THE INTE OK 7650 

CRNaTIONaL SYSTEM OF UNITS HILL BE ASSUMED FOR THIS CASE. ...yj. Dil_2AAa._ 

8005 format {/I3,2X,55H** ERROR INVALID SECTION SMaPE CODE waS SPEtl OK 7670 

AFyED A 8,i2,SaH (SHOULD B.E li 2 _a 8 . .il,.., .IH|,S._2E£1I5N. .MILL 3£ 8 KlPPE-.M_.mo_ 

80 ,/) OK 7690 

eooR format (/n ,2X ^l0 9 H«» ERRO R INPUT SECTION TYPE CODE fCARD COLUM. QK ..llfttL- 
ANS 1 AND 01 calls invalid section type, data card IUNORED, **/) DK 7710 
6005 FORMAT f / 1 3.2X,111 h*a ER ROR «»■ £R I TICA L OMlSSlONfS) IN SECT ION.. 

AUT DATA. SEE data TABULATION AT EnO OF OUTPUT FOR THIS CASE, *•/) OK 7730 

Sfo fe format t lHl//5y ,Al, l9A<l,A3i28X,qHPAGE,Ii) 1 DK ..17R0_ 

8007 format (120H NO, SECTION TYPE SMaPE Ml mI.DI?’' AHEAl OK 7750 

A Aj/An aR. CR s theta V^ MAMi_ LEN STH nP/jilL DP/OO pK TlfcO 

a / 30y, 3SH H 2 w2,D2 A,EA2 A2/A0 ,17X, OK HTO 

C—2.2J5! y 2 M A c H 2 ) D K , .7 7 8 0 . 

8 ooB pORMAT(26X,26H METeRS METERS 8 Q M, 17X, 32Md£0WEE5 M/SEC OK 7790 

A meTeRS / 1 2 qH A ♦•••♦ DK .7600 

R ^mmm^ ^mmm^ ^mmrnmm^ ♦•••+ +■■■•♦ ♦•••♦ OK 7810 

C J ... . DK...782Q 

6009 format (50X, 26H FEET FEET 80 FT , ISX , 31H0EGHEE8 FT/ pK 7850 

ASl£ FEET /120H »A A”.?.-.---?-?-* *•••* _* DK.746Q._ 


OK 7510 
OK 7530 


DK..17R0-. 
AREAl OK 7750 
DP/OO pK. 7IfcO.. 
OK 7770 

- DK ,7780 

M/SEC OK 7790 
DK .7600 
♦•••A OK 7810 


A»«»A A*"""A A"*"»A DK .7600 

.•■■A A"“"«A A*"»A OK 7810 

DK 7820 

, 1SX,31HDEGREE8 FT/ pK 7850 
A^M^A A.**A-.-. . -* DK .7JI.6Q-. 


eooR format (50X, 26H FEET FEET 80 FT , ISX , 31HPEGHEE8 FT/ pK 7850 

ASl£ FEET /120H AA A”?.--- -?-?•* A-M-A A.*-A..... _* DK.7 JIl6Q._ 

a***a a***a ♦•••♦ a"""*a a****a a**"*a a ok 7650 

C»»»A ) , DK .7860 . 

C, DK 7870 

c , , annotated Tabulation l ab ellin g and qaTa. formats. . dk 78tO- 

C,,. PK 7890 

el 00 fo rmat rlHl/ /Sx.Al.l9A a.A5.6x.l6H...CQNTINU£c...^,6.X,9HgAG£.1.3/ZA DK .7.9 Oil-. 

4 a«x, 52 m annotated input data Tabulation// ok 79io 

-S 1 2oH lEMP TYi IN OICtTES OPTIONaL, NON.RESUIREO INPUT PARAMETER M.DK-.7m_ 

cas seen omitted or Parameter may be intended as zero, / ok 7950 

r 1 20H tER ROR i indicate s H andaTQRY INPUT PaRaMETER HaS .SEEN . DMlTT._M_7A^^fl_ 

EEl), this must me corrected bEFO»E computation is possible, / OK 7950 

.. . . F 120H «EyT RA» indic ate s su perfluous I.nPUT. parameter HAS BEE.n UNN DK -7-96.6- 

GECESSaRILY included on input card and M,y bE REMOVED, / OK 7970 

H l2oH 'QPT'N‘ INDlCAlEa-OPTIONAL INPUT DaTa MAS.BEEN OMITTED AND .DK. -7166- 
I The PaHaMeTeR wRL DEFAULT TO a PREDETERMINED VALUE, //) OK 7990 

8lol format {a7X,j6MTllNNEL. Master CONTROL DaTA// ..QK..KDDD_ 


94 


PaTaCK 




A 20X.78H CASE TUNNEU UNITS SECT* SECT, HI »<l HOOEL VO PC OK ftOlO 

Bwffi. PI II P- ATH-..^ OK.-AX»M- 

C 20X,SIH 8EQ, TYPE INLET EXIT 01 BLK6E LE OK S030 

OVlL—l OK.lIUO 

eio2 FOR»>AT (2tX*4H i^0.,13X,12H SHARE SHApC, ISX, 17H P£P» HE6A*/ OK SOSO 

A . SQX.17H M H t£NJ..JU4£iL JtAlIS ATH DEC C ATM > OR-fiOtO 

SioS format (21X,«M N0,,13X,12H SHAPE SHAPE, 13X, V7H PEP* FT/ 10(3)/ OK SOTO 

A . S a X Ai TH FEET ■■FeEI__£ENj: 4EC HP aTW DEC P aTm ) DK-AOSO. 

8l0« format (1»H data F^ElD BEGiNS/ OK 8090 

A Iftfl ..IN .Card CQLUHN #Ji. l 1 a S S p nx s inft 

SIP 2l 28 3l 38 41 48 / OK 8110 

C._-APX,7S H M * A * * ». w » -* »*■ *- ♦ A«>8 ♦»»» ♦ 8»_QlL-ftilO- 

0»"A ♦•••♦ ♦•••♦ ♦•••♦/) OH suo 

_8lAS_FgRMA T f//44X,32HgA$E_TO_Hl N ATI QN CON DITIONS DATA// P K Slap 

A 27X.68H CASE TERHINaTION OCCURRED (DUE TO BLANKS IN CARO COLUMNS OK 81$0 

B 1 AND a^/ nK i«^A 


C 2 SX. 8 H after, 13 , 81 H INRuT SECTIONS, AND ACCOrOING TO THE FOLLOKZ OK 8170 


ONG-XO-NDIT 1 QNJ 3 . j.«// 

E 20 X. 81 H summary 

m oc se tu sn— 

6 20 X,S 1 HCHArACtEr 1 STXCS 


I 20 X. 81 H OUTPUT OF LENGTH TABUL 

J.XEA_P 0 H£Rl Case ^ 

K 17 H TERMINaTION.COOI/ 

L laH Data field ISyllH CONTaI N ED IN/ 

M 18 H CARO COLUMNS I OX, TOhSaiS 7 -B 

?ilO ^ LU 15 L 8 jl 2 JL( 

/l 

Slo 8 FORMAT ( 28 x,A 4 , 10 X,SA 4 , 5 x, 3 A 4 , 10 X,A 4 , 13 x,Aa/a 0 X, 3 Aa) 

eio 7 FORMAT /asx^aaHBECTION DESCRIPTION DAT Ayy 

A 120 H sect; sect, SECT, SECT, TOTAL ITEMS Ml w 

8 W 2 L/ 5 H. CHORP BLKGE PHI KEXP Pn RNREF. E 


PL 0 TTIN 8 AS 

A Function 

OF LENGTH 


C 120 H SEO'. type inlet EXIT NO, PfR D 

0 02 S/AL DHU 8 PRSTY KME 8 H K RUFNtS 

E ( 20 H INPUT shape SHAPE DUCTS OUCT 


INPUT 

DATA 

TABULATION 


eio 8 FORMAT ( 87 X, 42 H M/M, 
A 87 X.S 3 H 8 Q M/ 


B 24 X, 98 H NO, NO, 
CENT PEG 


SloY format ( 87 X, 42 HPT/FT, 
A 87 X. 83 H 8 Q FT/ 


B 24 X 98 H NO, NO, 




10 ( 8 ), / 

PvRa »l 


BO M 


10 ( 8 )# / 


DK 8390 
nw aana 


C DK B 410 


OK 8430 




FEET FEET FEET FEET FEET SO FT FEET C OK B 4 S 0 




i8i.vrr.Ti 




0 iSOH «« « 4 #•••« ♦•••4 DK 9510 

E ♦»••» ♦■•»» »»»«♦ 4»w* »#■■» 4»»«* 4»»»» *•■»* *«■ «» 4»w* /) DK iSlD 

em FORMAT (1H1//5 x,AI, 19A«, A), 9X,UH,,, CONTINUED, «M«9Xf4HFA6E>n//) OK «9}0 




-LT 


--IPiEfl. 


P^B£ _4 


SUt^ROUTlNE SPECO(«f*HACH,V) 


e THIS RO UTl»jE, A SUaROUTINg OP THg > 
C LOCAL HtCH numbers aNO VELOCITIES, 

C************************************! 

COMMON/BLOCK' 


80 10 

LftA*j»J-*i*4.t*AAAA.iuti[ *0 —20 

r******************************** 10 30 

:N PRQSgAM PEBrO BM, eOMPUTCg THE BO ,«0 

80 50 

>A*_« «J>A M->AAAAOJ> *A8 J AA-iLOAA ».» Ut H .80 80 


i 


I EMN • |mT.(A8TAP/A*(2 /(6«1’ ))»*((8A1 

J fl.»f6»l.W2.*gMT»«2lAAf fS *ltW2./>6«l .n»gi4TW 

2 ( ASTAP/A«CMT*(2«/(64l t ) }••( (3,*6)/2,/(6«l I ) )* 




X 4«AfX>lXXjn i 


MFrWllrffl 


IF (ABS(EMN«CNr)/CiiN .|.T, l»E-<») 6c TO 2 

EMl-i EHP 

r.o TO 1 

AW aCH ■ gWN 

ifrLbymb^Ui 


KM 


V m AMACH«a8L 

PgTUPM 

INo 



- -.OUIi>UT £AGC- 1 


SUgROUTjNE OUTPUTf OT 10 

C*«i*********!it*****(k********!t*]k *««««**«* Ol 2Ji 

C*******a******a*******a***a«******aa****«***a************a*«**«**«****« qT 30 

£ THIS Rnyri>£, i IMEiJyT.lJ!<.£_M JME-MaM, PROsRam pgRFOffM, han ol^es t he nr Jo 

C OUTPUT formatting OP CaLCULaTeO SECTION PeRFORMaNCC INFORMATION, OT 50 

C*AA.******l*****A*A.****t**J-4A±A*AAA*»tA**A*AAAiAAA*.tJA#*JA.A,tAt.AAAAtJlA-*.4a4.-0T &Q 

C*********************************************************************** nT 70 

COMMON/StOCXA/lSfeO.ISHlPl^JSHAPg^N ___ Ql _ 

COmmoN/RLOCKE/ amaCHi,aRaCH 2,AR,AI, MOA(,,A2,A2b40/Dl»O2,EK,EKO,cl, OT «0 
L_ Jjl«.H8«.T.M2«.jLl ai l-OiL 


OIMCNSION NSECTf i80),NSMaPE(3) 

bEFrNYfibN's 

tiTi^rs ^ rr Y) ^Tse 


~cTrTrs^6HoN->YPE N A mE 


data mSECTI 
_L- AHtElI 
UATA n 8ECT( ft),NSECT( 

.1. RITlSU-ftN-SEC^RHU 
OaTa NSeCTf ll),N3EcTf 
.4 . jtjj!lOfl£.»_,<AHk_lN,.-JLH-T. 

Data vsecti is)»!v,sECTf 
_L « HQPE N^ aH#THR , aHOAT 


7)»N8ECT( 
5-a— RHl-Pa-M^ 


3) »n8ECT ( 

?M A/ 


«) 


e)»NSECT( 


9) 


12>,NSECT( 13),N8ECT( l«) 

£!.♦ tt HT - 8E ,_..2HC I/. 

17)»»,sECT( l0)»^3ECT( 10) 


ithT- 3E^_2MC1/ 

DATA NSEcT( 2l),N8EcTf 2S)»NSEcT( 23),N3EcT( 2«) 

.1 . iHiON3_,__ilflT,ANl_,_iaM_.ARE,._iLHA-J^^ _ 

DATA N 8 ECT{ jfe)»N 8 CcTt 27 )»N 8 EcT( j8),N3ECT( 29 ) 

1 R^CQnT. R>< RACT. «Mis. 3, aHlNG. , 2HE / 

data NSECTf 31),N8ECT( 32),N3ECTf 33 ),nsIcT/ 3«) 
_1 gnTURN. ttHING , tiMVANE . . jiHi 


data MSeCTf 36),.v8BCTf 37)»M8ECTf 5 «),n8ECT( 39) 

A RHcQRN, UHER a, RMITM , «HVA NE 2M8 1 

DATA N3ECT( gi),N3ECT( « 2 )*N 8 EcT( 43)*'^8Ect( 44 ) 
J._ RHCqRm. 4HER. . 4HaiQ V. 4HA]siE8a 2H ./ 


data NSECT. 46),N6ECTf 47),NSECT( 4e),N8ECT( 49) 
4 4HCRNR, flH.OlP. gM 8 N t. flH VaN. aHES / 


52)»v8ECTt S3 ),msECT( 54) 

-t-JlH 

data v,SECT( 5S)»f^8ECT( 57 ),m 3ECT{ 5d),N8ECT( 59) 
!■ 4.HEXIT, oH KIN, rHETIC , 4H ENR, 2MGV/ 


data k,sECT{ 51),k; 8ECT( 
_1 4 HPIFF, 4MUSER, an 



N5ECT( 5)/ 


N8ECT( 10)/ 


OT 110 

OT 130 
_ai — uo. 

OT 150 
_£J — 140- 


OT 

- 01 - 


170 

AMU 


n8ECT( 20)/ 


OT 

-Hi. 


OT 

-Ol- 


190 

-iOJL 


210 

-120- 


n8ECT( 25)/ 


N8ECT( 30)/ 


NSECT( 35)/ 


n8ECT( 40)/ 


OT 230 

Ol 24JL. 

OT 250 

01 230- 

OT 270 
QJ 2iJl- 


OT 

-Ol- 


290 

lOIL 


N 8 ICT( 45)/ 


OT 

-HL 


ilO 

-Ufl. 


N8ECT( 

50)/ 

OT 
. OT 

330 

340 

N8ECT( 

55)/ 

OT 

350 



aT 

160 

N3ECT( 

60)/ 

OT 

370 



01- 

sao 

N8ECT( 

65)/ 

OT 

390 



nl_ 

500 


99 


o jO oin j n in r> 


nuieuT. 


»- 


04tA n8ECt( «n#N8ECT( 92)»N8 ECt( 93)»N8ECt( 9«),N8CCT( 95)/ 

^ gHSING, uHLg r, aHlXen, aH L08, 2M8 y 

data k 8ECT( 9fc),N8ECTf 97 ),n8ECT( 98),n8ECT( 99) , N8ECT{ 100) / 

-.J <lHHyLT.» »..8HCN8T , .2Jl..iu/ 

DATA VSECTf loll* NSEC T(lo2)»NSECT(lo3)»N8ECTcto9),N8ECTn05)/ 

I --A HMUyT, <Ltf D y c» . <t M.T C0,-9 HnT R A, ihCU 

data NSecTfl06),N3EcT(l07)»N8EcT(108),N8EcT(l09),N8BCT(llO)/ 

1 ah mult, an oue, < 1 HT T,, oh van, 2HE8 / ^ 

Data N3ECT(iU)»N8ECT(n2)»N8ECTni3)»N8ECT(ii<l),N8ECT(U5>/ 

_..J tMWyj T , «M .D. 1 ,. AM » h AL^flHL C R» 

data NSECr(itA)»N8£CTOj7)»A<8BCT(U8),N8ECT(ii9),N8ECTn8j?/ 

. L 4h P i, < t H» H AL|- 9 HL CR, 2 Hnr../ 

data N§EcT(121)»NSEcT(l22)»N8EcT(ia3),N8EcT(l2«),N8ECT(l25)/ 

1 qHw D uM PBN, <tHR, N, AtHQ VA, 2HNE/ 

data NSECTf j26)«N5ECT(i27)»N8ECT(i2e),N8ECT(i29),N8EcT(l30)/ 

1 UHM 0 , UH<»^A, <iMLL aHIP C, 2RNR / 

data NSECT(i3i},NSECT(t32>,N8ECT(iJ3),N8lCT(i3A|)#N8ECT(135)/ 

1 <1Hm D , « H?«k.A, 0. ; 

OATA NSECTf 136),N8pCT(l37),N8ECT(l38),N8eCT(l39),N8ECT(U0)/ 

1 « hhult. uH Due gHT Dl, <lHFF U8 , 2 HeR/ 

Data nseTt ( j « 1 )» nSECT (. 42), N8£cT( 443) »N8ECT( .«<»), N$ecT( 195)/ 

1 qHVtNE . 4H0 pi. 4HPFUS.^ OHgp l >j L 

data NSECT(t46)#NSECTe,47),N8ECT(i48)»N8 ECT( j49)#N8ECT(150)/ 

1 A HSUD . 9HEXP . 4Hrt 0 . aH« SN, 2H6 |l/ 

DATA N8pCTf )51),N8ECT(152),N8iCT(t53),N8ECT(l5«),N8ECT(155)/ 

1 KHSUr^ . 4HEXP 4HW n , 4 H, H 2Hp / 

Data n 8ECT(,56),^8ECT( j57),N8ECT(i58),N8ECTf jS9)#N8BCT(160)/ 

1 4 HFAN , 4MDUCT^ 4M A 8^ 4HTPUT, 2mS / 

DATA N8ECT(i6l)AN8ECT{,62),»M8EcT(i63)»N8eCTTi84)ANSECT(165)/ 

1 A HFAk, , 4HC n MT, UHRACTf 9H1QM , 2>t / 

data M88CT(l66),K.8ECT(167)»N8ECT(16e),N8ECT(189),N8ECT(170)/ 

I _ 4 HFan , 4HP I F8 , 4HP|CN ^ 9 H T R 2HDV/ 

data N8ECT( j7i)»N8ECTf i72),N8ECT(i73),N8ECT(t79),N8ECT(l75)/ 

J. 4NHULT, 4 H I^T, 4HRML , 4H a T RC-».-2MIJLZ 

DATA N8ECTf(76)>N8ECTo77)»N8ECTn78),N8ECT(i79),N8ECT(ie0)/ 
1 4 m hUi T, 4 H!P l , 4H FIXE^ _2Mii/ 


t • • 


qECT iQN END. shape NaHE DEFINI t lONg 


wl • WN*S»4 


9H CIRC»....9M SE CI.»-ttWit. Q/ 


IF flSHAPl ,k;E, I ,0R, I8HAP2 ,wE, 1) 60 To I 

7*wpff^'sTAH^i^T8*^'^P SECtIONj wmICH MavE CIpCul>AR CROSS-SEcTIUNS 

EN -D.8 


OT 510 
___OT — 520- 
OT 550 

O-T- 590- 

OT 550 

OT — 550- 

OT 570 

01-150- 

OT 590 

OT-- 500- 

OT 810 
aT_- MO- 
OT 930 

01-550- 

OT 650 

OT 670 

01 580 

OT 690 

. — Jll Too. 

07 710 

OT._I20. 

OT 730 

OT - 790 

OT 750 

OT 7-60 

OT 770 
Ol. —760 
OT 790 
OT 800 
OT 810 
OT 820 
OT 830 
. . -OT . 890 
OT 850 
-OT 860 
oT 870 
OT 680 
OT 690 
OT 900 
OT 910 
OT 920 
OT 930 
OT 990 
OT 950 
QT 960 


IF-CUYfE ,EQ->-l.). OT 960 

1 1^9176(6,^111) Nj,(ig8ECT(l),I»Nl,N5),w8MAPE(I8HARl),Dl,Al,Al0A0, OT 990 

. 2. Vl^-AMAC.ttl^EL,.E5>150+JJ8.T(AP.£.Xl.8K6.P-2>^52^2^52lLA.fi,V2,A‘tACH2 QT 1000 


100 


OUTPUT 


Pa&£ 3 


OT 
OT 
OT 
oT 
OT 
OT 
OT 
OT 
OT 
OT 
OT 
oT- 
OT 
OT- 
PT 
oT 
OT 
OT 
OT 
OT 
OT 
-OT. 
OT 
-QT 
OT 
.01. 
oT 
.0,1. 
fiT 
-OT. 
OT 
JU 

Cl I H SECTI0 >^8 f-HiCH HAVg.xXlfiCULAR £J?.aS8.<t8E^J0H-Al- OT 
C,,THE UPST»E*M END AND a NON.CIRCULaP CROSS.SECTION aT THE DOwHSTREah qT 

_C.^IND 0_I 

C t « » OT 

If (IIYbL ,fQi. U QI 

1 «RITE(6,9171) N,(N8ECT{I),l»Nl,NS),N3HAPE{18HAPl),Ol,Ai,AlOAO, OT 

iT (lfrBE~.EC, 2) ” ~ ' OT 

1 WRITE(P.91 81) N.t*» 8EC T/I)>l*Ni,HS).N 8HAPEfI8HA p 1 ) j OijAl^LQkQM - 0 T 

2 Vl,AHACHl.EK.EK0,\8HAP£(i8HAP2).H2,w2.A2, A20AO.V2. AMACH2 OT 

l£_i ITT P£_ FftlJ OT 

1 W(?ITE(6,9l9l) M,(\SECT(n,I«Nl,N5),»j8HAPE(I8hAPl),Cl,Al,A10AO, OT 

. .XLT. 
OT 
• QT 
OT 
OT 


IF (IryPE ,E0, 2) 

1 ''RITE (6» 9 1211 .Nt (NS£CT(I),I«Hi.N5)»HSHAPfetI5HAPl),Pl,Al,AlOAtt, 

2 Vl • AHACHl,EK.EK0»NSHAPE(tSHAB2) ,02. A2. A2oA0, V2, AMACH2 

IF (ITY^E ,£Q»-J)-. 

1 «-RITE(b,9l3l ) N, (N8 ECT(I),I«Ni,N5),N8haPE(ISHaP 1 ) , 01 , A 1 , A lOAO , 

2 4R,TH2,yi,4PA_C.Hl.ELtEK,CK0,f*SHAPEasoAP2)»D2»A2,A2QA0,V2,AOACH2 
IF (ITYPE ,eo, «) 

1 i»plTE(6,9132j__N, (NsECT(I), I*'^l,'“5),N3HApltISH4pJ) ,01,Al,AlOAQ, 

2 a«,TH 2.V, ,AHACHi.Et,H8MAPE(I3HAP2),02,A2, A20A0,V2, AHACH2,EK,EK0 

1 IF (ishaTi 

cI*NRlTE*STATEME^8*i^^8icflON5 'which have nonVcircuXaR CROS^^ 

C,,aT PO.TH E.NOS - - - - - - 

^ • 

IF (ITY.Pt_*EJJ.,-i4 

1 '►RITE(6,9i«i) H,(N8ECT(n,I«Ni,N5),NSHAPE(18HAPl),Hl,wl,Al, 

2 *lC*0,yUAf.ACHljeL,EKASMi-N8H4pEclSHAP21i.H2i>(2.,12,A2_aA4tVli. OT 

3 ahaChp 

IF (1tYPE._^EQ*_1) 

, w«ITE(6,9t5i) n,(NSEcT(I),I«Nj,N5),H8maPE{I5maP1),H1,(»1,a1, 

2 A 10A9^y,i^H4.CHIx6.'‘jJJ<i>..JHaHA.P.IlI8|lAP2)4.M2.,,«lU2^A2iLI^!^^ , 

IF (ITYPE ,fc0, 31 

t »««lTEY9r?.l9U-_J^t.Ua-ECILlljl*HltN9.JiNSHApet,laH4PUjMlj_HUil, 

2 AlOAU, aR,Tm2, VI, AHaCh1,EU,EK,EkO,NShAPE(I8hAP2),h2,w2, a2, A20A0, 

J V2, AHACH2 

IF (ITYPe EO. <t) 

1 . oaiTEf 6.91 62) H .rH8ECTrI).l»Nl.N5).RSHAPEfI8HAPn,Ml.iNl.Al, 

2 AiCAo.aR,Th2,Vi, A'^aCH^ ,EU,H8HaPE(I8HaP2),M2,h2, A2, A2UA0,V2, 

3 AMACHi,|,C.£Kil 

SO TO U 

_. 2 IF flSHAPl ,ME , 1 ,0R , I8HAP2 .EQ. 1) GO TO 3 

^ • 1 


.2 A H H2,l(.tjAnACHl D ».M8ttAPElXajlAP 21, H2i J*f2,.A.2 UZ CSASL^i^ - - 


5 AmACHp 

-If- (ii.y PE ...LiU-P-)- . . ... - - 

1 HBITI ( 6 , 9192 ) N, (N 8 ECT( 1 ),I«N 1 ,N 5 ),NSMAPE(ISHapI), 01 ,A 1 ,A 10 A 0 , 
-2 ARjJ>il,-y.|.tA.HACHl,eL,H8 HAPE{I 8 HAP 2 l,W 2 ,.H 2 , A 2 , A 2 DA 0 , V 2,.AMACM2, . . . 


1010 
1020 
1030 
1080 
1050 
1060 
1070 
1080 
IO9O 
1100 
1110 
1120 
1130 
1180 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
U85 
1250 
1260 
1270 
1.2.8 Q 
1290 
130(L 
1310 
1320. 
1330 
1340 
1350 
.1160. 
1370 
ISfiiL 
1390 
1400 
1410 
1420 
1430 
1880 
1450 
1460.. 
1870 
1480 
1490 
1500 
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c... oT isro 

3 IF (iTYpg .EB. 11 -at- l i a o 

r ««ITE(6,920i) N,/N8ECT(1),1«N1,N5),N8 MAPEj OT 1890 

2 aQAQ,Vl»l M AC H t, i L«gK,E t<Q ,Ji8.H4PE(IBMA.P2i>.D2,i2>82a8a».va,8MACH2 OT . U OXL 

If (ITYPE ,E0. ?) OT UlO 

i »DlTEf6.9 2 in N. f NyECTfn«HNl.N5^.Nl>HApEfmMAPl>.Hl.MlfAl f Ot- lUfl- 

2 «1OA0,V1,AhACh1,EK»EK0,NShAPE(18hAP2),D2|A2*A2oA0»V2»AHACH2 QT USO 

IF flTYPg , EQ, 1) OT^UOa- 

1 waiTE(6*422i) N»(N8ECT(n>I«N4,NS),N8HAPE(l8HAPn»HnMt|Al> OT U80 

2 tlOAO.ARfTH2rVlf**<ACHl,EL.EK.tKO,N8HAPE(iaMiPai.D2,A2^AgOAQ .Y2^ Ol-t^WL. 

J AmACMj OT U70 

. IF (ITYPg .EO , A ) Ot-UftO- 

I i»HtTE(6«9222) N.(M8ECT(I)ilaNl,N5)*NSHAPE(I8HAPl)»Hl|Ml,Al« OT U«0 

"'sHk 

a RETURN nT lT2n. 


performance CAWCULATION OUTpU T write P0RMAT8 

9111 FQbH>t t/n. 1 X. AAA, A2 . n, AA. 6«.F9. 2. P1U2 , ft. 2 . 16<. F8. I, P7. 8 , _■ 

* P’,2.2^’.S/25**A9*8X»P9,2.Fii 2,F7 2»16K,P8 1,FT,3) 

9121 F ORMAT (yiLt X, AAA, A 2nY«AAf8y^F9 dl J.LU2dLI^4 1 f-T^At 

* T<,a^«.8/23x,A(|,8x,F9,2,Fll,2,F7,2,lAx.^«tl,F7,S) 

1 ? I FObMAt f/n.n.OAA , A24llu*iU_8x4f?.*.2 *£11*14 lFJ*2*£.’*-2-»£®.l $ 1 UI$ 

A ra a ava Aii ra a via a a mm « mn 


OT U70 

Ot-lAftO- 

OT U«0 

, Ql_ilftO_ 

OT 1710 

01 - 1120 .. 

OT 1750 
_J}I_17«Q 
OT 1780 

0T-.X7P0. 

OT 1770 
_..OI_17«0. 

oT 1790 
.. QT 1800 


A F9 2.2F9 5/23X.Aa,F8 2,F9 2,FH 2,F7 2*16X»F8,1,F7,5) QT 1890 

. .?1 p2 .FON»Lil_ l/lii f il* 1 A * A. A I F 9 • 2* ^9 . 2 4 fUji^lEl* 24f.i*2* F« . I , F.7 , 3 , Q T 1 9 0 0 

A F9,2/25x,A4,Fe,2,F9.2,Fll,2,F7,2,l6x,F8,l-F7,l,9x»2F9,5i OT 1910 

,?171. FpRHAjL_f/n^lX,A>«,Al,lA^A«,8X,F9,2^F11^2-Ei^,10l4l®*l4F743,F9.2, OT 1920 
A sro (S^stw A<i.r«.3>o 3 r < i . 9 . *7 . a . t . vfi. i wjxs ni totn 


OT 1910 
OT 1920 
OT 1950 
QT 1990 
OT 1950 
OT 19fc0 
OT 1970 
OT 1980 
OT 1990 
OT 2000 
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Du.nii.1 


. -AA.be. . 5 


* F9,2,2F9,5/23x,*«,ax,P'9.2,FU,8,F7,«,i6x,Fe,l,F7,J) 

92 1 i i OB M A T. . u Al^X.,-*<l.>.E1.2, F9, 2,Flltl,F 7 ,24^Jt4g ft ^J 
A 2F9,5/25X,AU,0X,F9,2,FU#2»F7,2,ifcX,Fe,l,F7,3) 

9221 FgBHA.L..t/njJJ^j.tA9«A2«Uil9#F6i2if9il»Fn»2,2F7,E»g5.»2»£J 
A F9.2,2F9.5/23x,A4,6x,F9,2,|rU.2,F7.2,16x,FS.l«F7,S) 

9222 F 0 B M A T ( / n^l.X.» 4 A-9 , A 2 , j X ,2,F11,2»2F74 

/k F9.2/23X, A4,0X,F9,2,Fn.2,r7,2,UX,F8|l#Ft,3,9X*2F9,5) 
End_ _ 
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•UgROUTlNE PLOTlT(NioELl*»««UM|L,88UMKO,XU,lPLOT,ITITUe>TReTRN, 


C THIS ROUTIniEi a subroutine QF the hAIn program pERPORH, PLOTS HALL 
C pressure PIPPgRgNTIiL aNP/OR CUHHUL^TlVg. HONDlHgNRIONAL PRE BBURE _ 
C LOSSES aOaINST c^HMULaTIVE CIRCUIT cENTeBLINB LeNOTM, ThIS PLOT 


C HOrE,,, NHEN plotting in $I units* CENtIHEtEr scales hill result; 


C*****«********A**************«*«****«*««*****«***«*****************< 


DIMENSION IX(*),IXN(S),IXNM(6),1Y(S),ITN(*),IVNM(S) 

-C»«« 4 _. — 

C.. RLOT axis labels arrays 

data IXN(1),1XN(2),IXN(I)»IXNC«)#IXNC5),IXN(S)/ «HCIRC,«HUIT , 


PT 

10 

PT 

-10 

* pT 

A at 

30 

ittO 

t Pi 

PT 

so 

Bl— 

-PD 

PT 

70 

L._JLT_ 


PT 

90 

L> al 

i.00 

* pT 

110 

• T 

PT 

130 

Pi 

-lilP 


DATA IXNH(i),lXNH(2)|IXNH(3},IXNM(4),XXNM(S)«lXNM(S)/ 4HCIRC* 

1 AHUIT .aHLirNO <iHTH f aMMgTi OMRST / 

DATA iYN(l)7iYN(2)*TVN(3),tYN(<0»lYN(5)»iYN(6)/ «HwALL,*H PRE, 

1 AHSSURfttHE fLf4MB/80,AH PT)/ 

DATA IYNM(|),IYNM(2)ilYNM(3),IYNMC«)#lYNMt5)»IVNH(S)/ 4MHALL* 

V AH pRE.4Ha8UB,4ME .aHMT I 




G« • • •] 

.Ct n 

C... 


n2 ■ N42 

F/>a T m l . 

XLEn ■ IS. 

YL6N i . 

PYLEN • vCEN.,8 
YLA B ■ YLEN Y-1 
XnEXT ■ IT. 

Y>^A X 1 11. 

YmARG • .5 


DEFINITION OF AXIS LABELS 

m 

"*"* Do so I ■ 1,6 

lYfl) ■ lYNfT) 

lX(I) > IXN(I) 

RQ CONTINUE 


C.... TE8T F or TY PE QF UniITS 
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£Lflm 


PaBC 


-tND-f LOT-TINS OF . A -SHALL 


C 1 1 • 

- ig_f -&j- 

cIILIaIM4n_k£ 

C.t.tTHE ORISIm 

■C # # » — 

1000 continue 

- -C-ALl P A .C TQP .(F A , CT 1 

Caul 8mooe(0,1|) 

Call PLOTto.o.»VMAi(.«i^ 

CALL PLoT(0«0fyHAR6f>3) 

C't . ^tT* ! • • , *^*°*^ ^ ^ ' ^ * * ' ^ * * ^ ^ * ^ *^ Ma *4J 

C n t«.PI.l£aNINiTloN qF SCLaIE . F-|>ft-Cl atUU 1 .E NStM, ^>N_ x.- AXIS 

-111 CALL aCiLE (SS UNE La X LE N iNaI) 

C.*«i. 

c » » 

t IfLO-T- 0.0 

.C*kk-*x la (0.t0 1 o»j>.i I.x t ?? .gA><LiNi 0 tOAa auHELtNn«aAUMrL(N2n 

Call 8YH9OLf0,5,VLABMl»ITlTUEAO,0,79) 

C all S CaL E(88Uwk,->>VL EN.N.i) 

^All AXiS(0t0#0,0,l3HpRE88URE l 088» n,aVi IN»90 , /sSUMkOCNI) » 
a sEUNKOfNzn 

Call L1NE(S8uN8l#88UNK0»N, 1,1,0) 

-_lg. t lP LQT ,eO , 1 1 so T O 25 o 0 

c T7i 

C....DEP|NlylON OP location Of ORIGIN PQp Nt h PLOt 

■ Call PLOTf0.0.»YMAX.»ii 

Call PLOT(0«0, yNArG,*3) 

Call PLOTfXNEXT.n,o«-l> 

Call symbol to. o,o.o,,ooi,3,o,o,*i) 



cII6uftftuCAT!vS*^Rl88UpE LOSS PLOTTED AGAINbx CIRCUIT LENGTH 
^ ^ ^ ^ 

~ 2000 continue 

CALL A xItf0.0.0.0.lY.*2tt.XLEN.Q.0.88UMeLtNl).ffgUNgLfN2n 

CALL 8YMBOLf,5,YLAB,.l,ltlTLE,0,O,7^) 

Call SCALEtPgLP.YLEN.N.) ) 

Call AX18(0.0,0.0,1Y,2«,PYLEn,90,,DELP(i41),DELP(n2)) 

-OALL LINEfSSUMEL.DELP.N. 1.1.01 

c!!)Ioi^^lON._g£..L&CA.T lON O F Oal 61jLgg«_N^7 PLOT 

2500 CONtINm E 


PT 510 

RT sao 

»T 510 
A- AI- -PX_ -5«0 
PT 550 

PlT — 5A0 

PT 570 

ftl — SAa 

PT 590 
PJ — *A0. 


PT PIO 
-Itl— 820 


PT 810 

XT -8A0 
PT 850 
-XT — 880- 


PT 870 

_pt._8ao 

pT 890 
*I ._ TOO 
oT 7l0 
_.Pi_. 120 
t PT 730 
,XT . TOO 
, PT, , 750 
^1- 780 
PT , 770 
- Rl 780 
PT 790 
. rT- 800 
PT 810 
PT 820 
PT 830 
.XT- 880 
pT 850 
. XL J80 
pT 870 
xi __B80 
PT 890 
„XI ..900 
pT 910 
PL... 92H 
PT 910 
_ pT 980 
pT 950 
-XT. 980 
PT 970 
-XT 980 
PT 990 
__PL-W00 
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CALL PlOT(0«0r"rMAX|«3) 
CALL PLOTtn.fl.YHARO.,11 
call PlOr(XNEXTfO«0|«S) 


return 



CMtiPARAMETER Conversion to si units 


3000 CONTINUE 

J-i£.T P 2,/l,26S9S.- 

XLEN ■ |S,;i,2*999 
VLgN ■ 10,4l,2i.OR9 
VLAS ■ yLEN*,t 

Jj-NE XT R I 7,*i,2fc 99^ 

Y«A* ■ 11. *1.26949 

yw ARg » ,5 «1,S6?99 

00 3SOO I I I«6 

ITU) a ITNMfH 

3S00 IX(I) ■ IXNM(I) 
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APPENDIX D 


INPUT AND OUTPUT FOR SAMPLE CASES 


Six wind tunnels were used, in addition to the test case (fig. 11) , as 
sample cases to establish the reliability and accuracy of the computer program 
analysis technique for th^ various types of duct components and wind tunnel 
circuits. Each case included here is titled with the appropriate wind tunnel 
name and its pages are numbered. The performance analyses are presented on 
the first two to three pages of each case. The summary characteristics tabu- 
lations and the plotted information were omitted. The annotated tabulations 
of the input data were included for reference. 

The results of the performance analyses are summarized in table 6. They 
are discussed and critiqued in the Results and Evaluation sections of this 
report. 
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•EMPTY* INDICATES OPTIONAL* NON-REQUIRED INPUT PARAMETER HAS BEEN OMITTFO CR PARAMETER MAY BE INTENDED AS ZERO. 
•ERROR* INDICATES MANDATORY INPUT PARAMETER HAS BEEN OMITTED. THIS MUST BE CORRECTED BEFORE COMPUTATION IS POSSIBLE 

•EXTRA* INCICATEr SUPERFLUOUS INPUT PARAhETEh HAS 6 e 6 n uMKESSAftlLY ll^CLUDEO ON INPu'T CAftD A» 4 t) HAV BE ftEMOVCO. 

•OPT*h* INDICATES OPTIONAL INPUT DATA HAS BEEN OMITTED AND THE PARAMETER MILL DEFAULT TO A PREDETERMINED VALUE. 
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TABLE 1.- NUMERIC INPUT CODE DEFINITIONS 


Code 

type 

Code 

value 

Description of code meaning 

Tunnel type 

1 

Closed test section, single-return tunnel 



2 

Closed test section, double-return tunnel 



3 

Closed test section, non-return tunnel 



4 

Open- throat, single-return tunnel 



5 

Open- throat, double- return tunnel 



6 

Open-throat, non-return tunnel 

Units of measure 

1 

International System of Units (SI) 


- 

2 

U. S . Customary Units 

Section 

shape 

1 

Circular cross section 



2 

Rectangular cross section 



3 

Flat oval cross section (ceiling and floor 




parallel with semicircular sidewalls) 

Section 

type 


(See table 4) 

Plot type 

iO.O 

No plots 



1.0 

Cummulatlve pressure losses vs circuit 




length 



2.0 

Wall pressure differential vs circuit 




length 



>2.0 

Cummulatlve pressure losses and wall 




pressure differential vs circuit length 
(on separate plots) 
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TABLE 2.- TUNNEL MASTER CONTROL INPUT DATA DESCRIPTIONS 











4-1 


AJ 




0 







""" 


“““ 




03 









4-1 


4h 




0 




o 



p 







1 


1 

1 

1 

1 

1 


1 







H 


u 


01 



p 



Vl 



•H 

1 


1 

1 

1 


1 

Vl 


Vl 



44 

4J 4J 

cO 

01 


CO 






o 



C 

1 


1 

1 

1 

1 


1 

O 


O 



O 

CO O 

01 

CO 




CO 





:z3 

P 











0 




01 01 

Vl 


Vi 

4H 

o 

V 

o 

H 


o 

p 

H 










6 





4-1 (0 

cO 

O 

p 

P 

o 

p < 


o 

o 

< 







Cl 





Vl 




























o 

S3 


01 








/-V 











1— 1 

o 



1— 1 


A 

O 


p 






01 

p 

T3 












o 


01 

CO 


rH 

•H 


4-1 

01 



01 


p 

p 

01 











Cl 



p 

> 


CO 

4-1 



p 





p 

S— ' 

4-1 











tH 

Cl 


cd 

o 


> 

CJ 


0 



p 

cO 


cO 


CO 











Pu 




*0 

o 

01 


•H 

o 


o 

a 


S3 

T3 

01 

0 


0 







1— 1 

CO 

cO 


CO 

4-1 

C 


CO 



4-1 


p 



P 

01 

0 

O 


o 


01 






AJ 




CO 

01 

4J 



rH 



p 

01 


CO 

4J 

cr p 


p 


Vl 





Cl 


CO 


T3 

rH 


CO 

4J 


01 

CO 



p 


4-1 

CO 

01 

4-1 


4j 


S3 




u 

rH 

Cl 



a 

4H 


CO 



p 





4J 

P 

Vl 

CO 


CO 


CO 




Cl 

P 

Cl 



Cl 


(0 4H 

01 


Q 



Vl 

O 


CO 

S3 


S3 


S3 


CO 

CO 




CO 

CO 

P 


6 

Vl 

01 


4-1 


a 

rH 


o 

4-1 



O 

CO 

60 


60 


Cl 




6 

U 

'w' 

Cl 


o 

M 

Vl 



01 


p 



a 

P 

p 

CO 


cO 


Vl 

a 

u 


9 



B 


CO 


4J 

O 

M 


01 

T3 



CO 


0 

CO 


4-1 


AJ 


p 

o 

01 


c 

Cl 

Cl 


Cl 

Vl 

CO 


CO 


p 

Q 



P 


a 

O 

60 

CO 


CO 



•H 

•H 



Cl 


CO 


M 

CO 

p 

Vi 

rH 


4-1 

a 


4J 



p 


0 



>w<' 


CJ 

4J 

CM 


Cl 

(0 

o 

Cl 


AJ iH 

rH 


cO 

S3 




p 

0 


a 

01 

P 





p 


•H 


(0 


o 

u 


CO 

S3 


rH 

a 


44 

44 


o 

O 


p 

^3 

P 


p 


Vl 

•H 

M 


CO 



AJ 


p: Cl 

60 4-1 

0 

M 

*X3 

o 

p 


o 

P 


a 


P 

CO 


CO 


Cl 

U 

a 


o 

Cl 

Cl 

CO 


5 rH 

0 

CO 

60 *H 

S3 


>-✓ 


p 

4J 



o 

a 

4-1 


4-1 


Ti 

a 

Cl 




Vi 

p 


O 

CO 


0 

o 

Cl 

Vl 



Cl 

CO 


p 

4J 

p 

o 


o 


p 

CO 

•a 


u 

o 

P 

P rH 

•TJ CO 

4-1 

c: 

CO 


a 

O 

S3 


> 

p 


o 


4J 

4J 


4-1 


CO 

0) 

•H 


Cl 

u 

CO 



4J 

a 

o 

4J 

44 

4-1 

O 



0 


p 

CO 

P 




01 

o 

P 



(0 

Cl 

cO 

a 

Cl 

a 

Cl 

•H 

Cl 

o 

CO 

U 

P 


S3 

o 


p 

P 

O 

S3 


S3 

Vi 

a 




3 


o 

iH 

O Cl 

u 

4-1 

Cl 


Cl 

CO 

4-1 


o 

p 




1 

O 


O 

0 

4J 


M 



& 

B 

•H 


•H Cl 


O 

Vl 

Vl 

M 44 

Cl »d 


p 

cO 



S>^ 

P 

01 

P 

AJ 

CO 


CO 





4J 

cO 

4-1 CO 

mh 

01 


01 

AJ 


01 01 


4-1 

O 


M 

4J 

4J 

4-1 

V4 

4J 

CO 



a 


Vi 

M 

4-1 

U 

AJ 

Cl V-/ 

o 

CO 

4H 

4J 

CO 

01 

CO 


a 



O 

P 

P 

o 

S3 

a 

Vl 

P 




CO 


O 

Cl 


Cl 



O 

01 

P 

60 

S3 


01 

Vl 


P 

Cl 

O 

01 

CO 

01 

01 

cO 


M 


Vi 

iH 


CO 

Cl 

CO Cl 

4-1 

u 


a 

P 


-VlM 


CO 

01 



O 

O 

CO 

CO 

CO 

P 

S3 


Cl 


v» 

Cl 

CO 


Cl 

T3 ^ 

CO 

P 

CO 



cq , o 





Vl 

P 

P 


01 


a 

U 


4J 


•H 

C 

4J 

•M 

CO 

4J O 

60 

Cl 

4J 

•H 

4J 

a 

tt»! C- 


4-1 

o 


Cl 

01 

Cl 

4J 

Vl 

4-1 

01 

01 


CO 



C 

CO 


CO a 

•H 

4-1 

TS 

T3 

CO 

o 

4J P 


CO 

p 



> 

> 

CO 

p 

CO 

AJ 

4-1 





a 


Cl 


Cl 

Cl 


•H 



rH 

i • • 


01 



o 



01 


01 


X 




<1 

H 

p 

H 


H 

P 





P 



H 



p 



H 


H 


p 


§> 

•H 



























U Cl 

CO 


Vi 

U 

Vi 

M 


Vi 




















9 a. 



Cl 

Cl 

Cl 

Cl 


Cl 




















& 

09 


00 

60 60 

00 


00 




















C 4-1 

2 


Cl 

Cl 

Cl 

Cl 


Cl 

rH 


rH 



rH 



p 



p 



P 


p 


p 

M 

a 


AJ 

M 

XJ 

AJ 


AJ 

CO 


CO 



CO 



CO 



CO 



CO 


CO 


CO 




c 

C 

a 

0 


0 

Cl 


01 



01 



01 



01 



01 


01 


01 


s 


M 

M 

M 

M 


M 

Pi 


Pi 



Pi 



Pi 



Pi3 



Pi 


Pi 


Pi 

CO 



















































/— N 





4J 























o 




o 

C 









• 














• 




• 






•H 




P 














p 




p 

6 

•o 


iH 

iH 



rH 

Cl 





rH 



’T3 



P 



V-X 


*T3 


><-✓ 

Cl 

Cl 


CO 

CO 


Cl 


CO 

p 


01 



CO 



01 



CO 



4-1 


01 


4J 

M 

M 


a 

a 

•H 

U 


0 



M 



0 



Vl 



S3 



rH 


Vl 


rH 

•H 

•H 


o 

o 

P 

•H 


O 

• 


•H 



O 



p 



O 



0 


P 


S3 




•H 

•H 

CO 

0 


•H 

B 


S3 



•H 



0 



P 



cO 


0 


CO 

cr 

cr 


4J 

M 

4H 

cr 


4-1 

o 


CT 



4-1 



XT 



4-1 



P 


U* 


P 

Cl 

Cl 


a 

p* 

Cl 

Cl 


P 

Cl 


Cl 



p 



01 



P 



Cl 


Cl 


01 

Pi 

& 


o 

o 

P 

PS 


o 

o 


Pi 



O 



Pi 



O 



Q 


Pi 


o 





m 


H 














P 














U 


H 






a 






W 










• 






M 












> 








CO 


O' 


|h 










p 






Si3 










» 




M 








p 






P 








iH Cl 


CO 









rH 
















Cl i-l 


w 



CO 

• 


• S 



Q 



p 






pcS 







s 





H H 

SC 

H SC 






p 













H 



CO 


5 

M 

O 

CO 

O CO 






p 













<0 



< 


P 




« 

rH 


rH 



S 



o 



O 



H 


H 





u 


H 


CO 


CO 

P 


S2 



P 



> 



P 



P 


H 


P 





























CO 





































m 


o 



m 



o 



m 



O 


m 


o 

s 









rH 


CM 



CM 



ro 



CO 







m 

U B 
CO D 

rH 


CM 

CO 

<r 

m 


VO 

1 

rH 


VO 



1 

p 



1 

vD 



1 

P 



1 

P 


1 

P 


1 

p 

O rH 









rH 


rH 



CM 



CM 



CO 



cn 




'O' 

w 

o 





























145 


Default (X)" indicates the input is optional and defaults to X if omitted. 

Geom. Dep." Indicates the input requirement is dependent on section geometry. 

Optional" indicates the input may be selected and Included as desired. 

Required" indicates the input must be non-zero and Included for all cases or the case will terminate 
due to input error. 








TABLE 3.- SECTION INPUT DATA DESCRIPTIONS 
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TABLE 3.- SECTION INPUT DATA DESCRIPTIONS - Concluded. 
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**Geom. Dep." indicates the input requirement is dependent on section geometry. 

"Optional” indicates the input may be selected and included as desired. 

"Required” indicates the input must be non-zero and included for all sections or section will be 
skipped and case terminated due to input error. 

"Sect. Dep." indicates the input requirement is dependent on section type. 









TABLE 4.- ADDITIONAL, SECTION-DEPENDENT INPUT REQUIREMENTS 


Section 

Additional 


Type description 

Type 

code 

input 

title(s) 

Requirement 

Single ducts: 




Test section, closed, constant 

01 

— 


area, empty 

Test section, closed, constant 

02 

S/AL 

Required 

area with model 


BLKGE 

Optional 



CD 

Required 



D EPS 

Optional 

Test section, closed. 

03 

KEXP 

Default 

diffusing, empty 
Test section, closed 
diffusing, with model 

04 

S/AL 

BLKGE 

Required 

Optional 



KEXP 

Default 



CD 

Required 



D EPS 

Optional 

Test section, open- throat. 

05 

— 


empty 

Test section, open- throat. 

06 

S/AL 

Required 

with model 


BLKGE 

Optional 



CD 

Required 



D EPS 

Optional 

Constant-area duct 

10 

— 


Contraction 

20 

— 


Corner, constant-area, turning 

30 

CHORD 

Required 

vanes only 


PHI 

Required 



KT 90 

Default 



RNREF 

Default 

Corner, constant-area, with 

32 

CHORD 

Required 

turning vanes and walls 


PHI 

Required 



KT 90 

Default 



RNREF 

Default 

Corner, constant-area, with 

33 

PHI 

Required 

walls and without turning 


KT 90 

Default 

vanes 




Corner, diffusing, with 

34 

CHORD 

Required 

turning vanes and walls 


PHI 

Required 



KT 90 

Default 



RNREF 

Default 

Diffuser 

40 

KEXP 

Default 

Exit kinetic energy from flow 

45 

— 


dump 

Sudden expansion 

46 

— 



Card 

column(s) 


36-40 

46-50 

61-65 

76-80 

56-60 

36-40 

46-50 

56-60 

61-65 

76-80 


36-40 

46-50 

61-65 

76-80 


41-45 

51-55 

56-60 

66-70 

41-45 

51-55 

56-60 

66-70 

51-55 

56-60 

41-45 

51-55 

56-60 

66-70 

56-60 









TABLE A.- ADDITIONAL, SECTION-DEPENDENT INPUT REQUIREMENTS - Continued 


Section 

Additional 

input 

title(s) 

Requirement?^ 

Card 

column (s) 

Type description 

Type 

code 

Flow straighteners, thin 

51 

L/DH 

Required 

36-40 

honeycomb 


PRSTY 

Required 

46-50 



RUFNESS 

Default 

66-70 

Flow straighteners, thick 

52 

L/DH 

Default 

36-40 

airfoils 


PRSTY 

Required 

46-50 

Perforated plate with sharp- 

53 

PRSTY 

Required 

46-50 

edged orifices 





Woven mesh screen 

54 

DMESH 

Required 

41-45 



PRSTY 

Required 

46-50 



KMESH 

Default 

56-60 

Internal structure (drag 

56 

ITEMS 

Default 

9-10 

Item(s)) at upstream end of 


S/AL 

Required 

36-40 

section 


BLKGE 

Optional 

46-50 



CD 

Required 

61-65 



D EPS 

Optional 

75-80 

Fixed, known local loss item 

57 

K 

Required 

61-65 

at upstream end of section 





Multiple ducts: 





Constant-area ducts 

61 

DUCTS 

Required 

7-8 

Contractions 

62 

DUCTS 

Required 

7-8 

Corners, constant-area; 

70 

DUCTS 

Required 

7-8 

turning vanes only 


CHORD 

Required 

41-45 



PHI 

Required 

51-55 



KT 90 

Default 

56-60 



RNREF 

Default 

66-70 

Corners, constant-area, with 

71 

DUCTS 

Required 

7-8 

turning vanes and only one 


CHORD 

Required 

41-45 

side-wall each 


PHI 

Required 

51-55 



KT 90 

Default 

56-60 



BNREF 

Default 

66-70 

Corners, constant-area, with 

72 

DUCTS 

Required 

7-8 

turning vanes and walls 


CHORD 

Required 

41-45 



PHI 

Required 

51-55 



KT 90 

Default 

56-60 



RNREF 

Default 

66-70 

Corners, constant-area, with 

73 

DUCTS 

Required 

7-8 

walls and without turning 


PHI 

Required 

51-55 

vanes 


KT 90 

Default 

56-60 

Corners, diffusing, with 

74 

DUCTS 

Required 

7-8 

turning vanes and only one 


CHORD 

Required 

41-45 

side-wall each 


PHI 

Required 

51-55 



KT 90 

Default 

56-60 

i 


RNREF 

Default 

66-''0 


149 






TABLE 4.- ADDITIONAL, SECTION-DEPENDENT INPUT REQUIREMENTS - Concluded. 


Section 

— 

Additional 

— 

Card 

column (s) 

Type description 

Type 

code 

input 

title(s) 

Requirement? 

Corners, diffusing, with 

75 

DUCTS 

Required 

7-8 

turning vanes and walls 


CHORD 

Required 

41-45 



PHI 

Required 

51-55 



KT 90 

Default 

56-60 



RNREF 

Default 

66-70 

Diffusers 

84 

DUCTS 

Required 

7-8 

Vaned diffuser 

85 

KEXP 

Default 

56-60 




Sudden expansion from multiple 

86 

DUCTS 

Required 

7-8 

ducts to single duct 
Sudden expansion from multiple 

87 

DUCTS 

Required 

7-8 

ducts to multiple ducts 
Fan, constant-area annular 

91 

DUCTS 

Default 

7-8 

duct(s) with motor-support 


ITEMS 

Default 

9-10 

strut (s) 


S/AL 

Required 

36-40 



DHUB 

Required 

41-45 



BLKGE 

Optional 

46-50 



CD 

Required 

61-65 



ETA 

Default 

71-75 



D EPS 

Optional 

75-80 

Fan contraction(s) to annular 

92 

DUCTS 

Default 

7-8 

duct(s) with motor-support 


ITEMS 

Default 

9-10 

strut (s) 


DHUB 

Required 

41-45 



BLKGE 

Optional 

46-50 

Fan diffuser (s) from annular 

94 

DUCTS 

Default 

7-8 

duct(s), each with tapering. 


DHUB 

Required 

41-45 

cone-shaped centerbody 


BLKGE 

Optional 

46-50 



KEXP 

Default 

56-60 

Internal structure (drag 

96 

DUCTS 

Required 

7-8 

item(s)) at upstream end of 


ITEMS 

Default 

9-10 

each duct 


S/AL 

Required 

36-40 



BLKGE 

Optional 

46-50 



CD 

Required 

61-65 



D EPS 

Optional 

75-80 

Fixed, known local loss item 

97 

DUCTS 

Required 

7-8 

at upstream end of each duct 


K 

Required 

61-65 


"Default** indicates the input is optional and has a default value if omitted 
(see table 3) . 


'’Optional** indicates the input may be selected and included as desired. 
’’Required** indicates the input must be non-zero and included for all sections 
of the specified type or the section will be skipped and the case not 
completed due to input error. 
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TABLE 5.- CASE TERMINATION TASK DESCRIPTIONS 


Card 

column (s) 

Input 

type 

Input 

value 

Task description 

3-4 

Blanks 

Blanks 

Case termination card identification 

6 

Integer 


Summary characteristics page(s): 



0 

Non-'print 



^0 

Print 

7-8 

Real 


Plotting of summary information as a function 




of distance through circuit; 



^0.0 

No plots 



1. 

Cummulative pressure loss 



2. 

Wall pressure differential 



>2. 

Cummulative pressure loss and wall pressure 




differential 

9-10 

Real 


Complete, annotated tabulation of input 




values ; 



0.0 

No print unless internally forced by 




omission of required inputs 



^0.0 i 

"Chosen” tabulation 

11-15 

Real 


Power-matching (optimizing velocity for a 




specified power level) : 



0.0 

No velocity optimization 



ito.o 

Velocity optimization 

16-20 

Real 


Return to beginning for evaluation of another 




case: 



0.0 

No return, program termination 



jtO.O 

Return 
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®The quoted energy ratios are the best available and best achieved for each facility. The energy ratios of some facilities have dropped over the years 
due to deterioration, leaks, soot build-up, etc. 
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Ratio of diffusing length to equivalent radius at upstream end, L/R 




Centerline 



Figure 2.- Contraction design criteria. 



NASA-Ames 40-by 80-Foot 
NASA'Ames 7-by 10-Foot 
Lockheed-Georgia low-speed 
(V/STOL test section) 



(b) Contraction design curves from reference 













(b) Turns without vanes. 
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Figure 4.- Straight-walled diffuser expansion loss parameter variation with equivalent cone angle; 

w _ w 1 /AR + 1\ 

KeXP “ %ef . 9 - 8 sin 6 
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Figure 5.- Mesh screen Reynolds number sensitivity factor as a function of Reynolds number. 
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Figure 6.- Vaned diffuser loss coefficient as a function of equivalent cone angle. 


Main Program (PERFORM) 



(a) Main program. 

Figure 7.- Basic functional flow chart. 
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(a) Main program - Concluded. 



Data-Checking Subroutine (DATACK) 


Section 
Card Checking 




Determination of 
Required General 
Inputs 




rT’TTiTrTi 


Typical 

Section 




Definition of Section 
Input Requirements 


Checking and Setting 
Error Codes 


r Return J 



Tabulation of 
Input Data 


Master 

Card Checking 


Determination 
of Required 
Master Inputs 

i 


Checking and Setting 
Error Codes 


C 


Return 


Assignment of Fixed 
Format Parameters 


Checking 
Requirements and 
Filling Arrays with 
Data or Message 


Print Tabulation of 
of Data and Messages 



(b) Data-checking subroutine. 
Figure 7.- Continued. 
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Local Speed Subroutine (SPEED) 



Local Reynolds Number and Friction Coefficient Subroutine 

(FRICTN) 




Reynolds 

Number 



(c) Local speed and Reynolds number/ friction coefficient subroutines. 

Figure 7.- Continued. 
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Input deck setup 
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Figure 9.- Sample headings for input data sheets. 
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(b) Wind tunnel section and overall performance information. 
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Figure 11.- Continued 
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(e) Summary information plots 
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(e) Summary information plots - Concluded. 
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(f) Annotated tabulation of input data. 
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(f) Annotated tabulation of input data - Continued. 















(f) Annotated tabulation of input data - Concluded. 
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(b) Lockheed-Georgla Low-Speed Wind Tunnel. 





(d) Hawker-Siddeley Aviation 15-Foot V/STOL Wind Tunnel at Hatfield. 
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(e) University of Washington 8- by 12-Foot Wind Tunnel. 






